Radioactive decays at limits of nuclear stability 
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The last decades brought an impressive progress in synthesizing and studying properties 
of nuclides located very far from the beta stability line. Among the most fundamen- 
tal properties of such exotic nuclides, usually established first, is the half-life, possible 
radioactive decay modes, and their relative probabilities. When approaching limits of 
nuclear stability, new decay modes set in. First, beta decays become accompanied by 
emission of nucleons from highly excited states of daughter nuclei. Second, when the nu- 
cleon separation energy becomes negative, nucleons start to be emitted from the ground 
state. Here, we present a review of the decay modes occurring close to the limits of 
stability. The experimental methods used to produce, identify and detect new species 
and their radiation are discussed. The current theoretical understanding of these decay 
processes is overviewed. The theoretical description of the most recently discovered and 
most complex radioactive process — the two-proton radioactivity — is discussed in more 
detail. 



PACS numbers: 23.50.-l-z, 23.60.- 


fe, 23.90. -l-w, 23.40. He 






CONTENTS 






2. ISOL 








3. Future facilities 


I. INTRODUCTION 


m 


c. 


Detection 


A. Radioactivity and Nuclides 


m 




1. Charged particles 


B. Brief history 


m 




2. Neutrons 


C. Links and connections 


m 




3. Signal processing 


1. Nuclear structure 


m 






2. Nuclear astrophysics 


m 


IV. BETA DELAYED PARTICLE EMISSION 


3. Open quantum systems 


m 


A. 


Beta decay, general observations 


D. Outline 


m 




1. The beta strength 






2. Delayed particles 


II. LIMITS OF STABILITY 


m 


B. 


/3+ delayed emission of one particle 








1 . Occurrence of particle emission 


III. EXPERIMENTAL TECHNIQUES 






2. Fermi decays 


A. Production 






3. Gamow- Teller decays 


1 . Fusion-evaporation 


nni 




4. Selected spectroscopic tools 


2. Multi-nucleon transfer 




C. 


f}~ delayed emission of one particle 


3. Fragmentation 


m 




1. Occurrence of particle emission 


4. Spallation 


m 




2. Decays in different mass regions 


5. Fission 






3. Selected spectroscopic tools 


B. Separation 


[la 


D. 


Beta delayed emission of several light particles 


1. In-FUght 






1. The case of A = 12 






E. 


Beta delayed fission 



US 

CS] 
M 
El 
EE] 



' |pfutzner@fuw.edu.pl| 



V. SINGLE-PROTON RADIOACTIVITY 
A. Introduction 

1 . Fundamentals 



2 



2. Probability of proton emission 

3. Spectroscopic factor Sp 

4. Models of proton emission 
Seniority s < 2 proton emitters 

1. Odd-mass, s = 1 proton emitters 

2. Even-mass, s = 2 proton emitters 
Seniority s > 2 proton emitters 
Outlook 



EH 



VI. ALPHA DECAYS 

VII. TWO-PROTON RADIOACTIVITY 

A. Introduction 

1. Two-proton correlations 

2. Historical note 

B. Experimental results 

1. Democratic decays 

2. Two-proton emission from excited states 

3. Two-proton radioactivity 

C. Simplified theoretical models 

1. Direct decay model 

2. Simultaneous vs. sequential decay 

3. Diproton model 

D. Three-body model of 2p radioactivity 

1. Lifetimes 

2. Spatial correlations 

3. Momentum correlations and nuclear structure 

4. Long-range character of Coulomb interaction in 
the three-body continuum 

E. Two-proton radioactivity and many-body nuclear 
structure 

1. Spectroscopic information in R-matrix 
approaches 

2. Three-body model plus RMF amplitudes 

3. Shell model approaches 

4. Microscopic cluster models 

F. Three and four proton emission 

VIII. EMISSION OF NEUTRONS 
IX. CONCLUSIONS 
Acknowledgments 
References 



Ho] 

M 
m] 



113 



EH 
ED 



I. INTRODUCTION 

We consider the atomic nucleus as a quantum ob- 
ject composed of A nucleons (mass number): Z protons 
(atomic number) and N neutrons, held together mainly 
by strong nuclear forces. A neutral atom with the spec- 
ified numbers A and Z is called a nuclide. When using 
this term, however, we focus on the nuclear component 
of the atom. Such a system is stable only for certain 
combinations of numbers Z and TV. Presently, 256 stable 
nuclides are known. Systems different from stable con- 
figurations undergo spontaneous, radioactive decays until 
the stability is reached. A nucleus of such an unstable 
nuclide is considered as a well defined object if its half- 
life is much longer than 10~^^ s which is a characteristic 
timescale for processes governed by strong interaction. 
These nuclides are bound by nuclear forces and/or by 
Coulomb and centrifugal barriers. The number of un- 
stable nuclides synthesized in laboratories is constantly 
growing, and up to now more than 3000 were identified. 
In this review, we concentrate on radioactive processes 
observed for nuclides located at the limits of the nuclear 
chart. The emphasis is given on new decay processes 
and features of classical decay modes which do not take 
place among nuclides close to stability. We will refrain, 
however, from discussing the heavy frontier of the nu- 
clear chart. The quest for the superheavy elements was 



reviewed by 
recently by 


Hofmann and Miinzenberg 


(20001 and more 


Oganessian 


(2007 


1 and by 


Hofmann 


( 2009a 1. 



A. Radioactivity and Nuclides 

The notion of radioactivity is useful in making distinc- 
tion between emission of rays or particles by a highly 
unstable system (for example undergoing a nuclear reac- 
tion) from radiation emitted spontaneously by a system 
whose nuclear and atomic degrees of freedom are close 
to equilibrium. Such distinction, however, has to be ar- 
bitrary and usually a characteristic time scale is used as 
a criterion. Throughout this review we adopt the follow- 
ing definition. Radioactivity is a process of emission of 
particles by an atomic nucleus which occurs with char- 
acteristic time (half-life) much longer than the K-shell 
vacancy half-life in a carbon atom, which amounts to 
about 2 X 10 



-14 



s (Bambynek et al. 



19721 



A relativis- 
s a distance of a 



tic particle travels in the time of 10 
few micrometers which is close to the measurement limit 
in a nuclear emulsion. In addition, this value coincides 
with a decay width, defined as F = \n2h/Ti/2, of about 
0.03 eV which is roughly the thermal energy at room 
temperature. Thus, nuclear processes much slower than 
filling the K-vacancy, whose duration in principle can be 
measured directly, and with the width much smaller than 
the thermal energy at room temperature will be called ra- 
dioactive. This definition applies to both nuclear ground 
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FIG. 1 (Color online) The chart of nuclei. The stable nuclides and the radioactive ones which were experimentally identified 
are shown by black and light/yellow, respectively. The nuclides predicted to have positive nucleon separation energy according 
to the FRDM mass model ( [Moller et al.\fE)97\ , but not yet observed, are shown in dark/green. The lines indicate position of 
magic numbers corresponding to the closed neutron and proton shells (the numbers smaller than 20 are not shown). The insets 
show the location on the chart of the decay products of the parent nucleus which is indicated by a dark square. The observed 
decay channels of the proton-rich and the neutron-rich nuclei are shown on the left and on the right inset, respectively. 



states and to long-lived excited nuclear states (isomers). 

The definition of a nuclide relates to the definition of 
radioactivity. A nuclide is a neutral atom, specified by 
the numbers A and Z of its nucleus, which is either stable 
or lives long enough to be classified as radioactive. We 
say that a nuclide does not "exist" if its nucleus decays 
too fast to be called radioactive. All existing nuclides are 
represented on a chart of nuclides spanned by the atomic 
number Z and neutron number N (Figure [ij . In the last 
three decades their number was growing almost steadily 



from about 2200 in 1981 to about 3000 in 2006 ( [Pfennig 
et al. 2008| ), giving an average of about 30 new nuclides 
identified per year. Due to vigorous growth of nuclear 



facilities (Sec. III.B ) this trend is expected to continue in 
next decades. 

The domain of processes occurring on a time scale 
shorter than radioactivity is referred to as the resonant 
regime. The resonant phenomena are characterized by 
features having a directly measurable width in the en- 
ergy spectra. Typically, widths of order meV is taken 
as the lower end of the resonant regime. Thus, the dis- 
tinction between radioactive and resonant phenomena is 
that the former have a characteristic time which can be 
measured directly while the latter have a characteristic 
energy width which can be measured directly. 



Characteristic time scales for different radioactive de- 
cays are illustrated schematically in Figure [2j 



B. Brief history 



The discovery of radioactivity by Becquerel ( 1896 ) and 



the subsequent discovery of new radioactive elements, 



polonium (Curie and Curie 18981 and radium (Curie 



et al. 1898 1 initiated a scientific breakthrough into the 
world of subatomic structure of matter. Through early 
works on a- and /3-rays, distinguished and named by 



Rutherford ( 1899 ) and on 7-radiation, discovered by Vil- 
lard (19001, a manifold of nuclear phenomena started 
to emerge. The nature of the new rays was soon clar- 



ified ( 


Rutherford and Geiger 


1908 


1, ( 


Kaufmann 


19021, 


(Rutherford and da C. Andrade 




19141. Finally, two 



other milestones: the discovery of the atomic nucleus 



(Rutherford 1911 1 and of the neutron (Chadwick 19321 
gave birth to a new discipline — nuclear physics. A bril- 
liant and detailed account of these early years was given 



by Pais ( 1986 ) 



The full understanding of decay processes required 
quantum mechanics, as an illustrious explanation of a- 



decay by Gamow (19281 and Gurney and Condon (19281 
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FIG. 2 Characteristic time scales and decay widths of radioactive decays. Ranges of application of selected experimental 
techniques are sketched. 



demonstrated. Soon afterwards, early ideas of quantum 
field theory led Fermi to formulate the first theory of /3- 
decay (Fermi 19341. A novel type of radioactivity 
/S"*" decay 



It 



- was discovered by Curie and Joliot ( 1934 1 
was also the first instant of a radioactive nuclide synthe- 
sized in a laboratory in contrast to all previously known 
radioactive substances of natural origin. Wick (19341, 
and independently Bethe and Peierls ( 1934[ ) , realized that 



a process of orbital electron capture (EC) by an atomic 
nucleus is possible. Thus, the EC decay was the first 
type of radioactivity predicted theoretically. It was con- 



firmed by an observation of K-capture by Alvarez ( 1937 1 



Spontaneous fission, discovered by |Flerov and Petrzhak| 
(1940), completed the list of "classical" radioactive de- 
cay modes. Ever since, they have played a crucial role 
in learning about nuclear properties, in tagging reaction 
channels, and in identifying new nuclides. 

The beginning of the modern era of radioactivity can 



be marked by the work of Goldansky ( 1960 ) who in a 



systematical study considered properties of very neutron 
deficient nuclei and discussed possible decay modes like 
/3-delayed proton emission and proton radioactivity. He 
was the first to point out the possibility of the two-proton 
radioactivity and to describe its key features. The first 
emission of a proton following fi decay was reported by 
Karnaukhov et al. (19631. The first emission from an 



identified precursor was observed from an excited state 
of ^^Al populated in the /S"*" decay of ^^Si (Barton et al. 
19631. Such a variant of /3 decay — the /3-delayed emis- 
sion of particles from states of a daughter nucleus — 
become a field of study in its own and includes vari- 
ous decay channels with emission of protons, neutrons, 
a particles, deuterons, and tritons, see Figurtfl] This IS 
the subject of Sec. |IV[ The first direct emission of a pro- 
ton from a nuclear state was observed to proceed from an 
isomeric state in ^'^Co (Jackson et al. 1970 1. The ground- 



state prot on radioactivity was o bserved for the f irst time 
in i^^Lu (iHofmann et all |1982|) and in ^-^^Tm (iKlepper 



the et al. 19821. Since then, almost fifty proton emitters, 
including emission from isomeric states, were identified. 
This field is covered in Sec.|V] The prediction of the two- 
proton radioactivity {2p) had to wait much longer for 
experimental confirmation. Such a decay mode was ob- 
served first for ^^Fe by Pfiitzner et al. (2002 1 and in an in- 



dependent experiment by Giovinazzo et al. (2002). This 



freshly opened field of nuclear spectroscopy is covered in 
Sec.lflll 

For completeness, one should mention the cluster ra- 
dioactivity, in which a nuclear fragment heavier than 
the a-particle but lighter than fission fragments, is emit- 



ted. T his decay mode was discovered by ,Rose and Jones 
( 1984 ) who identified ^'^C ions emitted by ^^^Ra. Later, 



many similar decay channels were observed ( Bonetti and 



^Guglielmetti 


2007 


Poenaru et al. 


2002 


however, the dominant decay moc 


e is a 



ter emission being less probable by a factor of at least 
10^. In addition, this rare decay mode can be estab- 
lished only in long-living nuclei, not far from the nuclear 
stability. Thus, it is beyond the scope of this paper. 



C. Links and connections 

In general, information on nuclei come from different 
and complementary experimental approaches. Radioac- 
tive decay always was one of them and still remains a ma- 
jor tool in nuclear physics as well as in numerous branches 
of physics where nuclear degrees of freedom are relevant. 
An important attribute of the decay is the characteristic 
half-life. The very delay between the moment of produc- 
tion and the decay event, offers a filtering possibility, thus 
helping to increase the signal-to-background ratio. The 
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values of the half-life provide a way to identify processes 
and offer tests to nuclear models. Another important 
facet of spontaneous decays are selection rules, which re- 
flect conservation laws obeyed by underlying interactions. 
Effectively, they provide different filtering mechanisms, 
essential for the planning of measurements and for inter- 
pretation of their results. 

The quest for superheavy elements may exemplify the 



decays of 75 atoms (Miernik et al. 2007c ). The first infor- 
mation on excited states in '^'^Ni was deduced from a few 



tens of counts in a singles gamma spectrum ( Grzywacz 



et al. 1998). Here, the very clean selection was pro- 



vided by a delayed decay of a microsecond isomer. Such 
first information, in turn, becomes essential in later, more 
advanced experiments, allowing for coincidence measure- 
ments or for selection of events in high-background con- 



role of radioactive decays (Hofmann and Miinzenberg ditions, for example, by decay tagging (Jenkins et al. 



20001. The emission of a particles by synthesized heavy 2000 Seweryniak et al. 19971 



nuclides provides the very proof of their existence and a 
clean way for their identification. Moreover, the a decay 
is a powerful marker helping to establish the first chem- 
ical properties of a few-atom sample of a new element. 



Recently, the insights gained on nuclear structure near 
the proton drip-line by means of radioactive decay studies 



were reviewed by Blank and Borge ( 2008 1 



like in the case of ^^^Cn (Eichler et al. 20071 



Radioactivity plays a significant role in fundamental 
research, as can be illustrated by historical examples of 



the neutrino hypothesis by Pauli ( 1930 1 and the discovery 
of parity non-conservation in the /3 decay of ^°Co (Wu 
et al. 1957). A recent review of tests of the standard 



electroweak model by the beta decay (Severijns et al. 



2006 ) gives a comprehensive summary of this field. 



However, since the present paper is limited to radioac- 
tive studies at the limits of stability, below we mention 
briefly a few research areas where they are particularly 
important. 



1. Nuclear structure 

With increasing scope of studies on the chart of nuclei 
new patterns emerge and novel features are predicted. 
An observed anomaly in masses of neutron-rich sodium 



2. Nuclear astrophysics 

Several subjects of modern astrophysics are related to 
properties of atomic nuclei. They include questions on 
the origin of the elements, physics of compact objects 
like neutron staxs or white dwarfs, studies of stellar ex- 
plosions like supernovae, x-ray bursts, and many others. 
Although the most frequently requested nuclear data in 
astrophysics are reaction cross sections for some key pro- 
cesses, the radioactive decays provide indispensable in- 
formation for many applications. Nuclei very far from 
the beta stability play a particularly important role in 
nucleosynthesis of elements heavier than iron, as pointed 



out in the seminal paper of Burbidge et al. (19571. The 



isotopes (Thibault et al. 1975 1 is an early example which 



led to the notion of an "island of inversion" in vicinity of 
32]\/[g resulting from changes in sequence of single-paxticle 



rapid neutron capture (r-process), responsible for about 
half of the abundance of elements above iron, including 
all of uranium and thorium, passes through regions of 
very neutron-rich nuclei, mostly far beyond the reach 
of present experiments ( [Kratz et al. 20071. In turn, a 
number of neutron-deficient isotopes was produced in the 
rapid proton capture (rp - process) which involves very 



orbitals (Nummela et al. 2001 Yordanov et al. 20071. neutron-deficient areas on chart of nuclei (Schatz et al. 



Recently, a similar phenomenon was claimed to occur 1998). Both processes are expected to occur in explosive 



around ^^Ti (Flanagan et al. 



2009 



Tarasov et al. 



20091. 



It was realized that the classical shell gaps — the corner- 
stones of microscopic description of nuclei — do migrate 
in areas distant from stability. A detailed survey of the 



present situation is given by Sorlin and Porquet ( 2008 1 



A new feature of quenching of the shell-gaps due to influ- 
ence of continuum states was predicted close to the neu- 
tron drip line (Dobaczewski et al. 1994). A comprehen- 



sive review of novel aspects of nuclear structure emerging 
far from the beta stability is given by [Dobaczewski et al.\ 
(|2007l). 



The experimental reach for most exotic nuclides suf- 
fers from very low production yields. Radioactive de- 
cays allow not merely to identify a system of interest but 
often offer the only practical source of structural infor- 
mation. The first half-life determination for the doubly 



magic ^^Ni was done with 11 atoms ("Hosmer et al. 2005 1 



The first insight into the structure of ^Fe resulted from 



stellar conditions, details of which are still under debate. 
Theoretical reconstruction of these processes is addition- 
ally hindered by the lack of relevant nuclear data. To 
those of special interest belong masses, half-lives, branch- 
ing ratios, and beta-delayed particle emission probabili- 
ties. In principle, these values can be determined from 
radioactive decay studies. Presently, however, most of 
the nuclei important in this context are either difficult 
or impossible to synthesize in reasonable quantities. The 
theoretical models or empirical parameterizations have 
to be relied on, instead. Thus, the motivation to extend 
experimental decay studies to the limits of nuclear sta- 
bility is twofold. They deliver data directly needed by 
astrophysical calculations, and they help to test and im- 
prove models of nuclear structure, increasing reliability 
of theoretical extrapolations. 

A thorough and detailed discussion of the present sta- 
tus of the interplay between nuclear structure and as- 
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trophysics, with the emphasis on stellar evolution and the least known and its understanding is still in a status 



nucleosynthesis, can be found in Langanke and Martinez 



Pinedo (20031 and in Grawe et al. (20071. Both these re 



views reveal the necessity of exploring exotic nuclei and 
their decay properties in order to fully understand astro- 
physical aspects of our Universe. 



of development. In Section VIII the prospects of neutron 
radioactivity are examined. The main conclusions of the 
paper are shortly summarized in the final Section |IX[ 

Throughout this work we use the system of units in 
which h = c = 1. 



3. Open quantum systems 



Several fields of modern physics face the goal of de- 
scribing quantum many-body systems which are not iso- 
lated from its quantal environments. Examples of such 
open quantum systems include quantum dots, droplets 
of neutral atoms, microwave cavities, or weakly bound 



nuclei very far from the beta stability ( Okolowicz et al. 



2003 1 . Nuclear physics at the limits of stability appears 



as a particularly promising testing ground of new con- 
cepts. When the nuclear binding energy decreases, the 
conceptual separation of well localized, bound states from 
the continuum scattering states becomes artificial and, 
in fact, is hampering the correct description of various 
features, like neutron halos, Thomas-Ehrman shifts, or 



et al. 2010a I 



clustering phenomena (Dobaczewski et al. 2007 Michel 



The quest to formulate a unified description of nuclear 
structure and nuclear reactions resulted in extensions of 
nuclear shell model, like the Shell Model Embedded in 



the Continuum (Bennaceur et al. 1998 Okolowicz et al 



20031, or more recently the Gamow Shell Model (Michel 



et al. 2002 2009). New spectroscopic data on very exotic 



nuclei, and on complex decay modes, like two-proton ra- 
dioactivity, or beta-delayed multi-particle emission would 
stimulate further developments in this field. Thus, ra- 
dioactive decays at the limits of nuclear stability may be 
instrumental in improving our fundamental understand- 
ing of many-body quantum systems. 



D. Outline 

In Section |ll] we discuss the limits of nuclear stabil- 
ity using the concept of the drip-line based on nucleon 
separation energy. The experimental situation in access- 
ing both the proton and the neutron drip-line is briefly 
presented. The experimental techniques pertaining to ra- 
dioactivity studies far form stability are reviewed in Sec- 
tion |III[ Various reactions used to produce exotic nuclides 
and the main methods of their separation are shortly de- 
scribed. Finally, selected detection techniques of special 
importance for measurements of radioactive decays are 



presented. The following Sections IV VII are devoted to 
the main radioactive decay modes at the limits of stabil- 
ity: /3-delayed particle emission, proton radioactivity, a 
emission, and two-proton radioactivity. The latter decay 
mode is treated in considerably more detailed way, as it is 
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II. LIMITS OF STABILITY 

The limits of the nuclear world are determined by the 
nuclear binding energies. The limits relevant to this re- 
view are often characterized by the drip lines which sep- 
arate bound systems from the unbound ones. Although 
different definitions can be encountered in the literature, 
we adhere to the simplest and most common one which 
is based on the single-nucleon separation energy. The 
proton- and the neutron separation energy of a nuclide 
with numbers N and Z are given by: 

Sp{N,Z) ^ B{N,Z) - B{N,Z -I) (1) 
SniN,Z)^B{N,Z)-B{N-l,Z). (2) 

The B{N, Z) is the binding energy of the nuclide related 
to its mass M{N, Z) : 



M{N, Z) = ZAlH + Nm„ - B{N, Z), 



(3) 



where Mh and m„ are masses of the hydrogen atom and 
the neutron, respectively. 

When we move along the line of isotopes with the 
given atomic number Z, starting from stability towards 
neutron-deficient nuclides, the proton separation energy 
Sp decreases and at certain location it becomes negative. 
The proton drip-line is defined as the border between the 
last proton-bound isotope and the first one with the neg- 
ative value of the Sp. The typical situation, according to 



the predictions of a particular mass model ( MoUer et at. 
1997[ ) for the isotopes of iron and cobalt is presented in 



9 and 26Fe2o, while 
27>^022 and 27C023. 



Figure [3] It follows from this model that the proton drip- 
line for iron should lie between ||Fei9 and 
in case of cobalt it is located between 27C0: 
Generally, the proton drip-line for odd Z isotopes is closer 
to stability than in case of the neighboring even-^ which 
results from the proton pairing energy. 

In the fully analogous way, the neutron drip-line for a 
given neutron number N is defined as a border between 
the last neutron bound isotone, when counting from sta- 
bility, and the first one for which the neutron separation 
energy Sn is negative. The predicted separation energies 
for the = 26 and N = 27 isotones are shown in Fig- 
ure |4] Thus, for iV = 26 the neutron drip-line is expected 
to lie between g^F26 and foNe26. Similarly to the proton 
case, the neutron drip-line for the odd-A^ is closer to sta- 
bility than for the neighboring even-A^ which reflects the 
neutron pairing energy. 

The drip lines as defined above are very useful in iden- 
tifying and discussing limits of stability, but to some ex- 
tend they are arbitrary and they do not provide the un- 
ambiguous demarkation of nuclear stability. This can be 
seen by inspecting the two-nucleon separation energies: 

S2p{N,Z)^B{N,Z)-B{N,Z-2) (4) 
S2n{N,Z)^BiN,Z)-BiN-2,Z). (5) 

Due to the pairing interaction, in case of an even num- 
ber of nucleons the two-nucleon separation energy can 




19 20 21 

Neutron number, N 

FIG. 3 (Color online) The proton- and two-proton separation 
energies of iron and cobalt isotopes as predicted by the FRDM 
mass model (iMoller et aLl|1997l. 



be smaller than the single-nucleon value. For the cases 
discussed above, it is illustrated in Figures |3] and |4] Al- 
though the proton separation energy in |gFe2o is posi- 
tive, this nuclide is expected to be slightly two-proton 
unbound. Similar situation is observed in the A^ = 26 
isotones — the two-neutron instability develops first, be- 
fore the neutron drip- line is reached. The additional 
complication, which is essential on the neutron-deficient 
side, comes from the fact that the exact position of a 
nuclide with respect to the drip-line cannot determine 
alone its dominant decay mode. This is caused by the 
Coulomb and centrifugal barriers which hamper emis- 
sion of nucleons. Only when the nucleon penetration 
probability through the barrier, depending on the energy 
and the angular momentum of the initial state, is large 
enough, the particle radioactivity can compete with /3 
decay. Thus, although the mentioned 26Fe2o may be two- 
proton unbound, it is known to decay by /?"'" transition. 
In turn, ||Fei9 is sufficiently two-proton unbound to de- 
cay predominantly by the 2p radioactivity, although the 
/?"'" channel has a substantial branching (Miernik et al. 



20091. The exact position of this nucleus with respect 



to the proton drip-line turns out to be irrelevant for its 
radioactive decay. In case of odd-Z nuclides, the proton 
radioactivity can win the competition with /3 decay only 
when its proton separation energy is sufficiently negative. 
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10 11 12 

Proton number, Z 



FIG. 4 (Color online) The neutron- and two-neutron separa- 
tion energies for the A'' = 26 and N = 27 isotones, as predicted 
by the FROM mass model ||Moller et aE||1997|. 



Thus, the observation of the proton radioactivity proves 
that the nuclide is located beyond the proton drip-line 
but the exact position of this line cannot be determined 
from decay data alone. This can be achieved only by 
precise mass measurements of nuclides in the region of 
interest. 

We note, that on the neutron-deficient edge of the 
chart of nuclides above tungsten the dominant decay 
mode is the a emission which happens to proceed faster 
than /3 decay. Thus, beyond the proton drip-line in this 
region, the proton radioactivity competes actually with 
a decay. 

At the neutron drip-line the situation is different be- 
cause the unbound neutrons are not affected by the 
Coulomb barrier. The influence of the centrifugal poten- 
tial alone is much weaker as it decreases with radius ef- 
fectively as in contrast to the l/r dependence of the 
Coulomb potential. In consequence, the effect of the cen- 
trifugal barrier is expected to be observable only in rare 
cases of very low decay energies and large angular mo- 
mentum. The resulting possible neutron and two-neutron 



radioactivity is examined in more detail in Sec. VIII For 
practical purposes, any system with negative neutron- or 
two-neutron separation energy can be expected to live 
too short to be qualified as radioactive. Therefore, the 
limits of stability on the neutron-rich side could in prin- 
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FIG. 5 (Color online) The difference between the neutron 
number of the lightest experimentally observed isotope for a 
given atomic number Z and the corresponding prediction for 
the last isotope befo re the proton drip-line according to the 
FRDM m ass model (|Moller et fflZ.[ [T997)(line) and the HFB- 
17 model (Goriely et al. 2009[ ) (circles). The results for the 
even Z and the odd Z are shown in the bottom and in the 
top, respectively. The experimental values were taken from 



Magill et al. { 2009 1 with corrections contained in Baumann 

et aumori. 



ciple be established rather precisely by inspecting which 
is the lightest isotone still undergoing radioactive decay. 
The problem, however, is that it is very difficult to reach 
experimentally the neutron drip-line for N > 28. 

The general picture is presented in Figure [l] where all 
nuclides identified experimentally until now are superim- 
posed on the plot of all nuclides predicted to have pos- 
itive proton and neutron separation energy by the Fi- 



nite Range Dro plet Model (FRDM) developed by |M511er 
et aZ.| ( |1995| and |Moller et aZ.| ( |1997| . This model is a sue 



cessful represent ant of a class of macroscopic-microscopic 
mass formulae, combining the macroscopic liquid-drop 
parametrization with the microscopic shell and pairing 
corrections. The prediction reveals a few characteris- 
tic features, like the even-odd staggering for the neutron 
deficient isotopes and for the proton deficient isotones, 
or the strong influence of the N — 82 and N = 126 
shells on the neutron-rich side of the chart. To illustrate 
theoretical uncertainties we compare predictions of the 



FRDM model with results of the HFB-17 model (Goriely 



et al. 



2009 1 which represents a class of fully microscopic 
approaches based on Hartree-Fock-Bogoliubov formal- 
ism and Skyrme forces (Goriely et al. 2010). Figure [s] 
presents the current experimental situation and the com- 
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are positive and they increase rapidly with the increas- 
ing neutron number. This reflects the fact that except 
for the light nuclei, the neutron drip-line is far from the 
body of presently observed nuclides, as is dramatically 
evident also in Figure [T] In fact, the drip-line has been 
determined experimentally and unambiguously for even 
N only up to = 20 and for odd N up to TV = 27 



(Thoennessen 20041 



The significant expansion of the body of observed nu- 
clides, especially on the neutron-rich side of the chart, 
is expected only when the next generation of radioactive 
beam facilities will come into operation, see Sec. |III| 

The current status of the knowledge on atomic masses 
and of the global mass models can be found in |Lunney| 



et al. ( 2003 1 . The detailed discussion of the present ex- 



perimental knowledge of the limits of nuclear stability 



was presented by Thoennessen ( 2004 1 



FIG. 6 (Color online) The difference between the proton num- 
ber of the lightest observed isotone for a given neutron number 
and the prediction of the last stable isotone before the neu- 
tron drip-line according to the two theoretical mass models. 
The results for the even A'^ and the odd A'^ are shown in the 
bottom and in the top, respectively. The plot details are the 
same as in Figure [5] 



parison of the two models. It shows differences between 
the neutron number of the lightest observed isotope and 
the predicted values for the last proton-stable isotope be- 
fore the proton drip-line. We see that both predictions 
agree well with each other — they follow the same pat- 
tern and they differ by a few units at most. The negative 
values indicate those observed nuclides which are located 
beyond the predicted proton-drip line, in most cases they 
are proton emitters. The large group of such nuclides, 
seen for odd-Z values between 50 and 90 illustrates the 
strong impact of the Coulomb barrier on the heavy nu- 
clei. On the other hand for almost all even-Z elements, 
there are predicted bound isotopes which remain to be 
observed. A distinguished peak of positive values for the 
Z > 90 results from experimental difficulties to produce 
proton-rich nuclei in this region. 

The analogous information for the neutron-rich side 
is given in Figure |6] where the differences between the 
proton number of the lightest observed isotone and the 
predicted values for the last neutron-stable isotone before 
the neutron drip-line are plotted. Again, both models are 
consistent with each other. The largest difference is seen 
above the A^ = 80 where the FRDM model seems to ex- 
hibit large variations due to the neutron shell closure — 
an effect not pronounced in the HFB-17 model. In con- 
trast to the proton-rich side, however, almost all values 
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III. EXPERIMENTAL TECHNIQUES 



1. Fusion-evaporation 



Experimental studies of nuclei at the limits of stabil- 
ity belong to the front-line of physical research. A view 
on experimental techniques, given in this section, pro- 
vides a general perspective on the advanced methods of 
present-day low-energy nuclear physics. First, reactions 
used to produce radioactive nuclides will be mentioned 
followed by a short description of the main methods of 
their extraction and separation. Then, selected aspects 
of modern detection systems will be reviewed with an 
emphasis on recording manifestations of radioactivity. 



A. Production 

The methods of production of nuclides far from beta 
stability are almost exclusively based on nuclear reac- 
tions involving stable nuclides or their ions. In a simpli- 
fied view, a new nucleus is formed either by fusion of two 
other nuclei (projectile and target), by exchange of nucle- 
ons between the projectile and the target nuclei (trans- 
fer) , or in reaction leading to removal of nucleons either 
from the target or from the projectile nucleus (fragmen- 
tation, spallation, fission). 

In principle, a radioactive nucleus produced in one of 
these reaction, having sufficiently long half-life, can be 
used as a projectile to initiate a secondary reaction in 
which nuclei even further from stability are formed. This 
is the idea of radioactive beams which has been driv- 



ing many experimental developments (JPG 2011, Tani 



hata 2008). It is anticipated that reactions induced by 



radioactive beams will play a major role in the future 
expansion of the chart of nuclei. In addition, the sec- 
ondary reactions induced by radioactive projectiles rep- 
resent one of the main methods to produce radioactive 
nuclides with the shortest half-lives, in the nanosecond 
range or shorter. For completeness, we note that in some 
cases a radioactive target may be used for the production 
of exotic nuclei. For example, in the recent discovery of 
a new element with the atomic number Z = 117, the 
radioactive target of ^"^^Bk was used ( Oganessian et al. 



2010) 



Each production method has its own characteristics 
and a typical application range. In the following we men- 
tion briefly various reactions which are being used and we 
direct the reader to papers providing more detailed and 
broader presentations, as well as references to relevant 
technical contributions. Some aspects of reactions used 
to produce exotic nuclei were discussed by [Geissel et al.\ 



( 1995 ) and more recently by Schmidt et al 



(2002). 



In a central collision at low energy two nuclei can join 
together (fusion) to form a single heavier nucleus. In the 
second step, the resulting compound nucleus releases its 
excitation energy by emission of nucleons (evaporation) 
and radiation. The reaction cross section is very sensitive 
to the initial energy in the projectile-target system, which 
must be close to the Coulomb barrier. If the energy is 
too low, the probability of barrier penetration drops dra- 
matically, if it is too large, other channels start to dom- 
inate. This is the key reaction in the synthesis of super- 



heavy elements ( Hofmann 2009b Schmidt and Morawek 



1991). However, since the final nucleus tends to be lo- 



cated on the neutron-deficient side of stability, the fusion- 
evaporation is successfully used to produce very neutron- 
deficient systems. In fact, most of research on proton 



radioactivity employs this reaction (Ferreira and Aru- 



mugan 



2007 Woods and Davis 1997). The most com 



monly used tool for optimizing experimental conditions 
and for prediction of cross-sections is the statistical code 



HIVAP (Reisdorf 1981 Reisdorf and Schadel 1992), but 



other statistical codes like CASCADE (Piihlhofer 1977) 



or PACE (Gavron 1980) are also being used. The poten- 



tial of the fusion-evaporation reaction can be illustrated 
by an attempt to reach a-emitters above ^"^Sn ( [Korgul 



et al. 2008 1 and by an investigation to produce the light- 



est isotopes of bismuth and polonium (Andreyev et al. 



20051 



2. IVlulti-nucleon transfer 

Transfer reactions belong to the category of binary pro- 
cesses where instead of a fused system of two heavy ions 
a projectile-like and a target-like nucleus appear in the 
final state. This happens if the collision is not central. 
If it is also deep-inelastic (damped), a few nucleons can 
be exchanged between reaction partners leading to ra- 
dioactive products. Although a part of the energy of 
the relative motion goes into the excitation of the final 
fragments, which is released by evaporation of light parti- 
cles, still residual nuclei far from stability can be formed. 
Such multi-nucleon transfer reactions at Coulomb bar- 
rier energies has been used to produce unstable nuclides, 
including neutron rich ones ( Broda 2006 1 . The method 



is mainly used in combination with in-beam 7-ray spec- 
troscopy and isomeric spectroscopy in various regions of 



the chart on nuclei ( Cocks et al. 2000 Montanari et al. 



2011 1. The main advances on this field in the last decade 



and the summary of theoretical understanding of the re- 
action mechanism are given in the recent review by |Cor- 
radi et al. (20091. The current limits of nuclear stabil- 



ity cannot be reached by multi-nucleon transfer between 
stable projectile and target, but the importance of this 
reaction is increasing with developments of radioactive 
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beams. The transfer reactions are considered also as a 
tool to produce new isotopes in the region of superheavy 



nuclei (Zagrebaev and Greiner 2008). 



3. Fragmentation 

When the collision energy of two heavy nuclei is large 
compared to the Fermi energy of nucleons, the proba- 
bility that nucleons will be exchanged between the re- 
action partners becomes very small. Instead, violent in- 
teractions occur in the overlapping zone of the projectile 
and the target (participants), while their parts outside 
this zone (spectators) emanate as the projectile-, and 
target-like prefragments, respectively. After this abra- 
sion phase, the cooling of prefragments by evaporation 
of particles, by radiation, or by fission proceeds and the 
final fragments are formed. If the excitation energy of 
the prefragment is large, which happens in more central 
collisions, the multifragmentation takes place, i.e. the 
break-up into many intermediate-mass fragments. In the 
so called limiting fragmentation regime, for projectile en- 
ergies above 100 MeV/nucleon, a characteristic feature is 
observed that the total reaction cross section weakly de- 
pends on projectile energy and can be approximated by 
a simple geometric formula 



nrUA}/' + A]/'-c), 



(6) 



where At and Ap are the mass numbers of the target 
and the projectile, respectively, the radius parameter is 
ro = 1.1 fm, and a correction for nuclear transparency is 
introduced by a parameter c = 2 (Kox et al. 19871. 



In the present context, the fragmentation of the pro- 
jectile plays a special role. When high energy projectile 
ions collide with target nuclei, the projectile-like frag- 
ments surviving the abrasion phase continue moving with 
almost no change of velocity. Thus, the resulting unsta- 
ble nuclei form a secondary beam which can be trans- 
ported and filtered by means of ion-optical devices. This 
method is very fast and universal, since practically any 
nucleus with numbers N and Z smaller than those of the 
projectile can be produced. These features make frag- 
mentation one of the key reactions for radioactive beam 
facilities. The method of projectile fragmentation was 
pioneered at Bevalac facility at the Lawrence Berkeley 



Laboratory (Symons et al. 1979 Westfall et al. 19791 



Later the systematic studies of fragmentation cross sec- 
tions as a function of energy, projectile, and target were 
carried out by Webber et al. (19901. Recently, compre- 



hensive studies of projectile-like fragmentation are be- 
ing carried out at SIS-FRS system at GSI Darmstadt 



(Benlliure et al. 2008 Henzlova et al. 20081. The most 



advanced theoretical description of the fragmentation is 
currently achieved in a modern version of the abrasion- 
ablation model and is implemented as the Monte-Carlo 
code ABRABLA (iGaimard and Schmidtl 119911). The 



mechanism of prefragment excitation is understood in 
this model as a result of random creation of holes in the 
nucleonic Fermi distribution (Schmidt et al. 19931. The 



evaporation stage is modeled with the code ABLA (Ke 



lie et al. 2008). The simpler, analytical version COFRA 
] Benlliure et al. 1999 2000[ ) is applicable to the very- 
neutron rich fragments (cold fragmentation). For practi- 
cal estimates, the empirical parametrization of the frag- 
mentation cross sections is given by a simple analytical 



model EFAX (Siimmerer and Blank 20001. 



4. Spallation 

If in a high energy collision (above 100 MeV/nucleon) 
one of the reaction partners is a light ion, like proton, 
deuteron, or triton, the process is referred to as spalla- 
tion. From the perspective of production of exotic nu- 
clei, the main difference from projectile fragmentation is 
that in case of spallation usually the target nucleus be- 
ing the heavier partner is the source of radioactive nuclei. 
In addition, the mechanism of its primary excitation is 
different. The first step of the spallation is usually de- 
scribed as a series of collisions between nucleons in the 
target nucleus, induced by the projectile, which forms the 
base of intranuclear cascade models (INC). Such a cas- 
cade of collisions leads to a highly excited system which, 
in the second phase, deexcites in the same way as the hot 
target-like prefragment. In consequence, the target nu- 
cleus is destroyed, and in analogy to the fragmentation, 
practically any nucleus with numbers N and Z smaller 
than those of the target can be produced. Due to this 
universality, the spallation is the second main production 
process considered for radioactive beam facilities. In con- 
trast to fragmentation, however, the radioactive products 
have to be extracted from the target material. 

The example of a modern version of an INC approach 
with a discussion of the physics involved is given by 



Boudard et al. (2002). For the second stage of the re- 



action, various versions of statistical evaporation codes 
are being used (Le Gentil et al. 2008). The detailed 



experimental studies of the spallation are conveniently 
performed in the inverse kinematics, where a light tar- 
get (hydrogen, deuterium) is bombarded by heavy nu- 
clei accelerated to relativistic energy. Such studies are 
performed with use of the fragmentation facility which 
reflects a symmetry between these two reactions. This 
method was recently used at GSI Darmstadt in a com- 
prehensive study of residual fragments produced by the 
spallation of ^ssy by protons ( Ricciardi et al. 20061 Taieb 



et al. 2003 ) and by deuterons ( Casarejos et al. 2006 



Pereira et al 
■"^^XTand"^ 



2007). A similar work on spallation of 



Fe by protons was reported by Napolitani 



et al. (20071 and by Le Gentil et al. (20081, respectively. 



Residues of ^^^Pb+x and ^^sy+x at 1 A GeV 
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5. Fission 

Since heavier stable nuclei are more neutron-rich than 
the lighter ones, the process of fission of a heavy nucleus 
is a source of neutron-rich medium mass nuclei. In ad- 
dition, fission is one of the important decay channels of 
excited heavy nuclei. Thus, it plays a role as a direct 
source of exotic nuclei and as a process interfering with 
other reactions used for this purpose. Applications of 
fission to generate neutron-rich nuclei differ in methods 
used to excite a fissile nucleus and in the range of ex- 
citation energies imparted. On the low-excitation end 
is the spontaneous fission and the thermal-neutron in- 



get (Bertulani and Baur 19861. The higher excitation 



duced fission (Rochman ei aZ. 2004 Wahl 19881. To this 



class belongs also fission resulting from electromagnetic 



excitation (photofission) (Cetina et al. 20021. The pho- 
tons inducing fission may be produced directly, e.g. by 
converting an intense electron beam into bremsstrahlung 



(Diamond 19991 or they can be virtual, resulting from 



energies are achieved by bombarding fissile targets with 
beams of fast neutrons or charged particles. Low energy 
proton-induced fission is a frequent choice because of the 
relative technical simplicity. The main aspects of this 



method are discussed by Penttila et al. { 2010 1 who re 



cently developed a novel method to measure the particle- 
induced fission yields. The high energy reactions induced 
by light or heavy ions lead to high excitation energies 
and subsequent fission becomes one of the main deexci- 
tation channels infiuencing the outcome of the spallation 
and fragmentation reactions, respectively. High energy 
reactions in inverse kinematics, where a heavy fissile nu- 
cleus is the projectile, has been proven to be exception- 
ally fruitful. The pioneering experiments with relativistic 
^^^U beam at GSI Darmstadt revealed the properties and 



advantages of this approach (Bernas et al. 1994 19971 



a fast motion of a fissile system relative to a high-Z tar- 



When the target nucleus has a large Z number, the ex- 
citation of a fissile projectile-like fragment has a nuclear 
contribution (fragmentation) and an electromagnetic one 
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FIG. 8 (Color online) A general scheme of the two main meth- 
ods used to extract and separate radioactive nuclei. The 
dashed lines indicate connections which are considered in 
planning future facilities. 



(photofission). This situation was systematically investi- 
gated in the reaction of 1 GeV/nucleon 238-[j impinging 
on a lead target (Enqvist et al. 1999). A similar tech- 



nique was applied in a broad campaign dedicated to the 
comprehensive study of spallation of heavy nuclei. The 
contribution of fission was investigated in detail for the 
spallation of ^ssy by hydrogen ( Bernas et al. 2003 1 and 
by deuterium (Pereira et al. 20071, as well as of ^°*Pb 



by hydrogen ( Fernandez-Dominguez et al. 2005 1 



The emerging general picture of the production of nu- 
clei in both the spallation and the fragmentation reac- 
tions with inclusion of the fission channel was discussed 
by Schmidt et al. (20021. It is illustrated in Figure [t] 
which shows the production cross sections for beams of 
^°^Pb and ^38^ impinging at 1 A GeV on hydrogen, deu- 
terium, and lead targets. Although results of model cal- 
culations are shown in Figure [7] they represent very well 
features observed in experiments. 

Recently, in the new-generation RIBF facility at 
RIKEN Nishina Center, the in-flight fission of a 345 
MeV/nucleon ^ssy bg^m has been used to produce 45 



new neutron-rich isotopes (Ohnishi et al. 20101. 



B. Separation 

In reactions used to produce radioactive nuclei always 
a large number of different products is formed and some 
method of selection is necessary to filter nuclei of inter- 
est from an unwanted background. Very generally, one 
can classify all separation methods which are used into 
two distinct classes. The main difference is the target 
thickness relative to the range of products in the tar- 
get material. If the relative target thickness is small, 
the products emerge from the target with significant ki- 
netic energy and can be promptly manipulated by ion- 
optical devices. The separation techniques of this type 
are called In-Flight and the filtering devices are called 



recoil or fragment separators. On the other hand, if the 
production target is relatively thick, such that the prod- 
ucts are stopped in its volume, the nuclei of interest have 
to be extracted from the target for further filtering. For 
historical reasons such approach is referred to as ISOL 
(which stands for Isotope Separator On Line) technique. 
The principle of operation of In-Flight and ISOL type of 
facility is shown schematically in Figure [8] In the follow- 
ing we discuss briefly only basic features of both these 
methods and provide selected examples of laboratories 
in which they are implemented, referring the reader to 
technical papers with detailed information. A compari- 
son of the two separation methods with a discussion of 
some future prospects, including the hybrid combination. 



was made by Tanihata (20081. 



1. In-Flight 

The key feature of the In-Flight methods is that the 
kinetic energy of the reaction product is large enough to 
escape from the relatively thin production target. This 
method is applicable to reactions induced by heavy ions 
like fusion-evaporation, multi-nucleon transfer and pro- 
jectile fragmentation or fission in inverse kinematics. The 
products emerging from the target enter an ion-optical 
system of magnetic and electric fields where they are sep- 
arated from unwanted contaminants and delivered to the 
final experimental station. Usually, the main selection is 
applied to the mass-over-charge ratio A/q of particles by 
means of a uniform magnetic field, according to the for- 
mula relating the magnetic field B with the momentum 
p of a particle having charge q and moving in this field 
along a circular trajectory of the radius p: 



Bp = - ^ WP7— , 



(7) 



where u is the atomic mass unit, (3 is the particle velocity, 
and 7 is the Lorentz factor (7 — — /3^). Addition- 

ally, in some separators the crossed magnetic and electric 
fields are used to select the velocity of particles (Wien fil- 
ter). The In-Flight method is fast as the typical time of 
fiight through the separator is of the order of a microsec- 
ond or shorter. Another important feature is the lack of 
chemical sensitivity. 

At the low-energy end, the fusion reaction requires pro- 
jectiles with the energy of the order of 10 MeV/nucleon 
(Coulomb barrier). The target thickness is usually of 
the order of 1 mg/ cm^ . The average energy of the reac- 
tion product results simply from the momentum conser- 
vation in a system of collision partners. For this reason, 
the filtering devices in this case are called recoil separa- 
tors. The energy of the recoiling products is usually too 
small to allow for the in-fiight identification of ions. The 
general properties of recoil separators were reviewed by 



Davids (2003) 
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TABLE I The leading laboratories employing the In-Flight method to produce and study the radioactive decays very far from 
the stability. In the lower part of the table the facilities under construction are listed. 



Country Laboratory Driver 



Beams Max beam Separator 



Reference 



accelerator 






energy [AMeV] 


cyclotron 


Ne- 


-Kr 


~ 10 


linac 


H- 


-U 


11 


tandem 


H- 


-U 


~ 10 


tandem + linac 


H- 


-u 


17 


cyclotron 


Li- 


-Ar 


50 


cyclotron 


C- 


-U 


60 


2 cyclotrons 


C- 


-U 


95 


2 cyclotrons 


0- 


-U 


170 


4 cyclotrons 


H- 


-u 


350 


synchrotron 


H- 


-u 


1000 


linac 


H- 


-u 


500 


synchrotron 


H- 


-u 


1500 



Finland 

Germany 

USA 

USA 

Russia 

China 

France 

USA 

Japan 

Germany 

USA 

Germany 



Jyvaskyla 

GSI 

ORNL 

ANL 

FLNR 

HIRFL 

GANIL 

NSCL 

RIBF 

GSI 

FRIB 

FAIR 



RITU 
SHIP 
RMS 
FMA 

ACCULINNA 

RIBLL 

LISE 

A1900 

BigRIPS 

FRS 

Super-FRS 



Leino et al. (19951 



Miinzenberg et al 



Gross et al. 



(20001 



Davids et al. 



Rodin et al. 



Sun et al. 



(19791 



( 1992 1 



(20031 



(20031 



Mueller and Anne 



Morrissey et al. 



Sakurai 



(20081 



Geissel et al. 



Thoennessen 



Winkler et al. 



(19911 



(20031 



(19921 



(20101 



(20081 



By increasing the energy and shifting to the projectile 
fragmentation regime thicker targets can be used. At the 
energy of about 50 ^-MeV the typical targets have a few 
hundred mg/cm^ thickness, while at 1000 ^-MeV they 
reach the thickness of a few g/cm^. At the larger pro- 
jectile energies the kinematical focusing helps to achieve 
larger acceptance by the electro-magnetic fragment sepa- 
rator. Additional filtering, according to the atomic num- 
ber Z of the particle, is achieved by means of an energy 
degrader mounted in the middle of the fragment sepa- 
rator. The high energy of the reaction products allows 
for their full in-flight identification by means of time-of- 
flight and energy- loss measurements for individual ions. 
The resulting extreme sensitivity is one of the most im- 
portant advantages of the medium- and high-energy In- 
Flight technique. Indeed, this was the key factor allowing 
the first observation of doubly-magic nuc lei ^^''SnJLe 
witowicz et a?. [|1994| [Schneider et aZ.l| l994|, "^^N i (|EngeI 



jected into more sophisticated devices like storage rings 



(Nolden et al. 20081 or ion traps (Kluge et al. 20081, 



mann et al. ' 1995 1, and '^'^Ni (Blank et al.\ 20001, as well 



as the discovery of the two-proton radioactivity (Giov- 



inazzo et al. 2002 Pfiitzner et al. 20021. 



Although the transmission, separation, and in-flight 
identification of reaction products are easier at large pro- 
jectile energy, the maximum beam intensity which can be 
achieved decreases with increasing energy. Moreover, the 
quality of the beam of fragmentation products is rather 
poor. The energy spread of the fragments resulting from 
both the reaction kinematics and the energy-loss strag- 
gling in layers of matter hinders efficient stopping in the 
flnal detectors. Hence, the optimal conditions for produc- 
tion of a given nucleus result from a compromise between 
different factors. 

The special role of the high-energy In-Flight method 
comes from a relatively simple way of delivering beams 
of radioactive nuclei. Such radioactive beams may be in- 



or can be used to induce secondary reactions leading to 
a signiflcant expansion of the nuclear physics fleld of re- 
search. 

The leading facilities using the In-Flight method to 
study decays of very exotic nuclides together with their 
brief characteristics and the corresponding reference to 
the detailed information are collected in Table HI 



2. ISOL 

In an ISOL-type facility the nuclei of interest are pro- 
duced in a relatively thick target irradiated by a primary 
beam from a driver accelerator. If the products recoil 
from the target they are stopped by means of a catcher 
or a gas cell, otherwise they diffuse out of the target mate- 
rial. Subsequently they are transferred to the ion source 
and extracted, mostly as 1+ ions, by means of an HV 
potential of the order of 50 kV. Because of the constant 
charge of the extracted ions, the following separation in 
a uniform magnetic fleld corresponds to the mass separa- 
tion, see Eq. [?) In the flrst realizations of this technique, 
such mass separated, very low energy ions were deposited 
on a thin catcher foil in front of a detection system. In the 
modern variant of this technique the mass separated ions 
are postaccelerated and a high quality beam is formed 
allowing better manipulation of the ions and inducing 



secondary reactions ( jJPG 20111. 

Although all kinds of nuclear reactions can (and are) 
employed in this method, the most important are spalla- 
tion induced by protons, and flssion of target nuclei in- 
duced by protons, light ions or neutrons. The latter can 
come either from a reactor or from a beam of deuterons 
hitting a neutron converter in front of the target. An in- 
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TABLE II The main facilities based on the ISOL method for the radioactive decay studies far from the stability. In the lower 
part of the table the facilities under construction or planning are listed. 



Country Laboratory 



Finland 

Belgium 

Italy 

USA 

France 

Canada 



Jyvaskyla 

Louvain-La-Neuve 

INFN-LNS 

ORNL 

GANIL 

TRIUMF 



Switzerland CERN 



Facility 

IGISOL 

LISOL 

EXCYT 

HRIBF 

SPIRAL 

ISAC 

REX-ISOLDE 



accelerator 

cyclotron 

cyclotron 

cyclotron 

cyclotron 

2 cyclotrons 

cyclotron 

synchrotron 



Driver 

beam 
H-Xe, 130q^/A MeV 
H-Ni, 10 AMeV 
A < 48, 80 AMeV 
p, d, Q, 42-85 MeV 
H-U, 95 AMeV 
p, 500 MeV 
p, 1.4 GeV 



Postaccelerator, 
energy [ylMeV] 



tandem, 8 
tandem, ~ 10 
cyclotron, 25 
linac, 11 
linac, 3 



Reference 



Aysto ( 2001 



Kudryavtsev e t al. I 2008 



Cuttone et oZ.|j2008 
2003] 



Stracener 
Villari (2003] 
Shotter 

Voulot et al. (2008 



. (20 03J 
:r (20 ^ 



Italy LNL 
France GANIL 
Canada TRIUMF 
to be decided 



SPES cyclotron p, 70 MeV linac, 11 

SPIRAL-2 linac d, 40 MeV; HI, 14.5 A MeV cyclotron, 25 

ARIEL e-linac e, 50 MeV linac, 11 

EURISOL linac p, 1 GeV linac, 150 



Cinausero et al. I 2009 
Lewitowicz (2008 



http:/ / www.trium f.ca/a riel 
( 2009] 



Blumenfeld et al. 



teresting concept is to induce fission by bremsstrahlung 
photons originating from a very high intensity electron 



beam hitting the high-Z converter (Cheikh Mhamed 
et al.\ 12008). The main point is that the deHvered in- 



tensities of Hght beams, Hke protons, deuterons, or elec- 
trons, can be significantly larger than maximal intensi- 
ties of heavy ion beams. This, in combination with thick 
targets which can be used, results in the high yields of ra- 
dioactive nuclei which is the main advantage of the ISOL 
method. On the other hand, the transfer processes occur- 
ring in the target and in the ion source take time of the 
order of milliseconds (see Figure |2| which imposes limits 
on the half-lives which can be accessed by this method. In 
addition, some of these processes exhibit chemical sensi- 
tivity which for example hinders extraction of refractory 
elements in some implementations of this technique. 

In the last decades, a remarkable progress in ion-source 
techniques and in manipulating low-energy ions has been 



achieved ( Lecesne 2008 Wenander 2008 1 . A spectacular 



example is the application of resonant laser ionization al- 
lowing extremely efficient and clean extraction of selected 



elements from the source (Cheal and Flanagan 20101 



A selection of facilities employing the ISOL method 
with short characteristics and the reference to the corre- 
sponding technical information is presented in Table [IT] 



3. Future facilities 

In the lower part of both Table |l] and Table |ll] the 
future facilities which are being constructed or planned 
are listed. A new idea which is considered in these re- 
cent developments is to combine advantages of both the 
In-Flight and the ISOL methods into hybrid solutions 



(Tanihata 20081. Sufficient postacceleration of the ISOL 



secondary beam may enable taking advantage of instru- 
mentation developed for the In-Flight technique. Such 
an option is discussed in the EURISOL design study 



([Blumenfeld et al. 20091. On the other hand, the fast 
In-Flight fragment beam may be stopped in a gas cell 
and extracted at low energy with help of the ISOL tech- 
niques (Facina et al. 2008 1 . An example of such a solu- 



tion will be realized in the currently constructed FRIB 
facility (Thoennessen 2010). The possible connections 



between the two main approaches to extract and sepa- 
rate radioactive nuclei are marked in Figure |8] by dashed 
lines. 



C. Detection 

The detection of a radioactive decay requires detection 
of particle(s) emitted in the process. The large major- 
ity of nuclides decay by /3 transitions where primarily an 
electron or a positron is emitted (the presence of neu- 
trinos can be safely neglected in this context) and in the 
second step electromagnetic radiation follows in form of 7 
radiation if the daughter nucleus was formed in an excited 
state and/or of characteristic X-rays if the final atom 
was excited. This secondary electromagnetic deexcita- 
tions may proceed by emission of Auger electrons. When 
we move away from the stability line, however, the mass 
differences between isobars (and thus the beta decay en- 
ergy windows) increase and the particle unbound states 
become populated. On the neutron-deficient side this 
leads to the /3-delayed charged particle emission, mainly 
of protons. Beyond the proton-drip line the direct emis- 
sion of protons comes into play. That is why the detec- 
tion of particles like p and a plays the central role in 
the radioactivity studies at the neutron-deficient limit of 
stability. In turn, on the neutron-rich side emission of 
/3-delayed neutrons becomes important which makes the 
neutron spectroscopy necessary. However, also charged 
particles, like d, t, and a, are emitted following /3 decay of 
very neutron-rich nuclei. In the following we sketch the 
modern methods used to detect charged particles and 
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Some abbreviations used in Table 



|l]and|ll] 



ANL Argonne National Laboratory, USA 

FAIR Facility for Antiproton and Ion Research, Darm- 

stadt, Germany 

FLNR Flerov Laboratory of Nuclear Reactions at Joint 
Insitute for Nuclear Research, Dubna, Russia 

FRIB Facility for Rare Isotope Beams at Michigan State 

University, East Lansing, USA 

GANIL Grand Accelerateur National d'lons Lourds, Caen, 
France 

GSI Helmholtzzentrum fiir Schwerionenforschung 

GmbH, Darmstadt, Germany 
HRIBF Holifield Radioactive Ion Beam Facility at ORNL, 

USA 

HRIFL Heavy Ion Research Facility in Lanzhou, China 
INFN-LNS Istituto Nazionale di Fisica Nucleare, Laboratori 

Nazionali del Sud, Catania, Italy 
LNL Laboratori Nazionali di Legnaro, Legnaro, Italy 

NSCL National Superconducting Cyclotron Laboratory 

by Michigan State University, East Lansing, USA 

ORNL Oak Ridge National Laboratory, Oak Ridge, USA 

RIBF Radioactive Isotope Beam Factory at RIKEN lab- 

oratory, Wako, Saitama, Japan 

BigRIPS Big RIKEN Projectile Fragment Separator at 
RIBF 

EXCYT EXotics with CYclotron and Tandem at INFN- 
LNS 

FMA Fragment Mass Analyser at ANL 

FRS FRagment Separator at GSI Darmstadt 

IGISOL Ion Guide and Isotope Separator On-Line 
LISE Eigne d'lons Super Epluches at GANIL 

LISOL Leuven Isotope Separator On-Line 
RIBLL Radioactive Ion Beam Line in Lanzhou at HRIFL 
RITU Recoil Ion Transport Unit 
RMS Recoil Mass Spectrometer at HRIBF 

SHIP Separator for Heavy Ion reaction Products at GSI 

Darmstadt 

SPES Selective Production of Exotic Species at LNL 

SPIRAL Systeeme de Production d'lons Radioactifs avec 
Acceleration en Eigne at GANIL 



neutrons. Finally, we discuss briefly the digital signal 
processing techniques which represent a new development 
in the systems of nuclear data acquisition. 

The technique of 7 spectroscopy is the large subject in 
its own, as it is the main source of detailed data on the 
nuclear excited states. Here we refrain from discussing 
its development referring the reader to the rich literature 



1. Charged particles 

The most common devices for detecting charge parti- 
cles are based on silicon detectors which record the en- 
ergy deposit of a charge particle passing through its ma- 
terial. To increase the sensitivity of the measurement 
and to provide an additional information about the par- 
ticle, stacks of two or more detectors are frequently used. 
In such a telescope the energy loss information from thin 
transmission detectors is combined with the total kinetic 
energy from the final thick detector. A particle telescope 
may combine a thin gas chamber, acting as a transmis- 



sion counter with a thick silicon detector ( Axelsson et al. 



1998 Moltz et al. 19941. Additional information on the 



location of the particle's hit can be extracted from a po- 
sition sensitive detectors where the signal is read from 
two ends of a resistive electrode. A significant advance 
in detection technique was introduced with a concept 
of a silicon strip detector. Particularly successful was 
the development of a double sided silicon strip detector 
(DSSSD) with perpendicular sets of strip electrodes on 
its both sides ( Sellin et al. 1992 ) . The achieved granular- 



ity provides not only a simple measure of the position but 
allows to establish a position correlations between subse- 
quent events, like the implantation of an ion and its decay 
by particle emission in condition of the high total count- 
ing rate. In addition, it reduces the effect of energy sum- 
ming between /3-delayed particles and electrons which de- 



posit much less energy in a strip (a pixel) area (Biischer 



et all 2008). Most of results on /3-delayed proton emis- 



sion (Sec. IV I and on proton radioactivity (Sec.|V|) were 



obtained with help of such detectors. An example of a re- 
cent improvement is a novel design of a large-area DSSSD 
with an ultra-thin dead layer ( Tengblad et al. 2004 1 . The 



modern detection set-ups which require granularity but 
also a large angular coverage are usually constructed as 
arrays of silicon detectors. They may consist of a la rge 
number of simple Si diodes ( Fraile and Aysto 2003 1 , a 



box of DSSSD detector (Adimi et al. 20101, or a combi 



nation of DSSSD detectors with a gaseous multiwire pro- 



portional chamber and germanium detectors (Page et al. 



2003[ ). In the decay studies at projectile fragment sepa- 



rators (the In-Flight method) the large area of the final 
focal plane and the large range distributions of selected 
ions have to be taken into account. To meet the lat- 
ter challenge, the stacks of many DSSSD detectors are 
used. The example solution used for the /3-decay studies 
of ^°"Sn at the FRS separator (GSI Darmstadt) consisted 
of three large-area DSSSD detectors, providing in total 
3 X 60 X 40 = 7200 pixels, sandwiched between two sets 



of ten single-side silicon strip detector (Eppinger et al. 



20091. In another development for the FRS separator, a 



set of two rows consisting of three 5 cmx5 cm DSSSD 
detectors, each with 256 (16x16 strips) pixels was de- 



2006 Lee et al. 20031. 



on this topic (Farnea et al. 2010 Gelletly and Eberth veloped (Kumar et al. 20091. Even more ambitious Ad- 



vanced Implantation Detector Array (AIDA) , to be used 
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at the FAIR facility, will comprise twenty four 8x8 cm^ 



counting gas at atmospheric pressure 



DSSSD detectors with 128 x 128 strips (Davinson 20101 



The implantation of a nucleus which undergoes a mul- 
tiparticle decay into a silicon detector has a serious lim- 
itation that only the total decay energy can be mea- 
sured. The energy sharing between products and their 
momentum correlations cannot be accessed. Such prob- 
lem appeared in case of the two-proton radioactivity. The 
first observation of this decay mode was accomplished by 
implanting ions of ^^Fe into a stack of silicon detectors 
fPfiitzner aZ.[|2002| ) and into a DSSSD detector ( |Giov 



inazzo et al. 20021. The information on the decay time 



and on the total decay energy were sufficient to claim the 
new type of radioactivity, but for the detailed study of 
this process a novel experimental approach was necessary. 
To meet this challenge, two new developments were un- 
dertaken, both based on a principle of the time projection 
chamber (TPC) . The key idea is that such a gaseous ion- 
ization chamber can record tracks of charged particles, 
allowing their reconstruction in three dimensions. The 
radioactive ion stopped inside the active volume and the 
subsequently emitted particles ionize the counting gas. 
The primary ionization electrons drift in a uniform elec- 
tric field towards the charge amplification section produc- 
ing the two-dimensional representation of the particles' 
tracks. The drift time contains the position information 
along the electric field direction. In one solution, the 
amplified ionization charges are collected electronically 
by means of an anode plate with two sets of orthogonal 
strips (Blank et al. 20101. This detector rendered the 



first direct evidence for the two protons emitted in the 
decay of ^^Fe (Giovinazzo et al. 2007). In the second 



design, the idea of an optical readout (Charpak et al 



1988) was implemented. It is based on the observation 



that light is emitted in the final stage of charge amplifica- 
tion. In the Optical Time Projection Chamber (OTPC) 
this light is collected by a CCD camera and by a photo- 



multiplier (PMT) ( [Miernik e< am2007a| ). The construc- 
tion of this detector is shown schematically in Figure |9] 
The application of the OTPC yielded spectacular results, 
including the detailed proton-proton correlation picture 
for the 2p decay of '"^Fe (iMiernik et al.\ |2007c|) and the 



first observation of the /3-delayed three-proton emission 



channel (Miernik et al. 2007b). An example event of the 
two-proton radioactivity of *^Fe is shown in Figure 10 



In case of very short decay half-lives, in the sub- 
nanosecond range, the implantation technique generally 
cannot be used. A short-lived precursor decays in-flight 
very close to the place of its production. The identifica- 
tion of the nucleus and its properties can be deduced from 
the detection and tracking of all decay products. This 
approach was successfully appUed to the study the 2p 
decay of ^^Mg (Mukha et al. 2007) which also exempli- 
fies advantages of radioactive beams. Ions of ^^Mg were 
produced in a secondary target by a neutron knock-out 
reaction from a beam of ^°Mg delivered by the GSI FRS 
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charge 
amplification 
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FIG. 9 (Color online) A schematic view of the Optical Time 
Projection Chamber (OTPC). For each recorded event, the 
data consist of a 2D image taken by a CCD camera in a 
given exposure time and the total light intensity detected by 
a photomultiplier (PMT) as a function of time, sampled by a 
digital oscilloscope. The gating electrode is used to block the 
charge induced by incoming ions. 



separator ( Geissel et al. . 1992 1 . The tracking of emitted 



protons by means of silicon microstrip detectors ( Stanoiu 



et al. 2008) allowed to establish the longitudinal distri- 



bution of decay vertexes and to determine the half- life of 
^^Mg to be 4.0(15) ps. At the same time the information 
on correlations between emitted protons was collected. 
Since the beam impinging on the secondary target con- 
tains usually a mixture of different ions ("cocktail" beam) 
other reactions can be addressed simultaneously. For ex- 
ample, in the measurement of ^^Mg, the data on proton 
and two-proton decays from ^^F, ^^Ne, and ^^Na were ob- 
tained (Mukha et al. 2010). Similar technique has been 



applied t o study two-pro t on emission form excited state s 
of i^Ne (IChroinik et al. 



2002 



Zerguerras et al. 20041. 



The tracking method of the in-flight decay products is 
expected to provide information on several 2p emitters 



among light nuclei, see Sec. VII 



2. Neutrons 

Currently, two different methods for neutron detection 
are used in nuclear spectroscopy. The flrst one is based 
on thermal-neutron induced reactions, like ^He(n,p)'^H, 
^Li(n,Q;)'^H, or ^°B(n,a) ''Li, leading to charged particles 
which can be easily detected. The neutrons emitted by 
a radioactive source have to be flrst thermalized and 
this is achieved by means of a moderator — usually a 
large block of polyethylene surrounding the source. In 
the moderator cylindrical cavities, arranged in concen- 
tric rings, are drilled in which proportional counters are 



mounted (Mehren et al. 1996). In these counters which 
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FIG. 10 (Color online) An example of a registered two-proton 
decay event of *^Fe. Top: an image recorded by the CCD cam- 
era in a 25 ms exposure. A track of a *^Fe ion entering the 
chamber from left is seen. The two bright, short tracks are 
protons of approximately 0.6 MeV, emitted 535 after the 
implantation. Bottom: a part of the time profile of the total 
light intensity measured by the PMT (histogram) showing in 
detail the 2p emission. Lines show results of the reconstruc- 
tion procedure yielding the emission angles i9 with respect to 
the axis normal to the image. From Miernik et al. (2007c I. 



are filled with ^He or BF3 gas, neutron-capture reactions 
take place and are detected. Such construction allows to 
cover a large solid angle, approaching 47r and the large 
total effi-ciency of up to 30% can be achieved for a broad 
neutron energy-range from nieV to tens of MeV and al- 
most independent on the neutron energy due to thermal- 
ization. Since the information on energy is lost, such 
detector is used primarily for counting which makes it 
well suited for determination of branching ratios for vari- 
ous neutron emission channels. Another disadvantage of 
the moderation is that a neutron is detected up to about 
100 /is after the emission. Such a delay reduces the to- 
tal counting rate which can be accepted. An example 
of a modern version of such a 4 tt neutron counter is the 
NERO detector, recently built at the NSCL laboratory 
( Pereira et al. 2010 1 . Its layout is shown in Figure 
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FIG. 11 (Color online) Schematic drawings of the NERO de- 
tector. Left: side view showing the Beta Counting Station 
(BCS) chamber located inside of NERO with the DSSSD at 
the central position. Right: backside showing the cylindri- 
cal cavity to house the BCS and the three concentric rings of 
gas-filled proportional counters. The labels A, B, C and D 



designate the four quadrants. From Pereira et al. (20101 



actions of fast neutrons are detected, predominantly by 
recording elastic proton recoils. The neutron energy can 
be then determined by means of the time-of-flight (TOF) 
after a trigger signal, given for example by a /3 par- 
ticle. The panels containing liquid or plastic scintilla- 
tors are mounted at some distance from the radioactive 
source, which usually reduces the solid angle which can be 
achieved. The efficiency depends on the neutron energy 
and exhibits a low-energy threshold at about few hundred 
keV. In addition, such detector is sensitive also to 7 radi- 
ation and the pulse-shape analysis has to be performed 
for the n-7 discrimination ( Skeppstedt et al. 1999 1 . Such 



an approach to the neutron TOF spectroscopy can be ex- 
emplified by the detector TONNERRE developed at the 



GANIL laboratory (Buta e< aZ. 20001 



A different, and to a large degree a complementary 
solution, employs scintillation detectors in which inter- 



3. Signal processing 

In the conventional approach, signals from detectors 
of nuclear radiation are preamplified and then processed 
in analogue-electronics modules like shapers, amplifiers, 
discriminators, etc., to be finally converted to the digi- 
tal form in analogue-to-digital (ADC) and time-to-digital 
(TDC) converters, and stored in the electronic memory of 
the data acquisition system. With the increasing number 
of channels which have to be read, resulting from pixeliza- 
tion of the detectors (strips, segments, pads) the amount 
of necessary electronic units is growing and thus magni- 
fying the complexity and the cost of the instrumentation. 
In addition, by storing only the values of energy and time 
for an event, the information on the pulse shape is lost, 
which is very disadvantageous in some applications. 

A possible solution to these problems is offered by the 
technique of the digital signal processing (DSP) which 
since recently is taking over the conventional data acqui- 
sition systems in nuclear spectroscopy. Its basic princi- 
ple is that the output of a preamplifier is digitized first 
and all further manipulations are performed by numer- 
ical algorithms acting on this digital representation of 
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FIG. 12 Part of the preamplifier signal waveforms recorded 
by the front and the back strip of the 65-^m DSSSD during 
the ^*^Tm experiment. The recoil depositing about 14 MeV 
energy is followed after 0.55 pis by the 1.73 MeV signal. From 
Karny et al\^003^. 



the signal. These algorithms replace all functions of ana- 
logue electronics and additionally offer a choice of much 
more complex and flexible operations on the pulse. Orig- 
inally, the introduction of the DSP methods was mo- 
tivated by needs of segmented X-ray and 7-ray arrays 
and this sector is the main recipient of this technology 



(Crespi et al. 


2009 ICromaz et al. 


2008 IPietri et al.\ 


2007 


Starosta et al. 


2009 


1. One commercial develop- 



ment — the Digital Gamma Finder module (DGF-4C) 
by XIA LLC ([Hubbard-Nelson et al.\ |1999[)— proved to 



be particularly successful also in the domain of particle 
spectroscopy. Some applications of the DGF electronics 
in various decay studies of exotic nuclei were reviewed 



by Grzywacz ( 2003 1 . The trends and new products from 
this developer were presented by [Warburton and Grud 
beFglp006|. 



A good illustration of new possibilities provided by the 
DSP is the measurement of very short-lived proton ra- 
dioactivity in ^^^Tm ( Karny et al. 2003 ) . The technical 
challenge is to detect a low-energy proton 1.5 MeV) 
emitted very shortly (~ 1 /zs) after stopping of the parent 
nucleus which deposits up to about 35 MeV in the implan- 
tation DSSSD detector. Such sequence of events cannot 
be resolved easily when signals are passed through the 
analogue amplifiers. The solution offered by the DGF 
electronics is to store the whole waveforms of the sig- 
nals from the silicon strips and to analyze their shapes 
off-line. In addition, the special triggering mode was im- 
plemented to the DGF board which allowed to store only 
those events in which the pile-up of two pulses was de- 
tected. This feature leads to the large increase of sen- 
sitivity, which is especially important when many differ- 
ent ions are coming to the final detector and the decay 
investigated is rare. This technique was the key factor 



leading to the discovery of the fine structure in the decay 
of ^"^^Tm (Sec.|v|. An example of the recorded waveform 
representing a low-energy proton signal superimposed on 
the large implantation signal of the ^^^Tm ion is shown in 
Figure [12] The same method was instrumental in the ob- 
servation of the superallowed a-decay chain from ^°^Xe 



(Darby et al. 2010) which is discussed in Sec. VI 

Another development which is recently gaining impor- 
tance, especially in the domain of detectors with high 
granularity, is the technology of so called application- 
specific integrated circuit (ASIC) . It is based on a highly 
integrated circuit which is customized for a specific use 
rather than for a general-purpose application. Usually, 
one integrated circuit (IC) chip features several indepen- 
dent channels, each capable of handling energy and tim- 
ing of a single detector element (pixel or strip). In fact, 
the large silicon array AIDA ( Davinson] 2010) as well 
as one of the TPC detectors developed to study 2p ra- 



dioactivity (Blank et al. 2010), see Sec. III.C.l employ 
ASIC-type chips in their read-out electronic system. An 
example of multi-channel IC for the detection of nuclear 



radiation is described by Engel et al. (2011). In combi 



nation with an array of Si-strip detect ors ([Wallace et al. 
2007) it was used in a study of ^Be (Grigorenko et al. 



2009a I, presented in Sec. VII.B.l 
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IV. BETA DELAYED PARTICLE EMISSION 
A. Beta decay, general observations 

1. The beta strength 

The weak interactions and in particular their low- 
energy manifestation in nuclear beta decay are by now 



well understood ( 


Behrens and Biiring 


1982 


Grotz and 


Klapdor 


1990 


), see ( 


Severijns et al. 


2006 


1 for a recent 



survey of weak interaction tests in nuclear physics. The 
decay rate for allowed /?~ or /3+ decays can be trans- 
formed to give the known expression for the /t-value 



ft 



K 



9aBgt ' 



(8) 



where t is the partial halflife of the transition, K/g^ 



6144.2 ( 1.6) s and gA/gv = -1.2694(28) ( [Towner and 



Hardy, 20101, and Bp and Bgt are the reduced matrix 
elements squared for the Fermi and Gamow- Teller parts. 
(Note that some authors define Bqt to include the factor 



the strength distribution will be given in the following 
two subsections. 

For the lightest nuclei the Q-values are for a given 
mass number slightly higher on the proton-rich side than 
on the neutron-rich side of beta-stability and experimen- 
tally halflives are systematically shorter on the proton- 
rich side. The asymmetry is enhanced by the contribu- 
tion from the IAS transition in nuclei with N < Z. For 
masses above 100 the situation has changed: experimen- 
tally the halflives for nuclei more than 3-4 nucleons away 
from the beta-stability line are systematically shorter on 
the neutron-rich side than on the proton-rich side. For 
these nuclei the /3+ decay increases isospin and the sys- 
tematic difference can be understood from eq. ([9| since 
the summed f3~ strength is significantly higher than the 
/3+ strength. 

Several approaches have been employed to reproduce 
and predict the beta decay halflives in large parts of 



the nuclear chart (Borzov and Goriely 


2000 I 


iomma 


et al. 


1996 


Moller et al. 


2003 Nakata et al. 


19971. 



The increase in computing power has allowed the use 



(qa/qv) ■) Nuclear electron capture will also contribute, r- ■ , u- j.- j. j ™- • ™ j i /To 1 

\nAiyvj / f , 01 mcreasmgly sophisticated microscopic models ( Borzov 



but is mainly noticeable for low decay energies and in 
heavier elements. The phase space factor can be approx- 
imated roughly by / = (1 + Q / mf,(?)^ / ?iQ in terms of the 
decay energy (Q-value) for /3^-decays, more accurate de- 



et al. 2008 Borzov 2006 



methods (Navratil et al. 



Brown 20011, even ab initio 



2007 Pervin et al. 



the very light nuclei, and more extended RPA ( Toivanen 



terminations can be found through ( Behrens and Biiring 



20071 for 



et al. 2010 1 and shell-model calculations are underway 



1982 Wilkinson 1995|. Due to the parabolic mass sur- 



face the Q- values increase as one moves away from beta- 
stability which by itself enhances beta decay rates. Em- 
pirically the beta halflives fall off approximately expo- 



20071. 



nentially away from stability ( Zhang et al. 

There are only few beta decays where most of the beta 
strength is energetically accessible in the decay, and the 
detailed distribution in general plays a crucial role. The 
Fermi strength is concentrated around the Isobaric Ana- 
logue State (IAS). The summed strength fulfils a sum 
rule ^ B^ — ^ B^ — Z ~ N that involves Fermi tran- 
sitions in "both directions", this is relevant e.g. for the 
odd-odd N = Z nuclei from ^^Cl to ^®In where most 
ground states have T = 1 and can be fed as well as decay 
by Fermi transitions. 

The Gamow- Teller strength obeys the Ikeda sum rule: 



E^GT-E^GT = 3(iV-Z), 



(G. Martmez-Pinedo, private communication, 2010). 



2. Delayed particles 

Defining as usual the relative probability P5 for a given 
beta-delayed decay mode S (Ip, 2p, In, 2n, a etc), as 
the fraction of all decays that results in a final state con- 
taining S, one can find the average number of emitted 
neutrons as P„ = iPin with Pp defined in a similar 
way. The energetics for the different channels is sketched 
in figure [13] The figure implicitly assumes that decays 
take place through states in the emitter and that multi- 
particle emission happen sequentially; as will be argued 
below these assumptions may not hold for all cases. The 
Q-value for delayed emission of neutrons explicitly de- 
pends on the neutron separation energies, but also the 



(9) Q- value for some other modes can be rewritten (Jonson 



and much of this strength is collected in the several MeV 
wide so-called Gamow- Teller Giant Resonance (GTGR), 
although mixing with higher lying configurations is im- 
portant and removes a sizeable part of the strength to 



and Riisager 2001 1 to show that they depend on the 



higher energies (Ichimura et al. 20061; this is often re 



ferred to as the quenching of the GT strength, the key 
point being that the quenching factor, although depend- 
ing on what approximations are used when calculating 
the GT strength, is varying slowly as a function of N 
and Z. Note that the summed value of Bgt, even in- 
cluding quenching, is larger than Bp. More details on 



neutron separation energies of the precursor nucleus: 
Q0d = 3.OO7MeV-52„, g^t = 9.264MeV - 53„. (10) 

Most beta-delayed decay modes will therefore be en- 
hanced at the driplines since multi-nucleon separation 
energies will be low there: the "dripline" for emission of 
two or more neutrons will lie very close to the one neutron 
dripline. We shall return in section [TV.D[ in more detail to 
the particle emission processes, but can note already now 
that the probability for a delayed multi-particle emission 
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FIG. 13 Some of the possible energetically allowed beta- 
decay channels for a neutron rich nucleus. The precursor, 
^Z, beta decays to the emitter, ^{Z -\- 1), that has particle 
unbound excited states. All energies are given relative to the 
ground state of the emitter. 



may depend on the emission mechanism (simultaneous or 
sequential emission) as well as on the energy available. 

As we approach the driplines, the enhanced role played 
by beta-delayed particle emission implies that the physics 
problems investigated via beta decay will overlap more 
and more with the ones investigated via reaction studies, 
but the selection rules for beta decay will provide a spin 
selectivity that often is useful. We shall focus here on the 
general features of the beta decay processes and specific 
challenges met in decay studies, but shall give as well 
selected examples of structure questions that have been 
studied. 

A specific example of this is the population in beta 
decay of excited states that enter in the astrophysical 
rp-process (Wormer tt al. 19941. An even more direct 



need for beta-decay data in astrophysics arises in pro- 
cesses where weak interactions play a role, either directly 
as beta decay rates or indirectly where neutrino inter- 



actions are important, see ( 


Arnould et al. 


2007 


Borzov 


2006 


JGrawe et aL[|2007[|Langanke and Martuiez-Pinedo 


2003 





Quite apart from the general interest in the coupling to 
continuum degrees of freedom, cf. section [l.C.3| beta de- 
cay processes may provide specific information on isospin 
mixing that is expected to be enhanced for continuum 



states at low energy, see e.g. ( 


Garrido et al. 


2007 


Michel 


et al. 2010b 


Mitchell et al. 


2010 





The experimental considerations were covered in gen- 
eral in section [Till but a few specific comments may be rel- 
evant. Since beta halfiives are longer than about one ms, 
essentially all experiments make use of a stopped beam. 
This gives a source distribution that at in-flight facilities 
may have a considerable spatial extent. Experiments in 



storage rings (Litvinov and Bosch 2011 ) or ion and atom 



eral cases, but still present practical problems in partic- 
ular concerning efficient detection of all decay products. 
Complementary experiments at different type of facilities 
may overcome the disadvantages for a specific produc- 
tion method, one example being the study of ^^Ar (Bhat- 



tacharya et al. 2008 1 where a high-resolution spectrum 
obtained at ISOLDE/CERN was combined with an abso- 
lute intensity determination carried out at NSCL/MSU. 
ISOL facilities often have problems for very short halfiives 
and in determining absolute activities, whereas in-fiight 
facilities frequently employ composite beams (so-called 
cocktail beams) where special procedures may be needed 
in order to correct for background from decay of non- 



relevant isotopes, see e.g. (Dossat et al. 2007 Kurtukian- 



Nieto et al.\ 2008). As will be seen below many results 
from the past decade come from in-fiight facilities, of- 
ten through implantation of the radioactive ion into a Si 
detector. 

Beta-delayed particle emission has been the subject of 



several earlier reviews, both more general ones (Jonson 



and Nyman 1996 Jonson and Riisager 2001) and spe 



cific ones for proton-rich nuclei (Blank and Borge 2008 



Hardy and Hagberg 1988), neutron-rich nuclei (Hansen 



and Jonson 1988) and heavy nuclei (Hall and Hoffman 



1992). Since more detailed accounts can be found there 



the treatment here, in sections [l V . B [ [TVTC] and 



Wn will 



be somewhat brief. The remaining section |IV.D deals 



traps ( Severijns et al. 2006 1 have been undertaken in sev- 



with beta-delayed multi-particle emission and naturally 
has more interconnections to other parts of the present 
paper. 



B. 13'^ delayed emission of one particle 

1. Occurrence of particle emission 

Figure [14] shows Q ec and the beta-decay halfiives for 
the most proton-rich nuclei where beta-decay still plays 
a role. There is considerable scatter in the values, but 
also clear effects of the proton shells at Z = 50, enhanced 
by the fact that Fermi transitions contribute below this 
value and not above, and Z = 82 as well as the neu- 
tron shell at 82 corresponding to Z = 72: below this the 
competing decay mode is proton emission, above alpha 
decay takes over inside the proton dripline. For beta de- 
cays along the proton dripline (dashed lines in the figure) 
one still finds that the Q- values decrease towards 10 MeV 
in the heavier nuclei. The scatter indicates that local nu- 
clear structure still plays an important role in these de- 
cays. Even though protons and neutrons in many cases 
still are within the same major shell, forbidden decays 
will play an important role for the heavier nuclei. 

The Coulomb barrier plays a dominating role in beta- 
delayed particle emission in proton-rich nuclei, as illus- 
trated in figure [TS] and essentially limits the emitted par- 
ticles to protons. Delayed alpha emission is energetically 
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FIG. 14 Q_BC-value (upper figure) and beta-decay halflife 
(lower figure) as a function of proton number Z for the light- 
est nucleus for each element where beta-decay remains the 
dominating decay mode. Dashed lines (for 73 < Z < 82) are 
the values at the proton dripline. Experimental values from 



( Achouri et 0^120061 [Audi et al., 2003a"b': T) ossat et al. \ |2007 1 
are completed by estimates based on ( jMoller et aZ.[|1997| ^ 



allowed for many proton-rich nuclei and may seem en- 
ergetically favoured for nuclei above Z = 50 where the 
beta-daughter often has a positive Qa value (see figure 
151. However, this decay mode is mainly important in 
light nuclei: the only nucleus above mass 20 where the 
/?« branching ratio gets above 1% is -'^^"I. The /?« pro- 
cess has in heavier nuclei mostly been observed just above 
closed shells similar to what is seen for ground state al- 
pha decays. The competing process to delayed proton 
emission is therefore delayed gamma emission. The retar- 
dation from the Coulomb barrier will also be significant 
for protons, but the staggering of the proton dripline of 
course implies that there will be nuclei with sizable 
for most even Z-values. Turning to nuclei that lie within 
the "odd-Z dripline" it appears, with our present incom- 
plete experimental knowledge, that beta-delayed proton 
emission with Pp above 1% with few exceptions occur in 
nuclei that are at most one or two nucleons away from 
the line. One may thus regard significant beta-delayed 
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FIG. 15 The centre-of-mass energy that gives an s-wave pen- 
etrability of 10"^ (fuh hues), 10"^ (dashed hues) or 10"^° 
(dotted fines) for a beta-delayed proton or alpha particle are 
shown versus the charge Z for the precursors shown in figure 
|14| Gamma emission can be expected to compete for penetra- 
bilities be low 10~^ (c f. fi gure 21. For illustration the emitter 
Qa value ( [Audi et al.\\20Q^^ is shown for a few beta-decaying 
nuclei, see the text. 



proton emission as a dripline phenomenon. To give one 
example, the nucleus ^^^Ir has a ground state and an 
isomer that both decay by proton emission, alpha par- 



ticle emission and beta decay (Davids et al. 1997), but 
beta-delayed particle emission has not been reported even 
though the proton separation energy is below 2 MeV in 
the daughter nucleus ^^''Os (however, such events may 
be harder to see with the tagging technique employed in 
the experiment). 

The competition between proton and gamma emis- 
sion can lead to the occurrence of gamma-delayed proton 
emission. The angular momentum barrier for outgoing 
protons seems to make this happen frequently in high- 



spin physics (Rudolph 2002), but it may also happen 



after beta decays where angular momentum barriers are 
smaller. It has been suggested to take place in the decays 
of 2°Na dClifford et al\ |l989l iKirseboml |2010|) and ^^Ar 



( Wrede et al. 2009 1 , but may be expected also in other 
nuclei. Gamma emission preceding particle emission is 
well-known in the light nuclei and recent dedicated reac- 
tion experiments have now succeeded to observe it even 
for cases where one or both of the unbound sta tes are 
broad, namely for the 4+ to 2+ transition in ®Be (Datar 



et al. 20091. 



et al. 2005) and several transitions in ^^C (Kirsebom 
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FIG. 16 The total branching ratio of beta-decays to the IAS 
is shown as a function of relative proton excess for the light 
Ne and Ar isotopes. The dashed line gives the estimated 
branching ratio for Fermi beta-decays to the IAS for the Ni 
isotopes. Only the Fermi part of the transition is included, 
the partial halflife is assumed to scale inversely with Z — N 
and total halflives are taken from ( [Dossat et aL[[2"007[ ). 



2. Fermi decays 

For Z up to 50 the dripline nuclei have N < Z so the 
Fermi strength contribute to beta decay. The approxi- 
mate model-independence ol Bp makes the IAS transi- 
tion interesting even though it, as shown in figure |16| 
only dominates the decays close to N = Z where the IAS 
is at low excitation energy. The decay rate for the transi- 
tion to the IAS is proportional to {Z — N) f p{I^Ec) where 
the Coulomb energy shift /S.Ec depends only slowly on 



mass number for a given set of isotopes (Antony et at. 



1997). Even though this strength increases with Z — N 
the branching ratio to the IAS will decrease. Further- 
more since the IAS will be situated at higher and higher 
excitation energy its decay will become more fragmented 
and there will for the most proton-rich nuclei such as 
^''Ne and ■^^Ar not any longer be a dominant IAS peak 
in the final state spectrum. 

A first approximation of the wavefunction of the IAS 
will be \A >, the (normalized) state obtained by letting 
the isospin lowering operator work on the parent state. 
However, the two will not be exactly identical, a fact 
often referred to as isospin symmetry breaking. There 
are two aspects of this: the isospin of the IAS will not 
be pure (isospin mixing) and the radial wavefunctions 
may differ slightly (imperfect overlap, here isospin is in 
principle conserved) . The magnitude of the combined ef- 
fect has been an important issue in precision determina- 



matrix element Vud (Towner and Hardy 20101. Models 
predict incre ases from around 0.2 % for ^'Cto close to 
2% for ^^Rb dGrinyer et aZ.||2010|) and it has been mea- 



sured to 2.0(4)% in ^^Ar ( [Bhattacharya et al||2008| and 
can be expected to be of similar magnitude along the 
dripline. Isospin mixing is known from reaction studies 
to be small in the sense that the spreading width of an 



IAS typically is in the range 10 to 100 keV, see (Harney 



et al. 1986 Mitchell et al. 20101 and references therein 



Note that many of the present beta-delayed proton exper- 
iments may not be able to resolve the IAS from close-by 
levels that it mixes with; the average level spacing can 
be estimated from the mirror systems where proton scat- 
tering on even-even nuclei give about 10 keV for nuclei 



at mass 50 (Bilpuch et al. . 19761. 



To discuss the Fermi strength distribution in more de- 
tail one can start from the simple relation 



<A\H\A>^^E,\ <i\A> 



(11) 



where H is the Hamiltonian of the system and \i > a, com- 
plete set of states. Although |A > is not an eigenstate of 
the system, the expectation value on the left is the quan- 
tity that formally enters in the isobaric multiplet mass 
quation (IMME). However, the IMME is very resilient 



(Benenson and Kashy 1979 Bentley and Lenzi 20071 



and will hold in many cases even though isospin sym- 
metry breaking may be significant in the intermediate 
multiplet members. The Fermi strength that is spread 
out via isospin mixing will obviously remain close to the 
IAS, the important consequence of eq. (Ill being that 



this also holds on average for the Fermi strength spread 
out due to imperfect overlap between \A > and the IAS. 
The redistribution of Fermi strength has been checked 
experimentally in a few decays, e.g. -^^Na (Clifford et al. 



1989 1 , and it would be interesting to have thorough stud- 



ies in more nuclei where the effects are expected to be 
large such as when continuum effects become important. 

The particle emission from the IAS will in most cases 
be isospin forbidden (Auerbach 1983 Brown 1990} and 



the width will consequently be so narrow that gamma 
decays may have a substantial branching ratio. This is 
well established in light nuclei and must be kept in mind 
when detailed investigations of the Fermi strength be- 
come possible also in heavier nuclei. 



3. Ga mow- Teller decays 

The GTGR will for proton-rich nuclei lie above the 
IAS, but can be reached in beta decay e.g. for the hght- 



est Ar (Borge et al. 19891 and Ca (Trinder et al. 19971 



isotopes allowing for experimental tests of the predicted 
strength distribution. The GTGR is predicted to be ac- 
cessible even ioi N = Z nuclei above mass 64 ( Hamamoto 
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and Sagawa 19931. The experimental knowledge is still 



limited, but present data appear consistent with shell- 
model calculations ( Dossat et al. 2007 1 . For nuclei with 
N > Z the systematics of the Gamow- Teller strength 



Pinedo 20001. 



is given in (Batist et al. 2010 Langanke and Martinez 



There is a special interest in the nuclear structure 
around the doubly magic nucleus ^''^Sn. This is the last 
particle stable N = Z nucleus, the halflife has now been 



measured ( Bazin et al. 2008 1 for all of them. For nuclei 



with Z < 50 and 50 < iV allowed beta decay will mainly 
proceed via the Trggf2 ~^ ^97/2 transition and for nuclei 
approaching ^°°Sn all of the strength again appears to be 
accessible in beta decay. A comprehensive overview was 



given recently by Batist et al. (20101. 



4. Selected spectroscopic tools 

This subsection will present a few physics phenomena 
that can be employed to extract more detailed informa- 
tion on the states entering in beta-delayed particle emis- 
sion, namely recoil shifts, interference between levels and 
decays where individual levels are not resolved. 

To experimentally distinguish Fermi and Gamow- 
Teller transitions one may be guided by spin selection 
rules, but in general have to resort to beta recoil effects 



(Holstein 19741. The beta-neutrino angular correlation 



will give a significantly larger recoil shift in Fermi transi- 
tions than in Gamow- Teller transitions and can be stud- 



ied either as a function of beta-particle angle (Clifford 



et al. . 1989 1 or through measurement of the peak shape 



(Schardt and Riisager 19931. The size of the shift scales 



inversely with the mass number and is therefore easier to 
measure for light nuclei. It will depend on the spin se- 
quences in the de cay and has e.g . been used to determine 
the spin of Ar (iThaysen et al.\ 11999 1 . 



The level density of nuclei increases with excitation 
energy and with mass number. As it increases the lo- 
cal structure changes from rather regular to essentially 
chaotic, a transition well-studied theoretically but exper- 



imentally less understood in nuclei ( Weidenmiiller and 
[Mitchell] |2009| . In many nuclei around mass 100 the 
beta-delayed proton spectra will be dominated by un- 
resolved isolated resonances and fluctuation analysis is 
needed to extract information on the average spectral 



properties (Hansen et al. 1990 1 (see also Giovinazzo et al 



(20001 for later work around mass 70). The larger win- 
dows for beta-delayed protons in lighter nuclei close to 
the driplines will enable these studies to be continued to 
cases with different level density. In the decays with high- 
est Qec values one may reach excitation energies close 
to where Ericsson fluctuations have been observed in nu- 
clear reactions, i.e. the region where the level widths are 
larger than the average level distance. It will be experi- 
mentally challenging to look for such fluctuations in beta 



decay. 

Another aspect of spectroscopy at high level density 
is that "complete spectroscopy" will be very challenging 
to achieve, see the discussion in Hansen et al. (19901. A 



way of overcoming this challenge is the total absorption 
technique (Janas et al. 2005 Rubio et aZ.[|2005 l, where 
the aim is to measure the total emitted energy (apart 
from the emitted beta and neutrino particles) rather than 
the individual protons and gamma rays. This of course 
also holds for decays of neutron-rich nuclei where it, as 



demonstrated recently (Algora et al. 2010), is essential 



for a correct understanding of the decay heat in nuclear 
reactors. 

A final effect that can influence decay spectra signifl- 
cantly is interference due to overlapping levels of the same 
spin and parity. This will occur not only at high exci- 
tation energy, but also for otherwise well-resolved states 
whose tails overlap. Interference will be clearly promi- 
nent in light nuclei where broad states occur frequently, 
but it will certainly be an issue also for the broad states 
that can appear for heavier nuclei due to width collec- 
tivization once the level density is so large that levels 



start to overlap, see (Celardo et al. 2008 Zelevinsky 



1996 1 and references therein. Interference effects needs a 



more careful theoretical treatment, e.g. via the R-matrix 
formalism. The effects are often easy to identify once 
the statistics is sufficient and may range from slight dis- 
tortions, as in the beta-delayed proton spectrum from 
■^^Ar ( jSchardt and Riisager 19931, to considerable spec- 
tral modifi cations, as in the beta-de layed alpha spectrum 
from (Buchmann et al. 20071. However, interfer- 



ence effects are not always easily recognizable, as seen 
in the beta-delayed alpha -decays of ®B (Barker 19891, 
^^N (see subsectioiijl V . Dj) and to some extent also ^''N 
(Buchmann et al. 2009|, and when statistics is insuffi- 



cient spectral features arising from broad and interfering 
levels are easily misinterpreted as new weak transitions 
as demonstrated e.g. for the beta- delayed proton spec- 
trum from ^^Ne (Borge et al. 19881. This underlines the 



care that must be taken when interpreting decay spectra. 

C. /3~ delayed emission of one particle 

1. Occurrence of particle emission 



Figure 17 shows Qp and the beta-decay halflives along 



the neutron dripline. The values are theoretical estimates 
and will depend on the theoretical model chosen, in par- 
ticular on how the model predicts the nuclear shell struc- 



ture ( Sorlin and Porquet 2008 1 evolves. However, the fol- 
lowing general observations are most likely robust. (As 



shown by Moller et al. (19971 their theoretical halflives 



agree better with experimental value the larger the Q- 
value is.) 

The halflives for nuclei at the neutron dripline vary 
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FIG. 17 Q;3-value (upper figure) and beta-decay halflife 
(lower figure) as a function of neutron number A'^ for the light- 
est particle stable nucleus for a given A'^. The experimental 
dripline position is used for A'^ up to 30, all other values are 
taken from (iMoUer et a/.| 119971 120031). The dashed hne gives 



the Q-values from an estimate 
formula. 



based on the Weizsacker mass 



somewhat for the experimentally known ones with N < 
30, but are likely to mainly be in the range 1-3 ms once 
we get above N about 40. Deviations will be due to 
changes in the Q^-value rather than other structure ef- 
fects. The Q-values may be affected by shell structure, 
but decrease slowly towards higher masses. This overall 
trend is seen already from the simplest possible liquid 
drop formula, as also indicated in the figure. An even 
smoother dependence of halfiive with nucleon number is 
found in recent work based on the density functional ap- 



proach (Borzov et al. 2008), but it is clear that both 



Q-values and halflives vary much less for neutron-rich 
nuclei than for the proton-rich ones. 

Neutron emission will take place once it is energet- 
ically allowed and beta-delayed neutron emission will 
therefore be an important feature for neutron-rich nu- 
clei. The extent is illustrated in figure [TSj not only 
are beta-delayed multi-neutron decays energetically al- 
lowed shortly after beta-delayed one-neutron decay, the 



estimated beta strength distribution will soon give more 
than one emitted neutron on average per decay. As an 
example, for neutron dripline nuclei around mass 180 one 
expects more than 10 neutrons emitted in the decay chain 
towards beta-stability. Experimentally, we have today 
mainly reached this extended region of high P„-values 
for the light nuclei. 

Similar to what is observed for proton-rich nuclei the 
Coulomb barrier will limit significant beta-delayed alpha 
emission to the very light nuclei. However, the delayed 
emission of hydrogen isotopes, in particular deuterons 
and tritons, may also occur with small probabilities. 



Since their Q-values are limited, as seen from eq. ( 10 1 



the deuteron emission will be suppressed by three orders 
of magnitude at mass 100 whereas triton emission may 
still be possible to see up to mass 200. Their physics 
relevance will be discussed shortly. The major difference 
to the situation for proton-rich nuclei is therefore the 
prominent beta-delayed one-neutron and multi-neutron 
emission. 

The excitation energy of the Gamow- Teller Giant Res- 
onance, that carries the main part of the Gamow- Teller 
strength, is known to decrease linearly with respect to the 
Isobaric Analogue State as a function of (N — Z) /A 
ganke and Martinez-Pinedo| 



2000 Osterfeld 1992 



It 

was pointed out by [Sagawa et al.\ (|1993| that the GTGR 



Lan- 



for light neutron-rich nuclei (oxygen or below) even could 
move below the initial state so that a major part of the 
strength can be accessed in beta-decay. For heavier nu- 
clei, decays will take place to the tail of the GTGR. 



2. Decays in different mass regions 

To illustrate the present stage of the field, this section 
will present experimental results from several currently 
investigated mass regions, starting with the lightest nu- 
clei which is where the neutron dripline is reached and 



halo structures ( Jensen et al. 2004 1 have been studied 



Neutron halo nuclei must have low neutron separa- 
tion energy and have a "clustered" structure in the sense 
that the halo neutrons should decouple from the core to 
a large extent. It is obvious from eq. (10) that beta- 



delayed deuteron emission will be energetically favored 
in two-neutron halo nuclei, furthermore the component 
where the two halo neutrons decay to a deuteron (with 
the core as spectator) will give an important contribution 
to this decay mode. In fact, most theoretical calculations 
of the /3 d decay only includes decays of the halo neutrons 
directly to continuum deuteron states. The early work 



on this decay mode is reviewed in Nilsson et al. (2000). 
The decay has so far only been seen in ''He and ^^Li 
and the first experiments at ISOLDE have now been ex- 



tended at other laboratories both for ^He ( Anthony et al. 



2002 Raabe et al. 2009) and "Li (Raabe et aZ. "20081 



The branching ratio is for ^He now determined to be 
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FIG. 18 Nuclei with large beta-delayed neutron-emission probability are marked with an open square il the probability for 
emitting one or more neutrons is larger than 50% and with a filled square if the average number of emitted neutrons is larger 
than one. The P„ values are taken from experiment (Audi et al. 2003a Borge et al. 1997 [Yoneda et al. 20031 for A'^ < 20 and 



from Moller et al.\ ([2003 ) otherwise. The full lines indicate the line of beta-stability and the two driplines estimated from the 



Weizsacker mass formula and the broken lines the corresponding estimates for where beta-delayed one-, two- and three-neutron 
emission becomes energetically allowed. 



1.65(10) • 10^^ above a centre-of-mass energy of 525 keV. 
This very low value is understood to be due to cancella- 
tion in the matrix elements between contributions from 



small and large radii. The latest calculations (Tursunov 
et al. 2006a|b I reproduce both shape and intensity of 
the deuteron distribution, but it is not yet clear whether 
the theoretical and experimental maximum intensity po- 
sitions agree, so measurements at lower energy would still 
be valuable. For ^^Li the branching ratio is 1.30(13) TO^"* 
above a centre-of-mass energy of 200 keV and the spec- 
trum is again rather featureless ( Raabe et aZ.[|2008 l. The 



most recent theoretical calculations (Baye et al. 2006) 



give a qualitative agreement with data, but a real test of 
the theoretical understanding seems only possible once 
experimental data on the ^Li+d interaction at low en- 
ergy are available. 

The beta-delayed triton emission is again favoured at 
the neutron dripline and has been observed clearly in ^He 
and ^^Li and at the 10~^ level in ^''Be (Jeppesen et al 



nucleus is less well understood. Recent experiments on 
"Li (iMadurga et al.\ |2009| and ®He (iMianowski et al. 



20101 have confirmed the decay mode with new experi- 



mental procedures, see figure 19 but it seems that more 
experimental data is needed before one can determine e.g. 
whether the triton decays proceed through states in the 
daughter nucleus or, as the deuterons, directly to the con- 
tinuum. In the latter case the decay mode may depend 
on three- nucleon correlations in the decaying nucleus. 

The decays of A = 9 nuclei lead mainly to final states 
with two alpha particles and a nucleon. Complete kine- 
matics decay studies have been performed on these nuclei 
during the last decade and have resulted in the discov- 
ery of new decay branches and in spin determination of 



several intermediate levels (Prezado et al. 2003 20051. 



Strong Gamow- Teller branches in the mirror decays of 
^Li and go to states at around 12 MeV excitation en- 
ergy, but the deduced Bgt values are large and a factor 
4-5 larger for ^Li than for which is not possible to un- 



2002), but its relation to the structure of the emitting derstand from conventional theory ( |Kanada-En'yo 2010 
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FIG. 19 Beta-delayed multi-particle decays recorded with the optical tim e projection ch amber described in section |III.C.1| 

l) recorded so the incoming 



al. 



2007b 



The left panel shows beta-delayed three-proton emission from *^Fe (from (Miernik et 

track is not visible, the right panel the track of a *He ion entering from the right that after beta-decay breaks up into a triton 
(long weak track) , an alpha particle and an invisible neutron. 



Millener 20051. The reason for this is still unknown, the 



experimental strength ( Prezado et al. 2003 1 may perhaps 



investigated in the heavier neutron-rich nuclei can be 
mentioned the stability of the = 28 shell that has 



involve more than one level, and a proper theoretical in- 
vestigation of the three-body continuum may help resolve 
the puzzle. 

It would also be valuable to have calculations of the 
decays of the halo nuclei ^^Li and ^^Be that take the 
continuum into account explicitly. A more complete de- 
cay scheme for ^^Li is now available both at low (Mat- 



been probed by extensive halflife measurements (Grevy 



et al.\ |2004[ ) as well as the = 32 and possible = 34 
subshells probed in decays of Sc and Ti isotopes (iCraw- 



ford et al. 20101. The observed isotopic anomalies in 
some meteorites is known to depend on decay properties 
of very neutron-rich nuclei and motivated new measure- 



toon et aZ.[[2009|) and high ([Madurga et aL[|2008| excita- ( [Sorlin et aL[|2003[ ) isotopes. 



ments on the heavy Ar ( Weissman et al. . 2003 1 and Sc-Co 



tion energies in the daughter, but the deduced strength is 
still significantly less than that predicted by recent the- 



Recent experiments (Hosmer et al. 2010 Hosmer 



ory (Kanada-En'yo 20101, in particular it has still not 
been possible to experimentally check the above men- 
tioned prediction of the GTGR being placed below the 
initial state. A similar situation seems to be present for 
^^Be where the experimental decay strength distribution 



( Jeppesen et al. 2002 1 at high excitation energy is signif- 
icantly lower than theoretical predictions. A better de- 
termination of the beta-delayed neutron branches could 
alleviate the problem, but may not suffice to solve it. 

Much less is known about the decay of heavier dripline 
nuclei, but at least major beta-delayed neutron and / o r 
gam ma lines are known ou t to ^'' B (|Raimann et al. 1996 1, 



et al. 2005j Winger et al. . 2009 1 have succeeded in deter- 
mining halflives and P„ values for ^^Ni and nuclei around 
it. Apart from the interest in settling the properties of 
this doubly magic nucleus the information is also needed 
to fine-tune calculations of the astrophysical r-process in 
this mass range where there is sensitivity in particular 
to the halflife of ^®Ni itself (Hosmer et al. 



2010J . 



, At 

higher masses the nuclei participating in the r-process 
have been reached experimentally at iV = 82, see (Lan- 



ganke and Martinez-Pinedo 2003 PfeiS'er et al. 20011 



for more details. 
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Ozawa et al. . 1995), ' 


2010 


and (Reed et al. 



N ( Sumithrarachchi et al. 3. Selected spectroscopic tools 



1999) and halflives and P„ 



values are known for the heavier B, C and N isotopes 



(Yoneda et al. 2003). In the region above oxygen the 
halflives are not known for the most neutron-rich iso- 
topes of any element. The major decay branches are 
established for nuclei at a similar distance fro m the line 
of stabiUt y, e.g. for ^^Ne (|Tripathi et a?.||2006| and ^^Na 
(i Nummela et aL[ |200l] |Radivojevic et aL[|2002| ), but the 



Beta-recoil effects may play a role for beta-delayed neu- 
trons similar to that discussed above for beta-delayed 
protons, but has only been explored in a few cases such 
as ^Li (Nyman et al. 1990). If the nucleus recoiling 



dripline from here on is significantly further out (see sec- 
tion |ll]). 

Among the different physics questions that have been 



from neutron emission emits a gamma ray the latter 
will also be Doppler broadened provided the gamma- 
emitting state is sufficiently short-lived. This has been 
used to improve the decay scheme of ^^Li (Fynbo et al. 



2004 Mattoon et al. 2009), and is a valuable way to 



cross-check results from the sometimes complex neutron 
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FIG. 20 (Color online) The branching r atio for beta-delayed alpha decay of ^^N (filled triangles) is shown as a function of 



the total energy (Hyldegaard et al. 20101. The solid line is a fit to the feeding to and 2"*" states and does not include the 



contribution to the 1^ state at 12.7 MeV. The filled circles (open squares) give the contribution from decays that do (do not) 
proceed through the *Be ground state, the dashed lines are the corresponding fits. See the text for details. 



spectra and neutron-gamma coincidence measurements, 


et al. 


1983), 


see (Hirayama et al. 12005*) and references therein for 
^^Li. Analogous neutron-gamma experiments have al- 


few years ago 
19 1 and only 


ready been performed e.g. for ^^N (Li et al. 20091, '^'^Mg 


et al. 


2011al. 



45t 



and ^'^Al ("Angelique et aZ.I|2006l 



Fe (Miernik et al. 2007b ) (see figure 



19 1 and only recently also reported in^^^Cr (Pomorski 



The analysis of decays through regions of high level 
density proceeds similarly to the case for the proton rich 
nuclei, except that the experimental challenges are higher 
due to the neutron detection. Fluctuation analysis will 
again be an important tool in order to extract reliable in- 
terpretations from beta-delayed neutron spectra (Hardy' 



et al. 



19781. For lighter nuclei where the level density 
is smaller one should in principle in the analysis of neu- 
tron spectra worry about exact lineshapes, interference 
effects etc. as was the case for the corresponding delayed 
proton spectra. However, most experiments presently de- 
termine neutron energies through time-of-flight and as- 
sume (at least implicitly) that their resolution will smear 
out such effects so that peaks in the spectra can be fit- 
ted with Gaussians. This practice could lead to wrong 
assignments. 



D. Beta delayed emission of several light particles 

Apart from decays through *Be and states in ^^C above 
the triple-a threshold (and a few weak transitions involv- 
ing an a particle and a nucleon such as occurring in the 
decay of ^^Ne (Chow et al. 2002)) beta-delayed emis- 
sion of several light particles involves only nucleons. The 
first beta-delayed multi-nucleon decays, /32n and /33n, 



were discovered about 30 years ago (Azuma et al. 1980 



1979). The /32p process followed shortly after (Cable 



As shown in figure[T8]beta-delayed multi-neutron emis- 
sion will become dominant in the decays of very neutron- 
rich nuclei, whereas the other processes only occur with 
small to moderate intensity (with the exception of the 
A = 8, 9 decays). Somewhat ironically, the multi-neutron 
process are the least studied ones, partly for experimental 
reasons due to the difficulty of neutron detection, partly 
due to the quite few cases of beta-delayed multi-neutron 
emission known today. The one case, ^^B, where beta- 



delayed four neutron emission has been reported ( Dufour 



et al. 1988 ) needs to be confirmed since other multi- 



neutron branches reported in the same work has since 



been shown to be too large (Bergmann et al. 1999). 



The question of the particle emission mechanism is of 
prime importance. Although, as discussed in more detail 
in section |VII[ calculations of multi-particle final states 
in principle are becoming feasible now, it is still of in- 
terest to know whether simpler decay mechanisms, such 
as sequential decay, can describe a process or whether 
break-up directly into multi-particle continuum states 
takes place. Currently the only experimental information 
on beta-delayed multi-neutron emission comes from sin- 



gle neutron spectra (Azuma et al. 19791. There is more 



and Borge| |2008|. 

that of ""^Ar (iFynbo et al. 



knowledge on /32p decays, as recently reviewed in (Blank 
The most th oroughly studied case, 
seems to display only 



20001 



sequential emission and other /32p cases are consistent 
with this. Turning to other beta-del ayed multi-particles 
modes, the /3ap decays from ^^Ne (Chow et al. 20021 
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FIG. 21 (Color online) The beta-delayed alpha decay data for N shown in figure 20 are displayed corrected for the beta-decay 



phase space factor and alpha particle penetrability factors (from ( |Hyldegaard[ |2010[ )). The total fit is divided into contributions 



from 0+ and 2"' 
energy. 



states in C. Note the clear interference between 0^ states at low energy and the enhanced decay rate at high 



could also be analyzed assuming only sequential decays. 
For the light nuclei, a recent study ( [Madurga et al. 2008 ) 
of ^^Li indicated that most of the decays in the three- 
body (n-l-a+^He) and the five-body (2a + 3n) channels 
are sequential and to a large extent proceeds through var- 
ious He isotopes; a smaller direct break-up component 
is, however, still possible. Finally, for the A — 9 decays 
mentioned above in section |IV.C| most decay branches 
can be approximately described in a sequential picture. 



but there are indications ( Prezado et al. 2005 1 that there 
is a direct break-up component from the 5/2~ level at 
2.43 MeV in ^Be. The hypothetical sequential branches 
through ^HCgs and ^Be(2+), which give energy distri- 
butions that for the higher-lying broader states can be 
distinguished and seem to be observed, would give over- 
lapping energy distributions for the 2.43 MeV level. The 
sequential picture anyway does not really make sense in 
this situation as discussed in section |VlT] For the specific 
case of the 2.43 MeV level the break-up mechanism has 



fits. 

It would be interesting to have more detailed data 
on the decay mechanism for /33p decays. On a longer 
timescale it is without doubt the multi-neutron detection 
capabilities that constitute the key challenge for future 
progress in this field. 



1. The case of yl = 12 

The complications that may arise in beta-delayed de- 
cays can be illustrated with the case of ^^N (and ^^B) 
whose decay into the 3a continuum has recently been 
studied in detail fDiget et al. 



2009 Hyldegaard et al.\ 



2010 2009 1 motivated by the importance of this contin- 



uum for the astrophysical triple-a process ( Fynbo et al. 
2005 1 . The decay goes through narrow states as 



been investigated in three-body calculations (Alvarez 
Rodriguez et al. 2008a 



well as through several 0+ and 2+ states that all couple 
strongly to the continuum. The experimental spectrum 
is shown in figure [20] as a function of the excitation en- 



where the experimental energy ergy in ^^C. One can experimentally identify decays that 



and angular distributions could be reproduced. 

The analysis of such data often makes use of the R- 
matrix formalism since this allows for level parameters to 
be fitted to experiment. The adaptation of the formalism 
to beta decay is described e.g. in [Barker and Warburton] 
( 1988 ). It is a priori applicable only for two-body decays, 



but has been employed in practice also for sequential de- 



cays due to the lack of better approaches. Robson ( 1975 ) 



proceed through the narrow ground state in "Be, these 
decays are also marked in the figure and obviously cor- 
respond to sequential break-up. Apart from the Hoyle 
state at 7.65 MeV all other states in the fit are broad. 
The resulting interference may be easier to see in figure 
21 where the data have been corrected for the beta phase- 
space factor and the alpha particle penetrabilities. The 
interference actually also involves the upper tail of the 



has shown how to formally make sense of extensions of 
R-matrix to multi-particle situations, but this has not 
been implemented in data analysis. One can therefore 
not rely fully on results derived from present R-matrix 



Hoyle state (its "ghost", cf. Barker and Treacy (1962l) 
that owes its narrow width to a small value for the pen- 
etrability. 

The increase in strength at higher energies seen for the 
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2+ states is very hard to reconcile with the sum rule, eq. 
([9]); if the decays are a ssumed to proceed throu gh levels 



in C. This was taken (Hyldegaard et al. 2010 1 as an in- 



dication that some of the ^•^N decays proceed directly to 
continuum states. As mentioned when discussing figure 
[13] the standard assumption is that beta decays proceed 
through states in the emitter rather than directly into 
continuum states. Corresponding transitions directly to 
the continuum have been known for a long time for strong 
and electromagnetic processes (direct reactions and di- 
rect radiative capture, respectively), but are not gener- 
ally recognized to occur also in weak decays. However, 
as mentioned in section IIV.CI it is the most natural ex- 
planation also for /?d decays. 



E. Beta delayed fission 

Beta-delayed emission of particles heavier than alpha 
particles has only been seen as beta-delayed fission. An 
overview of this phenomenon with references to the early 



work can be found through 


Hall and Hoffman ( 


19921; 


Kuznetsov and Skobelev 


(1999 


I; and 


Shaughnessy et al. 



(2002 1 . The probability for such decays depends both 



on the beta strength at high excitation energy and on 



the fission barriers (Moller et al. 20091, and the decay 



mode may therefore provide experimental information on 
fission in regions with high Qjj values that is hard to ob- 
tain otherwise. Such information will enhance our un- 
derstanding of the fission process and can help to deter- 



mine better the role of fission in the r-process ( Martinez 



Pinedo et al. 2007 1 



Much of the recent activity has been driven by the 
continuous developments in radioactive beam production 
capabilities and has focussed on EC delayed fission in the 
light mass region. Experiments have been carried out on 



At at GSI and on at ISOLDE (Andreyev et al 



2010). The latter experiment showed a surprising asym- 
metry in the mass distribution of the fragments. More 
detailed information should become available in the com- 
ing years. 



V. SINGLE-PROTON RADIOACTIVITY 
A. Introduction 

The proton radioactivity is the process occurring in 
odd-Z nuclei located beyond the proton-drip line. Due 
to the potential barrier (Coulomb and centrifugal) the 
emission of a proton from an unbound nucleus success- 
fully competes with other forms of decays (/3+, a) only 
when the Qp value for the decay is sufficiently large, see 
Eq.p^ Th e proton radioactivity was discovered by |Jack- 



son et al. (19701 who observed protons emitted from 
an isomeric state in ^^Co at the excitation energy of 
3.2 MeV. First observation of the ground-state proton 
radioactivity was reported 12 years later by Hofmann] 
et al\ (|1982| for ^^iLu and by [Klepper et ol] (1982) for 
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Tm. Presently, more then 40 proton emitters (from 
g^^I to gg^Bi), including emission from long-lived isomeric 
states, have been established experimentally. Six of them 
("iEu,i4i™Ho,"i9^Ho,"4Tm, i45xm,i46Tm) have tran- 
sitions (so called fine structure) to the excited states in 
the respective daughter nuclei. 

The importance of proton radioactivity follows from 
the fact that the knowledge of the decay energy and 
the half-life (width) , combined with the relatively simple 
model of the potential barrier penetration, yields infor- 
mation on the nuclear wave function. Thus, relatively 
simple observables provide constrains on nuclear mod- 
els for exotic nuclei, located beyond the proton drip-line. 
Since nuclear structure information is usually interpreted 
with help of the shell- model, it is convenient to divide 
proton emitters into two groups: those of the combined 
seniority one or two (s < 2) and others with the com- 
bined seniority larger then two (s > 2). The combined 
seniority is defined as the number of unpaired nucleons 
(protons and neutrons). In the first case, apart from 
the odd-proton, and possibly an odd-neutron, no proton 
and neutron pairs are broken. Such decays are typical 
for ground states and for low lying isomers. In the s = 1 
case the odd proton can be pictured as moving in a single- 
particle orbital in the nuclear potential of the even-even 
daughter nucleus, while the s ~ 2 case corresponds to 
even-odd daughter with an odd neutron acting as a spec- 
tator. The majority of known proton emitters belong to 
this category, they will be discussed in Sec. |V.B[ The 
situation of s > 2 corresponds to the proton emission 
from highly excited isomers having multiparticle charac- 
ter, which involves breaking additional pairs of protons 
or neutrons. In Sec. IV. CI we discuss a few known cases. 

In addition to the information extracted directly from 
proton emission observables, nuclear structure informa- 
tion has been gathered by using the emitted protons as 



a tag in the recoil decay tagging (RDT) studies (Sewery- 



niak et al. 2007b 2001 Yu et al. 19981 



An overview of nuclear structure studies at the pro- 
ton drip-line by means of proton radioactivity was given 



31 



recently by Blank and Borge (2008). More detailed dis- 2. Probability of proton emission 



cussion of proton radioactivity was given by [Woods and 



Davis (19971 and Hofmann (19951. The work of Son- 



zogni ( 2002 1 contains a compilation of results on proton 



emitters known in 2001. 



1. Fundamentals 



Relatively simple calculations of the proton emission 
lifetimes are based on the result obtained by |Gurvitz and] 
Kalbermann ( 1987 1 who analyzed the decay widths and 



The necessary condition for a nuclide to decay by pro- 
ton radioactivity is a positive decay energy Qp defined as 
the difference between binding energies (Eq. [T2| ) of the 
parent and the daughter atoms: 



shifts of quasistationary states in the quantum mechan- 
ical two-potential approach. By investigating the quasi- 
classical limit they provided simple formulae which are 
similar to, but more general than, those achieved in the 



framework of WKB approximations (Brink et al. 19831 



In this approach the width of the proton-emitting state 
is given by: 



Qp = B{N, Z-1)- B{N, Z) = -Sp 



(12) 



N 

= Sp — exp[— 2 / k(r) dr] 

V Jr2 



To separate the contribution from the atomic electrons, 
the decay energy is expressed in the form: 



k{r) = {2^,\Q;-^-V{r)\) 



1/2 



(20) 
(21) 



?p — Qp 



ES. 



(13) 



where the normalization factor N has to satisfy the equa- 
tion: 



where Q™'^ is the nuclear part of the decay energy. It is 
determined by the nuclear masses: 



f''^ dr 



(22) 



qnuc ^ M^-^^iN, Z) - M""'''{N, Z - 1) - TOp , (14) 5^ jg ^]^g spectroscopic factor described later and fj. is the 

reduced mass of the proton and the daughter nucleus. 
Integration limits are the classical turning points, de- 
fined by V{ri) = Qp"'', where V{r) is the radial part of 
the nucleus-proton potential, see Fig. [22] 

To simplify calculations, some authors replace the fac- 



where nip is the proton mass. The ES is the electron 
screening correction defined as the difference between to- 
tal electron binding energies in the parent and the daugh- 
ter nuclides: 



ES = Be{N, Z) ~ Be[N, Z-1) 



(15) tor ^ by the so called frequency of assaults factor v ( Hof 



In the above the electron binding energy in the hydrogen 
atom has been neglected. The value of the screening cor- 
rection ES can be calculated from the tabulated electron 



mann 



1996 ) calculated for the case of an s-wave proton 



leaving the square well plus Coulomb potential (Bethe 
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binding energies (Huang et al. 19761, or estimated by a 
simple formula: 



ES = 0.49 + 0.0144 • Z^-^ keV. 



(16) 



^i^RlJ{Z-l)e:^/R,-Q^^- 



(23) 



The accuracy of this parametrization is better than 0.5% 
for 42 < Z < 75 and drops to 1.6% for Z=83. No influ- 
ence of the neutron number (isotopic effect) is taken here 
into account. 

The decay energy is shared between the proton Ep and 
the recoiling atom. Therefore, the measured kinetic en- 
ergy of the emitted proton is given by: 

(Af (iV, Z + nie) 



where Rc is the channel radius Rc — r„„c — ^oi^^ 1)^/^ 
with ro — 1.21 fm. For example, in the decay of ^^^Lu i' 
equals 4.1 • 10^^ s^^ 

Then, the proton emission decay constant is calculated 
as: 



\ = SpV exp[— 2 / k{r)dr] 

Jr2 



(24) 



E„ — 



nip + M{N,Z - 1) + TO, 



P ■ 



(17) 



In addition, the angular momentum as well as parity con- 
servation laws have to be satisfied: 



The potential V{r) is taken as a superposition of nu- 
clear {Vn), coulomb (Vc), centrifugal Vl, and spin-orbit 
Vis terms. The nuclear part is usually described by the 
Woods-Saxon form with various parametrization: 



TTf 



If + {1 + S) 

{-ly 



(18) 
(19) 



VN{r) = -Vr 



1 



1 + cxp^- 



(25) 



where, li and If are spins of the initial and final nuclear 
states respectively, I is the angular momentum of the 
emitted proton, s is the spin of the proton, tt^, tt^ are 
parities of the initial and final states, respectively. 



Comparison between different parametrizations was done 



in the work of Ferreira et al. (20021. Although not sup- 



ported by the work of Ferreira et al. ( 2002 1 , the most 



frequently used parametrization is that of Becchetti and| 
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FIG. 22 Schematic view of the radial part of the nucleus- 
proton potential. The classical turning points for a particle 
with energy Er are marked. The nuclear contribution turns 
to zero around r = Tnuc- 



|Greenleesj (19691. Detailed potential descriptions and 



calculations for some of the emitters can be found in 



Aber^ 


; et al. 


(1995 


). 



Buck et al. 



(19921; and 



Hofmann 



As an example, the half- life for the proton emission 
from ^^^Lu as a function of the decay energy, for three val- 
ues of the orbital angular momentum is shown in Fig. [23] 
The characteristic strong dependence on the available en- 
ergy and the angular momentum is clearly seen. 



3. Spectroscopic factor Sp 

The spectroscopic factor Sp is a measure of the single- 
particle purity of the initial wave function. Within the 
BCS theory the spectroscopic factor is given by S'*'* — uj, 
where the vacancy factor is the probability that the 
spherical shell-model orbital with (n, l,j) quantum num- 
bers is empty in the daughter nucleus. For some pro- 
ton emitters the factors can be found in the work of 



Aberg et al. (19971. This theoretical value is compared 



with the experimental value Sp^^ derived as the ratio of 
the measured partial decay constant and the calculated 
one assuming Sp = 1. The agreement between the two 
values indicates that the correct assumption about the 
initial wave function has been taken. For example, in the 
case of ^^^Lu Sp'^P = 0.5 and 5*^ = 0.54 assuming the 
proton was emitted from the [7r/iii/2]ll/2~ state. The 
good agreement supports such an interpretation. In con- 



trast, for the case of "Sxni S^'^p 



0.48 and S: 



'th 



0.65, 



which suggests that the spherical potential used in the 
calculation may not be applicable. Indeed, the deformed 



et al. 2003 


Seweryniak et al. 


(Arumugam et al. 


20081. 
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FIG. 23 The half-life for the proton emission as a function of 
nuclear decay energy Qp"^ and the orbital angular momentum 
carried away by the proton. Calculations were done for the 
case of ^^^Lu with Sp = 0.54. The measured values of the 
decay energy and the half-life, indicated by the black square, 
suggest the transition with L = 5. 



4. Models of proton emission 

The usefulness as well as limitations of simple models, 
introduced in previous sections, may be illustrated with 
the example of ^^^Tm. In this nucleus two proton tran- 



et al. 


2003 


and E 





1.40 MeV 
half-lives are 3.4 ^s and 32 /is 



sitions from the same state have been observed (Karny 
The proton energies are Ep = 1.73 MeV 
while the corresponding partial 
The first transition is in- 
terpreted as a decay to the 0+ ground state of -'^'''^Er, 
while the second goes to the first excited 2+ state in 
this nucleus. In the frame of the spherical quasi-classical 
approach we may assume that the first transition orig- 
inates from the 7r/iii/2 orbital (/ = 5) while the second 
from the 77/7/2 = 3) component of the initial wave 
function. The calculations, including theoretical spectro- 
scopic factors (w^(/iii/2) = 0.647 and u'^{fj/2) = 0.985) 
yield the partial half-lives of 2.29/j,s and 1.28/is, for the 
two transitions, respectively. Thus, by comparing with 
the experimental values, we may conclude that the emit- 
ter wave function consists of 67%= • 100% of the 
I = 5, 7r/iii/2 ^ 0+ and 3.7% oi I ~ 3, 77/7/2 ^ 2+ compo- 
nents. The remaining 29% of the wave function does not 
participate in proton emission and therefore can not be 
determined within this simple model. The more elabo- 
rate coupled channel model which takes into account the 



dynamic deformation (Hagino 2001 Karny et al. 20031 
yields values of 56% for 7r/iii/2 (81 0+ state and ~ 3% for 
7r/7/2 ^ 2+ state. 

Although simple models are useful for the first order 
approximations, they cannot be expected to yield cor- 
rect results for highly deformed nuclides. For those cases 
more elaborate theoretical approaches have to be ap- 
plied. Examples of such approaches are given by Aberg 



2007a I and theoretically et al. (19971; Esbensen and Davids (20011; Fiorin et al 



(2003); and Kruppa and Nazarewicz (20041. Calcula 
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TABLE III Summary of literature data for s = 1 proton emitters. Empty place means no data available. For isotopes 
with ground and isomeric state emission a combined literature is given. In 'Reference' column additional letter E - refers to 
experimental papers, while letter T - points to the theoretical papers where properties of the referred nuclei are explicitly 
calculated. 



Emitter 



55 "-^S 



'La 



Cross sec. 



40fih 



0.5/ib 



30^ib 



240nb 



0.3nb 



9nb 



90nb 



6nb 



13nb 



1.4/ib 



lOnb 



Ep(MeV) 



0.8126(40) 



0.807(7) 



0.806(5) 



0.882(10) 



1.020(15) 



0.932(7 
811(7 



s 



1.179(7) 



1086(10) 



1169(8) 
968(8)" 
1235(9) 
1031(11)" 



1700(16^ 
1430(25 



Qg""'(MeV) 



0.829(4) 



0.815(7) 



0.9771(37) 



0.813(5) 



0.900(10) 



1.028(15) 



1.188(7) 



Ti/2 



93.5(3)/is 



0.5(l)ms 



16.7(7)/is 



24(3)ms 



lOl^ ms 



0.90 



+0.49 



17.8(19)r 



0.94 



+0.33^ 



6(3)ms 



4.1(1) 

7.4(3)^is 
7.4(3)^is 



Aug. mom. 



1=2 



1=2 



1=2 



1=2 



=4 or 5 



1=2 



1=2 
1=2 



1=5 



1=5 



1=5 
1=3 
1=0 
1=2 



1=5 
1=2'' 



References 



Facstcrman et al. 



Gillitzer et al. 



Heine et al. 



jl984 l (E) 



(19871 (E) 



(1991|_(E) 



Barmore et al. 



Sellin et al. 



(20001 (T) 



( 1993 1 



Mazzo cchi et al. (2007) (E) 
Page et (1994; 



(E) 



Ferreira and Maglione (20011 (T) 



Faesterman et al. ( 1984 1 (E) 



Gillitzer et al. 



Page et al. ( 1994 1 (E) 



(19871 (E) 



Batchelder et al. 



Maglione et al. { 1998 1 (T) 
Barmore et^[200o| (T) 



(19981 (E) 



Soramel et al. (2001) (E) 
Mahmud et aZ.| |20Cri|) (E) 



Bogdanov et al.\ ( 1973 1 (E) 



Robinson et al. ( 2005 I (E) 



Mahmud et al. ( 2002 1 (E) 
^^^^jj — . \ 



Davids et al. (2004) (E) 



Maglione and Ferreira ( 2000 1 (T) 



Davids et al. 



Sonzogni et al. 



1998HE) 
( T999I (E) 



Maglione et al. 



jl999 l (T) 
Davids and Esbensen] ( 2000| (T) 



Kruppa et al. 



(20001 (T) 



Esbensen and Davids 



Ferreira et al. 



I2OO2RT} 



(20011 (T) 



Ferreira and Maglione ( 2005 1 (T) 



Woods et a/.|(|2004[) (E) 

Rykaczewski et al"ll%m\ (E) 
Maglione and Ferreira^ (^2002J (T) 



Ferreira and Maglione 



Davids et al. ( 1998 1 (E) 
Ryk aczewski et oL J 1999| (E) 



Maglione et al. 



Barmore et al. 



(20011 (T) 



(|1999| (T) 
( 2000| (T) 



Esbensen and Davids' (2001^) (T) 
Seweryniak et a/. ( 2001) (E) 



Karny et al. ( 20081 (E) 
Arumugam et al. (2009) 



(T) 



Grzywacz et al. X2005^ (E) 
Bingham et al. (20051 (E) 



^ transitions to the excited states 

^ evidence of the transition to the excited state based on two counts 
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TABLE IV Summary of literature data for s = 1 proton emitters. Continuation of Table [Till 



Emitter 



?Tm 



73 -L<* 
156mrp^ 



me 



i^s^^Re 



Cross sec. 



O.S^ib 



Ifih 



30Atb 



3/ib 



70Atb 



50nb 



20nb 



l^ib 



150nb 



Ep(MeV) 



1728(10), 
1398(10) 



1191(1) 
1016(4)" 



1120(1)" 



1110.8(3.9) 



1261(4) 
1286(6) 



1232.9(2.0) 
1310(10) 



1444(15) 



1007(5) 
1108(8) 



927(7) 



1805(20) 



1261(6) 



1192(2) 
1315(7) 



(MeV) 



1071.4(3.3) 
1139.3(5.3) 



Ti/2 



3.1(2)^is 



66ms 



200ms 



560(40)ms 
360(40)ms 



46(5)ms 



80(2)ms 
16(l)^is 



144(24)ms 
375(54)ms 



21(4)^is 



0.79(16)ms 



440(l)/is 
14.7(3)ms 



Ang.mom. References 



1= 



Livingston et oLj 1993 1 (E) 



Batchelder et al. 



Rykaczewski etoL (2001b I (E) 



Karny et a/.,()2003l (E) 



Seweryniak et al. 



Seweryniak et al. 



Arumugam et 



Rykaczewski (20021 (E) 



1998ME) 



(2007a 
2008, 



2005a | (E) 

I (E) 
(T) 



Tantawy et al. (20061 (E) 



Rykaczewski et aVhmih) (E) 
Rykaczewski et a/.^ ( ^2001a^ (E) 



Ginter et al. 



(20031 (E) 



Seweryniak et al. ( 2005a I (E) 



Klepper et aZ.Jj 19821 ) (E) 
Larsson et a l. ( 1983^ (E) 



Sellin et al. 



Toth et al. 



(|1993| (E) 
( 1993t (E) 



Seweryniak et al. ( 1997 1 (E) 
Seweryniak et al. (2005a) (E) 



Sellin et 



g/j (il993|) (E) 



Woods et aZ. ( 1993 ) (E) 



Ginter et al. 



( 20001 (E ) 



Ferreira and Ma glione' ('2001 ) (T) 



Maglione and Ferreira^ ( ^2002| (T) 



Ginter et al. 



(2003 i(E) 



Robinson et al. (2003)(E) 



Yu et al. 



Hofmann et a/.|()1982| ) (E) 



Sellin et al. 



(1993) (E) 



Bingham et al. 



( 1998 | (E) 

(fl999|) (E) 



Ferreira andllaglione| ( 200o| (T) 



Uusitalo et al. J 1999 1 (E) 



Page et al ( 2007 ) (E) 



Joss et aV (12006 1 (E) 



Page et al. (1992 ) (E) 
Livingston ef al. ( 1993| (E) 
Page et ffl/.'('l996') (E) 



Page et a/., (1996) (E) 



Irvine et al. (1997) (E) 



Joss et al. ( 2006 1 (E) 



Page et al. ( 2007 ) (E) 



Page et al. 



1992 



(1996 



Page et al. 
Hagino '( '2"00 1 ) (T) 



(E) 
(E) 



Irvine et a/.|(,1997| (E) 
Lagergr en et a l. " 2006 ) (E) 



Hagino ( 20011 (T ) 



Arumugam et aV ^2007) (T) 



see tlie discussion in | |Page et aZ.[ |2007[ ) about contradictory results from the work of | [Uusitalo et al.^ jl999| ) 



35 



TABLE V Summary of literature data for s = 1 proton emitters. Continuation of Table [Till 



Emitter 


Cross sec. 


£L/pt^iviev ) 




J- 1/2 


Ang.moiQ. 


References 


164m T„ 

77 


n/a 


1807(14) 




58l«MS 


1=5 


Kettunen et al 


(|2001 


(E) 


Mahmud et al. 


'.'I Ml.-' 


(E) 


^77 ""Ir 


0.2^b 


1707(7) 




0.30(6)ms 


1=5 


Davids et al.^ ( 


L997 


\ (E) 


166 T 

77 ir 
7*7 Ir 




1145(8) 
1316(8) 




10.5(2.2)ms 
15.1(9) 


1=2 
1=5 


Davids et al. { 


L997 


1 (E) 






167t 

77 ir 

167mT 
77 ir 


lOfih 


1064(6) 
1238(7) 




35.2(20)ms 
30.0(6)ms 


1=0 
1=5 


Davids et al. { 


L997 


1(E) 


Davids and Esbensen 


(20001 (T) 


Scholey et al, ( 


200f 


).(E) 


1 7n * 
JgAu 

170m A 




1463(12) 
1735(9) 




286l4jl^is 

^+50 

617l4o^ims 


1=2 
1=5 


Mahmud et al. 


(j20( 


m 


(E) 


Kettunen et al 


(2004 


(E) 


"V^Au 


2^ib 


1437(12) 
1694(6) 




22t^^s 
1.09(3)ms 


1=0 
1=5 


Davids et 


al. ( 


L997 


1 (E) 


Poli et al. 


Il99 


91 (E) 


Back et al 


[(20 


031 (E) 


(E) 


Kettunen 


et al 


TT2OO4 


176rpl 

81 




1258(18) 




5.2t?°ms 


1=0 


Kettunen et al 


(2004) (E) 


81 
177mrp] 
81 


30nb 


1156(20) 
1958(10) 




18(5)ms 
230(40)^is 


1=0 
1=5 


Poli et al.\ 


(jl99 


9() (E) 


Davids et al}\'. 


>001 


04 


3) 
(E) 


Kettunen et al 


(20 


185T5; 
83lJl 


eonil^] 


1594(9) 




60(4)Ats 


1=0 


Davids et al. | 


L996| (E) 


Poli et al. 


(200 


11 (E) 




Andreyev 


et al 


(2004 


(E) 
(E) 


Andreyev et al 


. (2005 


Arumugam et 


al. C 


MO 


71 (T) 



6-iOnb for 3n evaporation channel (Andreyev:2004) 



tions with tri-axially deformed potential are considered 1. Odd-mass, s = 1 proton emitters 



by Arumugam et al. (2009); Davids and Esbensen (2004); 



ajid Kruppa and Nazarewicz (2004). An extended de 



scription of the theoretical models used for proton emis- 
sion can be found in chapters 4,5, and 6 of the Lecture 



et al. (2006b). 



notes in Physics by Delion (2010), as well as in Delion 



B. Seniority s < 2 proton emitters 

Seniority s < 2 proton emitters are the "classical" pro- 
ton emitters in which an unpaired proton leaves the nu- 
cleus from the ground state or the isomeric state. The 
isomeric state has to be low enough to allow only for 
unpaired particle excitation. The half-life of those emit- 
ters span from T1/2 = l-^tol /is for "''Tm to 0.560(40)s 
for ^^^Tm. They were all but one produced in fusion- 
evaporation reactions with the exit channel containing a 
proton and from 1 to 6 neutrons. ^^^Bi was also pro- 
duced in the 3n evaporation channel ( Andreyev et al\ 
2004). Typical cross sections range from 0.3 nb (Ip, 6n) 
foP^iPr to 70 /ib (Ip, 2n) for ^^iLu. Combined proper- 
ties of "classical" proton emitters are presented in table 



As liientioned above, due to the strong dependence of 
the half-life on the Qp-value and on the angular momen- 
tum of the emitted proton (see figure 23 ) , proton emis- 



sion is a valuable tool for nuclear structure study beyond 
the proton drip line. Measured proton energy and decay 
half-life in most of the cases directly point to the config- 
uration of the decaying orbital. This is especially true 
for the odd-mass seniority s — 1 cases where the emit- 
ting nuclei can be described as even-even 0+ core cou- 
pled to the unpaired proton. With the O"*" ground state 
of the daughter nuclei, establishing the angular momen- 
tum leaves only two possibilities for the total angular 
momentum j = Z ± 1/2 of the emitting state. This value 
can then be used to calculate the components of the nu- 
clear wave function of the emitter. Depending on the 
nuclear shape either spherical single-particle orbitals or 
Nilsson type deformed orbitals can be used. Evolution of 
proton emitting states starts with {Z = 53, / = 2) 
proton radioactivity. In this region just above Z — 50 the 
7rd5/2 and TTg'j/2 orbitals are close to each other, neverthe- 
less due to the lower / of the former configuration {TTd^/2) 
most of the emitted protons carry away two units of angu- 
lar momentum. There is though one case discussed in the 
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literature namely ^^^Pr, for which Robinson et al. (2005) 
suggests two possible configurations 3/2+ [422] from the 
spherical TTgr/2 orbital and 3/2~[541] from the spherical 
"■^11/2 orbital. While in case of a high deformation the 
3/2+ [422] state can contain admixtures of the spherical 
77^3/2 and 7rc?5/2 orbitals, making the I = 2 proton transi- 
tion possible, the negative parity 3/2~[541] state cannot 
contain an Z = 2 contribution. In the Delion et al. (2006b ) 
systematics I = 2 is assigned to the proton emission from 
i2ipj, jjiaking it "compatible" with other proton emitters 
in this region. 

Although for proton emitters above Z — 64 three or- 
bitals ((^3/2, ^11/2 and S1/2) should be considered as pro- 
ton emitters, an interesting phenomenon occurs for ^■^^Tb 
(Z = 65) and for "iRo {Z = 67) due to the high de- 
formation. The 7/2~ ground states of both emitters are 
dominated by the [523] component from the /111/2 orbital, 
but the proton emission is driven by a small admixture 
of the /7/2 orbital to the wave function. The small I = 3 
component (<2% in case of ^^^Ho) wins o ver the larger 
= 5 component (>78% in case of ^^^Ho) (Karny et al. 



2008 ). The presence of the 1^/7/2 component in the wave 
function of nuclei in this region is also confirmed by the 
analysis of ^''^Tm fine structure data. The 9.6% proton 
branching to the first excited 2+ level in ^^"^Er can only 
be explained by the presence of the I = 3 component in 
the wave function ( Karny et al. 2003 ) . The small energy 



difference between /iii/2 and ^3/2 related levels manifests 
itself through the presence of I = 2 proton emitting low 
lying isomers in -'^'*^™Tm and ^^^""Lu. 

Odd-mass proton emitters above Z ^ 72 are char- 
acterized by the presence of two proton emitting states 
1/2+ ground state and the ll/ 2~ isomer. T here are two 
cases ^^^Re ( |Joss et al.\ |2006| and ^^^^Xa (I Page et al. 
|2007[ ) where only one I = 5 proton emission has been 
observed. In the case of ^^^Re, the expected half-life for 
I = emission is below l^s, which was beyond the ca- 
pability of applied experimental technique. The ^^^Ta 
has been observed as a second generation decay after 
^^^Re ion implantation and its subsequent a decay to 
^^^Ta. Although the expected half-life of Z = proton 
emission from ^^^Ta is long enough to be observed, the 
combination of the low production cross section for the 
1/2+ state with a small detection efficiency for the sec- 
ond generation decays, may explain the non-observation 
of the I = proton channel. Further studies have still to 
confirm that observed decays in both cases do originate 
from the 11/2^ isomeric states. The heaviest known pro- 
ton emitter ^^^Bi {Z = 83) decays with the I = proton 
emission from the 1/2+ [400] intruder state pushed by the 
deformation above the Z = 82 shell. 

It is worth noting that the s = 1 ^'^^Ho proton emit- 
ter is among the most extensively studied and under- 
stood isotopes beyond the proton-drip line. We know 
proton emission from both ground and isomeric states in 
^''^Ho to the 0+ ground state as well as first excited 2+ 



TABLE VI Proton energies Ep and branching ratios Jp^*" 
measured for proton emission channels from ^''^Tm together 
with the calculated values based on the particle-core vibration 
couphng model ( |Hagino[ |200T| [Tantawy et fflH|2006[ ). E/ de- 
note excitation energy of the final state in ^''^Er. The spin and 
parity for states in "^Er is 1/2+ , 3/2+, 11/2", and 13/2" for 
the excitation energy of 0, 175, 253, and 484 keV, respectively. 
All energies are in keV. 



Ep 


Ip^^ (%) Wave function composition 


real 
P 


Al 


Ef 






ground st&tG 












r ^5-, Ti/a = 68(5)ms 








938(4) 


13.8(9) 


2% 7rSi/2®!//!.il/2®0+ 







253 


1016(4) 


18.3(11) 


4% 7r/7/2 ® i/Si/2 ®2+ 


15 


3 


175 






41% 7r/iii/2 ® vsi/2 ® 2+ 


0.003 


5 


175 


1191(1) 68.1(19) 


53% 7r/iii/2 ® 1^31/2 ® 0+ 


70 


5 









isomeric state 












r = 10+, Ti/2 = 198(3)ms 








889(8) 


1.0(4) 


2.5% 7r/7/2 ® iyhii/2 ® 2+ 


1.2 


3 


484 






41% 7r/iii/2 ® iyhii/2 ® 2+ 


0.04 


5 


484 


1120(1) 


100(1) 


55% Tvhii/2 ® i^hii/2 ® 0+ 


98.6 


5 


253 






0.1% 7r/l9/2(8)!//lil/2®0+ 


0.2 


5 


253 






0.4% tt{1 > 5) ® i//iii/2 









Value based on the experimental intensity ratios not predicted 
by the particle-core vibration coupling model 



State in ^^^Dy ( 


Davids et al. 1998 


Karny et al. 


2008 


Rykaczewski et al. 


1999|). In these decays three differ- 



ent angular momenta (Z = 0, 2, and 3) are involved 
Proton emission was also used in a recoil decay tagging 
study of this isotope allowing observation and interpre- 
tation of rotational bands up to I'^ = 35/2" (Seweryniak 
Theoretical works focused on ^""^^Ho include 



et al. 



2001). 



Arumugam et al. (2009); Barmore et al. (2000); Davids 
and Esbensen' f 2000; '2004 );'Esbensen and Davids (2001"); 



Fiorin et al.^ (^2003); Kadmensky and Sonzogni, ( ,2000( ); 
and Kruppa and Nazarewicz ( 2004 1 



2. Even-mass, s — 2 proton emitters 

There are 17 even-mass s = 2 proton emitters known. 
The lightest known is 5PCS57 and gpTlgs is the heaviest. 
In these odd-odd nuclides the wave function compositions 
results from the interaction of the unpaired proton and 
neutron. For nuclei with the neutron number iV < 82 
the valence neutron has no significant infiuence on the 
proton emission, i.e. proton emission from the odd-odd 
emitters follows the pattern of the odd-even emitters of 
the same element. For example, ^^°Lu has two states 
decaying via proton emission: the ground state which 
emits I = 5 proton and an isomeric state emitting proton 
with I = 2. The same pattern is found in ^^^Lu, where 
the ground- and isomeric states decay with I — 5 and I = 
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TABLE VII Properties of s > 2 proton emitters. 



Emitter 


E*(MeV) 


Ep(MeV) 


Ti/2 


Configuration 


Ang.mom. 


References 


■"'"Co 

94m^g 


3.197(29) 

6.457(1) 
5.780(30) 


1.57(3) 

1.28(5) 
0.79(3|^ 


242(15) ms 

152(4) ns 
0.39(4) s 


[nfy\(»iyfyl]19/2- 

k(/r>3/2)®'^/772]10+ 


1=9 

1=5 
1=4 


Jackson et al. 


( 1970 1 


Cerny et al. ( 


19701 


Rudolph et a 


. (20081 


Mukha et al. 
Cerny et al. ( 


(20051 
20091 



^ A second transition of 1.01 MeV reported by [Mukha et a/!] | |2005[ l has not been confirmed by the work of [Cerny et al.\^2009\ and 
therefore is not included in the table. 



2 proton emission, respectively. The situation changes 
when neutrons start fiUing 1^/7/2 orbital above N = 82. 
In these cases an attractive interaction due to the tensor 
force between 77^3/2 and 1^/7/2 orbitals pushes the former 
above the 7rsx/2 orbital leading to ^ = 2 proton emission 
from the ground state. In other words, the Si/2 ground 
state in odd-mass emitters is replaced by the ^3/2 ground 
state in odd-odd emitters of the same element. The I = 2 
emission from the even-mass isotopes has been observed 
for Ta83, Regs, Irgg, ^^"Augi. The ^I^Tlgs emits 
axil = proton from its ground state just like its neighbor 

177T1. 

The ^^^Tm, s = 2, proton emitter is the richest proton 
emitter known. There are 5 proton transition known in 
this case (Tantawy et al. 20061. Three transitions are 



coming from the ground state and two were assigned to 
its isomeric state. Table |VI| shows experimental results 
obtained for these transitions together with the calcu- 
lated wave functions components, based on the work of 



Tantawy et al. (20061 and Hagino (20011 



The ground state emits protons with Z = 0, 3 and 5, 
while the emission from the isomeric state has mainly 
the I = 5 component. It is worth noting that in the cited 
calculation I = emission is due to the 7rsi/2(8iJ^ft.ii/2<8)0'^ 
2% component, which is the isospin symmetric to the 



dominant (53%) tt/i 



11/2 



Sl/2 



.0+. 



not present in the proposed configuration of the emitting 
state. Both cases can not be described by the core plus 
proton model used in case of s < 2 proton emitters. In 
the case of ^'*™Ag / = 4 emission is assumed to originate 
from TTgg/2 orbital. The 21+ isomeric state is created by 
three proton holes on the gg/2 orbital coupled to three 
neutron holes on the gg/2 orbital. Table VII shows the 



combined information on these high spin proton emitters. 



D. Outlook 

The wealth of nuclear data established by proton ra- 
dioactivity studies is impressive and indicates that this 
field of research is mature and the applied experimental 
techniques are well advanced. They appear, however, to 
be still not sufficient to address potential proton emitters 
with atomic numbers below 50. The low production rates 
and short half-lives, expected for these nuclei, present a 
challenge to the experimentalists. Observation of pro- 
ton radioactivity in nuclei with Z < 50 and establishing 
their properties will be important for calculations of the 
astrophysical rp-process. Of special interest are nuclei 
ar ound the waiting po ints, like ^*Se (see a recent article 
of Rogers et al. (2011 1), and the region just below ^"°Sn, 



at the expected end of the rp-process path. 

The prospects for experimental studies of proton emis- 
sion in the region between N = 82 and Z — 82 were 



C. Seniority s > 2 proton emitters 



discussed recently by Page (20111. 



In this category there are 3 proton emitters to be men 
tionedj_Jf^Co dJackson et d] |1970|, ^'^"'N i ( [Rudolph 
et al.\ |2008| and ^^'"Ag ([Mukha et a^|2005|. They are 



all high spin, high excitation isomers with a multi-particle 

configuration of the wave function. The ^^™Co was the 
first proton emitter discovered. T1/2 — 242(15) ms. 
Its wave function is best described by the [^fy\ <8i 



configuration. The transition goes to the 

0+ ground state of ^^F e, thus the proton ca rries 9 units 
of angular momentum (Jackson et al. 19701. The ^^"Ni 



is the first and so far the only proton emitter produced in 
fragmentation reaction. The angular momentum of the 
emitted proton equals 5, although the 7r/i]^]^y2 orbital is 
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VI. ALPHA DECAYS 

Emission of a particles belongs to the oldest known 
(together with /?-decay) types of radioactivity. Its first 



4000 



5000 



6000 



theoretical description by Gamow ( 1928 ) and indepen 



dently by Gurney and Condon ( 1928 1 was one of the ^5 



early triumphs of quantum mechanics applied for the first 
time to the atomic nucleus. In particular the empirical 



law of Geiger and Nutall (19121 could be successfully 
explained. Presently, the calculations of the a-decay life- 
times are performed in analogy to proton radioactivity 



by using Eq. 21 (Gurvitz and Kalbermann 19871 where 



the proton spectroscopic factor is replaced by the a 
preformation factor Sa. The latter measures the prob- 
ability that the a particle is formed inside the mother 
nucleus. Combining the shell model with the cluster 
model proved to be successful in calculating the abso- 
lute a decay width of 2i2po ( |Varga et al| |1992[ ). The 
result, r — 1.45 • 10~^^MeV, agrees very well with the 
experimental value of 1.5 • 10~^^MeV. The large body 
of experimental and theoretical findings about a decay 
mode is covered extensively in a number of books and 



reviews (Delion 2010 Rasmussen 1966 Roeckll 19961. 



For a compilation of even-even a-decay data see [Akovali 
(19981. An extended version of the Geiger-Nutall rule 
has been recently proposed by Qi et al. ( 2009a|b l. Since 
this paper is devoted to the decays at the limits of stabil- 
ity, here we focus mainly on the latest studies of a decay 
close to the proton drip- line. 

One of the regions which attracts attention is located 
above ^""Sn, where due to the proximity of = 50 and 
Z = 50 shell closures, the energy available for a decay is 
large enough to overcome the Coulomb barrier. This re- 
sults in an island of a radioactivity for 52 < Z < 56 and 
the neutron number A^ up to 60. Apart from the energy 
factor one has to note that the nuclei in this region are 
among the heaviest with protons and neutrons occupying 
the same type of single-particle orbitals. For these nuclei 
the active single-particle orbitals are 37/2 and ^5/2 which 
differ in excitation energy by only a few hundred keV. 
In the case of protons and neutrons occupying the same 
orbitals, their spatial overlap is maximized leading to the 
maximal preformation factor. For this reason alpha de- 



cays ^"*Te 



100 



Sn and ^°^Te 



102 



Sn are expected to 



be the best examples of the superallowed alpha decay 



( [Macfarlane and Siivola[|1965[|Roeckl[|1995 l. 
While the search for the superallowed ^°^Te 



100 



Sn 



decay is still an ongoing effort, the successful measure 
ments of ^°^Te — >■ ^°^Sn have been reported (Liddick 
et al. 2006| Seweryniak et al. 2006). In the work of 
Seweryniak et aZ. ~ 2006 l the decay of ^°^Te was mea- 
sured directly. Ions of ^'^^Te were produced in the fusion- 
evaporation reaction of a ^^Ni beam impinging on a ^'^Cr 
target. The products were separated by means of the 
Fragment Mass Analyzer of Argonne National Labora- 




50 100 150 200 250 300 

Energy (keV) 



350 



FIG. 24 (Color online) (a) Energy spectrum of the first 
(^"^Xe) and the second (^''^Te) a pulses obtained from the 
a — a pileup traces (inset). The lines at 3910(10) and 
4063(4) keV are assigned to the ^°^Xe-i>^°'^Te transitions, 
while the lines at 4711(3) and 4880(20) keV are assigned to 
the ^"^Xe— >^''^Sn decay, (b) 7 spectrum in coincidence with 



the analyzed a — a traces. From (Darby et al. 20101 



detector. Alpha decays events were measured and cor- 
related with the implanted ^°^Te ions. Thirteen counts 
were identified as representing the decay ^"^Te ^'^^Sn. 
As a result, an a decay energy — 4720(50) keV, cor- 
responding to = 4900(50), and a half-life of — 
0- 7^n ^7 us were establis hed. In the different experiment 

of ■ ■ 
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Liddick et al.\ (|2006| the ^^Ni beam impinging on a 
Fe target was used. With the Recoil Mass Separator 
(RMS) of the HRIBF facihty in Oak Ridge, the ^°^Xe 
— > ^"'^Te — > ^"^^Sn decay chain was analyzed. The rela- 
tively long half-life of -"^"^Xe (~ 13 ms) helped to over- 
come the inevitable losses due to the finite flight time 
through the RMS (~ 2^s) in case ^°^Te was studied di- 
rectly. Pulses of correlated a — a decays were analyzed 
with help of digital electronics programmed to trigger and 
collect only the signals of high energy decay events. Anal- 
ysis revealed two branches of -'^^^Xe a decay with energies 
Eai = 3918(9) keV and Ea,2 = 4062(7) keV, followed by 
the a decay of i°^Te with the energy E^ = 4703(5) keV. 
The half-lives of 13(2) ms and 0.62(7)/j,s were determined 
for ^''^Xe and ^'^^Te, respectively. The preformation fac- 
tors derived from these experiments were found to be 
larger by about a factor of three from the values in the 
well studied region of the doubly magic ^"^Pb (Liddick 



et al. 2006 Mohr 2007 1 



tory (Davids et al. 19921 and implanted into a DSSSD 



In the next experiment at the HRIBF laboratory, using 
a similar technique, the fine structure in the a decay of 
lo^Te was found (iDarby et aL||2010|. The setup used by 
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Liddick et al. { 2006 1 was additionally equipped with four 



germanium clover detectors placed around the DSSSD 
implantation detector. Double pulses of a-a events from 
the ^°^Xe — ^ ^"^Te ^"^Sn decay chain were stored and 
analyzed. The result is shown in Fig. |24] The observed 
a lines at 4711 keV and 4880 keV are assigned to the 
decays of ^"^Te leading to the first excited and to the 
ground state of ^'^^Sn, respectively. The 7 transition be- 
tween these two states (172 keV), coincident with dou- 
ble a pulses was also detected (Fig. p4](b)). This re- 
sult confirmed the previous evidence for the first excited 
state in ^°^Sn obtained by Seweryniak et al. (2007b I who 



employed the recoil-decay tagging method (RDT) (Paul 



half-life all three alpha transitions were assigned to the 
decay of the ^^°Po ground state. Analysis of coincidences 
between a particles, electrons, X- and 7-rays suggested 
that the spin of the three final states is 0"'". The analysis 
of the preformation factor lead to the conclusion that the 
presence of the three 0"'" states in ^®^Pb similar energy, 
is a manifestation of shape coexistence where the ground 
state is spherical, the first excited 0+ state at 532 keV is 
oblate and second excited 0+ state at 650 keV is prolate. 

It is worth noting that in recent years a decay served as 
a tagging signal in recoil decay tagging studies of heavy 
nuclei providing valuable nuclear structure information. 
Recent highlights from the RITU spectrometer at Uni- 



et al. 1995 1 by combining 7 spectroscopy with /3-delayed versify of Jyvaskyla were presented by Julin (2010 1. The 



proton detection. 

The stronger a line in the decay of ^°^Te, at 4711 keV 
with the intensity of 89(4)%, is interpreted as correspond- 
ing to the decay with no change of the orbital angu- 
lar momentum {Al — 0), while the line at 4880 keV is 
assigned to the Al = 2 channel (Darby et al. 
In addition, the strong Al = transition goes to the 
first excited state in ^"-'^Sn in contrast to the a decay 
of ^"^Te where the strong Al = decay connects nu- 
clear ground states. Thus, the level inversion occurs 
between ^°^Sn and ^"■^Sn. The 5/2+ ground state in 
^°'^Sn becomes the first excited state in ^"^Sn, while the 
7/2+ excited state in ^°'^Sn becomes the ground state 
of "'^°^Sn. This phenomenon is interpreted as a result of 
the interplay between the pairing on the 1^57/2 orbital 
(yP"''^'^ {§7 /2) = 1.4 MeV) being much stronger than for 
the iydr,/2 orbital (^'"'''■(^5/2) = 0.56 MeV), and a small 
(0.17 MeV) energy difference between these two orbitals 
( Darby et al. |2010[|. This interpretation contradicts the 



conclusions of Seweryniak et al. (2007b I who assigned 
spin and parity 5/2+ to the ground state of ^°^Sn. We 
note that both experiments agree on their common ex- 
perimental finding but differ in theoretical interpretation. 
Further experiments are required to firmly establish the 
^5/2 ~ 97/2 order in ^"^Sn. For example, an observation 
(or exclusion) of the Gamow- Teller /3-decay between the 
ground state of ^"^^Sn and the 9/2+ ground state of ^"^In 
should settle the controversy. 

In the region of very neutron-deficient lead isotopes, re- 
cent a decay studies provided information on the shape 
coexistence in ^®^Pb (Andreyev et al. 20001. The states 



of °°Pb were populated in the a decay of ^^Po, pro- 
duced in the fusion-evaporation reaction of ^^Cr beam 
impinging on a 



^^^NdFa target. The products were sep- 



arated by means of the SHIP velocity filter (Miinzen- 



berg et al. 1979 ) at GSI Darmstadt and implanted into 



a position sensitive silicon detector, backed by a germa- 
nium clover detector for X-ray measurements. A set of 
silicon detectors was mounted for measurement of con- 
version electrons. In addition to the ground-state-to- 
ground-state decay, two other channels were observed in 
coincidence with conversion electrons. Due to the similar 



RDT experiments with the FMA separator coupled to 



germanium detector arrays were reported by Carpenter 
et al.\ (|1997|), Reiter et al. (1999), and by Seweryniak 



et al. (1998 



2005b I . Finally, it should be mentioned that 



most of the discoveries of new elements rely on a decay 



2OIOI. (Hofmann 2009a) 



40 



VII. TWO-PROTON RADIOACTIVITY 
A. Introduction 

The two-proton radioactivity (2p) is the most recently 
observed type of decay and thus the least known. The 
experimental studies are still in the early stage. The 
detailed understanding of its mechanism requires novel 
theoretical approaches which in certain aspects are still 
under development. A very early look can be found in 
the book by Baz' et al. (19721. The current experimental 



and theoretical status of the 2p decay were summarized 
recently by Blank and Ploszajczak f2008) and by Grig- 



orenko (20091 with focus on specific theoretical methods. 
Because of its exceptional status, we discuss here this 
decay mode in more detail, emphasizing the major qual- 
itative features of the phenomenon. The illustrations are 
provided mainly by the examples of ^Be, ^^Mg, and ^^Fe. 
These nuclei belong to p, s-d, and p-f shells respectively 
and their lifetimes span about 18 orders of the magni- 
tude, providing support for universality of the currently 
achieved understanding of the two-proton decay. 

The emission of two protons from a nuclear state is 
in principle possible in various decay scheme situations 
which are sketched in Fig. [25] We introduce here the fol- 
lowing notation: Et is the system energy relative to the 
nearest three-body breakup threshold, while £'2,- is the 
lowest two-body resonance energy relative to this thresh- 
old. The 2p decay in the pure form, which we will call the 
true 2p decay (or true three-body decay) is represented 
in Fig. 25 ^c) . In this case sequential emission of protons 
is energetically prohibited and all final-state fragments 
are emitted simultaneously. Such a situation is common 
among even-Z nuclei at the proton-drip-line and results 
from pairing interactions, see Sec. [IT] The decay dynamics 
of true 2p decay is not reducible to the conventional two- 
body dynamics and should be addressed by the methods 
of few-body physics. 

A somewhat special situation, represented 
Fig. 25 d), occurs when the ground state of 



m 
the 

subsystem is so broad that the emission of the first pro- 
ton becomes energetically possible (although E2r > Et) 
which opens a way for a sequential transition. Similarly, 
the decay may formally proceed in a sequential manner 
(i?2r < Et) but the ground state of the subsystem is 
so broad that no strong correlation between outgoing 
fragments at given resonance energy can be formed, see 
Fig. 25 e). We refer to such scenarios as democratic 



decays and discuss them briefly in Sec. |VII.B.1[ 

The three-body character of the 2p radioactivity places 
it in the broader context of nuclear processes exhibiting 
essential many-body features. This includes studies of 
the broad states in continuum and excitation modes, like 
the soft dipole mode (Aumann 2005 1 . Another topic. 



romean property that none of the three two-body subsys- 
tems is bound. The 2p decay can be seen as an analogue 
of the two-neutron halo, requiring similar ingredients in 
the proper many-body description of its properties. The 
illustration of this point is provided by the isobaric mir- 
ror partners ^He and ^Be: the first is the "classical" Bor- 
romean halo nucleus and the second is the lightest true 2p 
emitter. The crucial difference, however, comes from the 
fact that the 2p decays involve charged particles in the 
continuum which significantly complicates the theoreti- 
cal description. Another example: ^^Ne is a Borromean 
two-proton halo nucleus, while the first excited state of 
^^Ne and the less bound ^^Ne are true 2p emitters. 

All ground-state two-proton emitters studied experi- 
mentally up to now are collected in Table |VIII[ 



1. Two-proton correlations 

The two-body decay of a resonance is characterized 
only by the energy and the width of the state. The three- 
body decay is much more "rich" as complex information 
about momentum correlations becomes available. 

For decays with three particles in the final state there 
are nine degrees of freedom (spins are not counted). 
Three of them describe the center-of-mass motion and 
three describe the Euler rotation of the decay plane. 
Therefore, for a fixed decay energy Et there are two pa- 
rameters representing the complete correlation picture. 
It is convenient to choose the energy distribution param- 
eter £ and an angle Ok between the Jacobi momenta k^, 
and k,,: 



Et — -El- 



pursued actively in the last decades, is the phenomenon 



of two-neutron halo (Zhukov et al. 19931 with its Bor- 



^ ^ iAi + A2)kl {Ai+A2 + A3)kl 
2MA1A2 2M{Ai+A2)A3 ' 

_ A2ki - Aik2 _ A3(ki+k2)-(Al+A2)k3 

"~ A,+A2 ' ^1+^2+^3 

e^EJEt , cos(0fc) = (k, • k„)/(/c, fcj,) , (26) 

where M is "scaling" nucleon mass, Mnucieus = M(Ai + 
A2 + A3). The Jacobi momenta for two-proton emitters 
(protons are indistinguishable) can be defined in two "ir- 
reducible" Jacobi systems, called "T" and "Y", see Fig. 
[26) In the "T" Jacobi system, the core is the particle 
A3 and the parameter e describes the energy distribu- 
tion between the two protons. In the "Y" Jacobi system, 
the core is the particle A2 and e corresponds to the core 
and proton subsystem. The Jacobi momentum k^, is the 
momentum of particle 1 in the cm. of particles 1 and 2, 
hy is the cm. momentum of particles 1 and 2 in the cm. 
of the whole system (particles 1, 2, and 3). Distribu- 
tions constructed in the different Jacobi systems are just 
different representations of the same physical picture. 

A more general (5-dimensional) correlation pattern be- 
comes available for systems with total spin J > 1/2. 
Manifestation of such correlations requires existence of 
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(A-l)+p {A-2) + 2p 



(A-\)+p (A-l) + 2p 
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{A-\)+p {A-2) + 2p 



FIG. 25 Energy conditions for different modes of the two-proton emission: (a) typical situation for decays of excited states 
(both Ip and 2p decays are possible), (b) sequential decay via narrow intermediate resonance, (c) true 2p decay. The cases 
(d) and (e) represent "democratic" decays. The gray dotted arrows in (c) and (d) indicate the "decay path" through the states 
available only as virtual excitations. 



a selected direction in space and spin alignment, which is 
naturally achieved for short-lived states populated in nu- 
clear reactions. The only example of such detailed stud- 
ies is so far the two-neutron decay of broad states in 
2005| ) 



(Golovkov et al. 



2. Historical note 



The possibility of a true two-proton emission was men- 



tioned for the first time by Zeldovich (19601. This work 
comprises the dripline prediction for light systems. Af- 
ter predictions about existence of ^'^O and ^°Mg isotopes 
Zeldovich notes: "The existence of ^^O, ^^Ne, and ^^Mg 
is not excluded [...] These nuclei could appear to be 
unstable with respect to emission of two protons simul- 
taneously." The explicit and detailed statement of the 
two-proton radioactivity phenomenon was given by V.I. 
Goldansky a bit later ( Goldansky 1960[ p" While the pro- 
ton and cluster radioactivity are quite straightforward 
generalizations of a-radioactivity, the few-body decays 
are qualitatively different and required ingenuity to fore- 
see. The pioneering work of Goldansky contained several 
important insights which we illustrate by the following 



citation (Goldansky 19601 



"Thus the simplest approach to the theory of two- 
proton decay would consist in using the product of two 
usual barrier factors, that is, in introducing an exponen- 



^ Zeldovich and Goldansky lived next door to each other. The 
problem is known to be a subject of many of their informal dis- 
cussions (which are acknowledged in the paper of Goldansky). 
Later, on occasion of the priority discussion raised by some peo- 
ple, Zeldovich rejected any credits for the idea. Zeldovich was 
famous for providing in his works insights important for later 
development of physics in a very compact form and without at- 
tempt of further elaboration. 



tial factor of the type 



w{e) = exp 




1 

7i 



(27) 

where is the sum of the energies of the two protons 
(energy of emitted diproton), e and (1 — e) are the frac- 
tions of energy referring to each of the protons. 

It can easily be seen that the total barrier factor 'w{x) 
is maximum for e = 0.5, i.e., when the proton energies are 
equal. It will be noted that the value in the exponent is 
just the same as for the sub-barrier emission of a diproton 
with the energy E'j' BS ci whole.'Q 

The general character of the energy distribution pre- 



dicted by Eq. (27 1 has proven to be correct and is now 
confirmed also experimentally. The idea of emission of 
a "diproton particle" turned to be an attractive concept 
but finally appeared to be misleading. 

Later, significant theoretical work was devoted to iden- 
tifying the best candidates for the observation of the 2p 
radioactivity. Due to the extreme sensitivity of the 2p 
decay probability to the width of the Coulomb barrier, 
the decay energy of a candidate must fall into a rather 
narrow window ( Nazarewicz et al. 1996| . The resulting 



2p partial half-life should be long enough for an efficient 
separation in the spectrometer (typically a fraction of 
a microsecond) but short enough to compete with the 
/3+ decay channel (~ 10 ms). Thus accurate mass pre- 
dictions for nuclei beyond the drip line were necessary. 
One of the most exact methods was the application of 



the isobaric multiplet mass equation (IMME) (Benenson 
and Kashy 1979 ) combined with the experimentally mea- 



sured mass of the neutron-rich member of the multiplet. 



^ In Eq. |27| we have modified the notation of Goldansky to make 
it consistent with the notation of this work. 
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TABLE VIII Ground-state 2p emitters investigated experi- 
mentally. The indicated half-life corresponds to the partial 
value for the 2p decay. 



«Be 

12q 



"Ne 



19 



Mg 

«Fe 



E keV 
1371(5) 

1820(120) 
1790(40) 

1800(400) 
1350(80) 
1400(20) 
1350(80) 
750(50) 

1100(100) 
1140(50) 
1154(16) 



r or Ti/a 
92(6) keV 

400(250 

580(200 

600(500 

200(100 

110(4o]^l 

<200 keV 

4.0(15) ps 
3_2+2.6 

5.71 



Reference 



( Whaling 19661 



(KeKelis et al. 19781 



( Kryger et al. 1995 1 



( Suzuki et al. 2009| 



(KeKelis et al. 19781 



(Woodward et al. 19831 



(Mukha et al. |2008bl 



( Mukha et al. 20071 



ms (Pfiitzner et al. 20021 



( Giovinazzo et al. 2002 1 



( Dossat et al. 20051 



■^^Ni 1350(20) 
^■*Zn 1480(20) 



'0 



jMiernik et al. 2007c I 



(iD ossat et al. 
(^Pomorski et al 



(Blank et al. 2005 \ 



20051 



2011b I 



^ According to theoretical calculations much smaller widths are 

expected | |Barker| [l999l |2001[ [Grigorenko et a;.[|2002| |. 
^ Only one decay event observed! 

Coefficients of the IMME can be calculated within the 
shell model or deduced from the Coulomb-energy system- 
atics. Both approaches were undertaken (Brown 1991 



Cole 1996 Ormand 1996) and the choice of the best 



candidates was narrowed down to three cases: "^^Fe, ''^Ni, 
and ^"'Zn. These predictions played an essential role in 
motivating experimental search for the 2p radioactivity. 



B. Experimental results 

While the discovery of the true 2p radioactivity awaited 
more than 40 years from its prediction, interesting infor- 
mation about related phenomena has been accumulated. 
In particular, this include 2p democratic decays and the 
2p emission from excited nuclear states populated in 13 
decay and in reactions. In the following we overview the 
main steps in the experimental progress. 



1. Democratic decays 

Since it is relatively easier to reach experimentally light 
proton-drip line nuclei, the first attempts to search for 2p 
emission phenomena started with the light 2p unbound 
systems: ^Be (Bochkarev et al. 



1984 



1978 



Geesaman et al. 



Kryger et al. 19951, 



1977 ), (KeKelis et al. 

and '^'^Ne ( jKeKelis et aZ][l978[ ), see also Table |VIII[ Due 
to a very low Coulomb barrier these cases have half-lives 
much too short to classify as being radioactive, more ap- 
propriately they should be referred to as 2p resonances. 



(a) Jacobi "T" (b) Jacobi "Y" (c) Shell model "V" 
O 



core 




FIG. 26 (Color online) Coordinate systems for two-nucleon 
plus core problem. In the Jacobi "T" system (a), the "dipro- 
ton" and the core are explicitly in configurations with definite 
angular momenta l^ and ly. For a heavy core the Jacobi "Y" 
system (b) is close to the single-particle "V"-system (c) used 
typically in many-body approaches. 



Already Geesaman et al. (1977) could not fit the re- 
sults on ^Be using simple decay scenarios (phase volume, 
diproton decay, simultaneous emission of p- wave protons) 
and concluded that a full three-body computation is nec- 
essary to understand the measured spectra. Later the in- 
terpretation of correlations observed in the decay of the 
^Be ground state have led to the concept of "democratic 
decay" (Bochkarev et al. 1989). In such a decay there 



is no strong energy focusing of the particles and they 
are smoothly distributed in the momentum space. The 
system is demonstrating a kind of equal rights among dif- 
ferent parts of the kinematical spac^ The democratic 
decay is not a phenomenon on its own but rather a name 
for the experimental fingerprint of a true three-body de- 
cay in light systems with the relatively small Coulomb 
force or in two-neutron emitters. 

The study of opening angle between protons emitted 
from ^^O was motivated by the search for diproton cor- 
relation. The measured spectrum, however, was found to 
be consistent with the sequential emission via intermedi- 
ate "'^"'^N state. Later it was found that indeed, the ground 
state energy of ^-"^N is below that of ^^O and the decay of 



the latter belongs to the class shown in Fig.[25[e) (Azhari 
[erZ| [l998| ). 

Recently, the full correlation picture for protons emit- 
ted by ^Be has been experimentally established (Grig- 



orenko et al. 2009b I and was found to be in very good 



agreement with the predictions of the three-body model 
which will be discussed below. Both experimental and 
theoretical distributions in the "T" and "Y" Jacobi co- 
ordinates are presented in Fig. |27j Similarly, the new 
results of p-p correlations in the decay of ^^Ne could be 



well described by the three-body model (Mukha et al. 



^ Now the term democratic decay has become accepted for descrip- 
tion of the mentioned class of phenomena. There exist, however, 
an anecdotic story that when it was used for the first time at 
some conference at Soviet Union in the end of 80-s the authors 
were heavily criticized because "there can not be any democracy 
in nuclear physics". 
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FIG. 27 (Colour online) Complete correlation picture for ®Be 
g.s. decay, presented in "T" (a,c) and "Y" (b,d) Jacobi systems. 
The (a,b) is theory, (c,d) is experimental data. Qualitative 
illustration of the meaning of different kinematical regions is 
provided above the panels. Data and calculations are from 



by 7 radiation was found to dominate. In a different 
approach, however, higher excited states in ^^Ne were 
populated by In stripping reaction fr om the radioactive 
^*Ne beam (Zerguerras et al. 2004). The 2p angular 



correlation was found to be peaked at small angles in- 
dicating a contribution from a seemingly non-sequential 
mechanism. Due to small statistics and limited infor- 
mation on the identity of the decaying states, however, 
no definite conclusions could be reached. The 1~ reso- 
nance at 6.15 MeV in ^*Ne was also thought to provide 
opportunity for the true 2p emission, as no states are 
known in ^^F to be located within the decay energy win- 
dow. The 2p emission from this state, populated in the 
^^F+^H reaction, was measured and a diproton-like con- 
tribution to the p-p correlation spectra was claimed ( del 
Campo et al. 20011. Similar evidence was announced in 



another work, in which the same state was populated by 
the Coulomb excitation of the ^^Ne beam (Raciti et al. 



2008 1. The statistical significance of both claims, how- 
ever, is weak. The excited states of ^^Ne, populated in 
the inelastic scattering of ^^Ne on hydrogen target, were 



found to emit two protons sequentially (F. de Oliveira 



Santos et al. 2005 1 



An interesting case is the claimed 2p emission from the 



high-spin 21+ isomeric state in ^''Ag ( Mukha et al. 2006 1 



The authors explain an anomalously high 2p decay rate 
and observed p-p correlations by assuming a very high 
deformation of the isomer. This work has caused a lot of 



controversies 


(|Cerny et al. 


2009 


Kankainen et al. 


2008 


Mukha et al. 


|2008a 


Pechenaya et al. 20071 which call 



( Grigorenko et al. 2009b I 



for further experimental investigations. 



2008b). 



3. Two-proton radioactivity 



2. Two-proton emission from excited states 

Several cases of /3-delayed 2p emission are known, as 
mentioned in Sec. |IV.D| They are discussed in more detail 
in Blank and Borge, ( ,2008 1 . It is believed that in all cases 



studied, the mechanism of the emission is sequential, i.e. 
it can be described as a sequence of two two-body decays. 

In a few other cases, the 2p emitting excited states 
were populated by nuclear reactions. Like in the case 
of ground-state resonances, the main motivation of these 
studies was the search for p-p correlations going beyond 
the sequential mechanism. The 2+ state at 7.77 MeV in 
^"^O was excited by the two-proton transfer reaction from 
■^He impinging on a '^^C target (Bain et al. 1996). The 



The first case of 2p radioactivity was found in 2002 in 
the decay of '*^Fe measured at GSI Darmstadt (jPfiitzner] 



et al.\\2002\ and at GANIL ( [Giovinazzo et aZ.[|2002[ ). In 
both experiments ions of ^^Fe were produced by the frag- 
mentation reaction of a ^®Ni beam and separated using 
the in-flight technique (Sec. III.B.ll. The selected ions 



were implanted into silicon detectors and the only ob- 
servables measured were the decay time and energy. It 
sufficed to claim the observation of a new decay mode be- 
cause only the 2p emission hypothesis was consistent with 
the measured data. Later, in another GANIL experiment 
the information on ""^Fe was obtained with larger statis- 



tics (30 atoms) and improved accuracy (Dossat et al. 
2005[ ), see Table |VIII| and Fig. [28] With the same exper 



imental technique, applied again at GANIL laboratory, 
^■*Zn has been identified to decay by the 2p radioactivity 



(Blank et al. 2005). One decay event of ''^Ni was found 



2p emission from this state was found to proceed sequen- 
tially through a 1/2+ state in ^'^N. The first excited state 
in ^"^Ne (3/2~) was reached by the Coulomb excitation to coinc ide with the 2p dec ay energy predicted for this 
of the radioactive ^^Ne beam to search for its 2p decay 
(Chromik et al.\ 2002). Although the conditions for the 



nucleus (Dossat et al. 2005) 



true 2p decay were fulfilled in this case, the de-excitation 



In the next step, gaseous detectors, based on the prin- 
ciple of the time projection chamber (TPC), were devel- 
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FIG. 28 (Color online) Lifetime of '*^Fe vs. decay energy cal- 
culated in different models, (a) Simplified models of 2p decay. 
All spectroscopic factors are taken as unity, (b) Three-body 



model (Grigorenko et al. |2003a Grigorenko and Zhukov 



2007a I and continuum shell model results (Rotureau et al. 
2006p . The experimental points demonstrate the rapid im- 



provement of the data: circle ( Pfiitzner et al. 2002 1 , square 
flGiovmazzo ef 'fflTI |2002| ) trian gle ( |Dossat ef a/.[ |2C)05[ ), and 
diamond (IMiernik et al. |2007c|). 



oped to directly record emitted protons and to establish 
the correlations between them. The first direct observa- 
tion of the two protons ejected from ''^Fe was achieved 



by Giovinazzo et al. (20071 who recorded projections of 
protons' tracks on the anode plane of the TPC. Later, 
this detector was used to directly de monstrate the two - 
protons emitted in the decay of ^'*Zn (Blank et al 



20111 



Miernik et al. (2007c I applied a novel type of detector, 
utilizing the opt ical readout of the TPC signals (OTPC, 
see Sec. III.C.ll, to the detailed decay study of ^^Fe at 



the NSCL/MSU laboratory and succeeded to fully re- 
construct tracks of emitted protons in three dimensions. 
The full correlation picture for the 2p decay of ^^Fe es- 
tablished in this experiment is shown in Fig. |29[ Re- 
cently, the OTPC detector was used to the decay study 



of ''^Ni at the NSCL/MSU laboratory and provided a di- 
rect and unambiguous evidence for the 2p radioactivity 



of this nuclide (Pomorski et al. 2011b I. From the six de 



cays recorded, four corresponded to the 2p emission and 
two were interpreted as /3-delayed proton emission. An 
example of a 2p decay event of ''^Ni is shown in Fig, 
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An application of the different technique based on de- 



cay in-flight and particle tracking (Sec. III.C.ll is nec 



essary to extend decay studies to very short half-lives in 
the subnanosecond range. This technique was used by 



Mukha et al. (20071 to investigate the 2p radioactivity 
of ^'^Mg whose half-life was found to be 4.0(15) ps. The 
p-p momentum distributions, projected on the transverse 
detector plane were obtained in ( Mukha et al.\ 2008b I, 
see Fig. pi] 



C. Simplified theoretical models 

Simplified theoretical models of three-body decays are 
typically based on the factorization of the decay ampli- 
tude into the product of two-body terms. For 2p emission 
this factorization can be done either in the "T" or in the 
"Y" Jacobi systems, see Fig. [26) This results in formula- 
tion of diproton model or the direct decay model respec- 
tively. Factorization of amplitudes becomes possible only 



for specific forms of the underlying Haniiltonians (Grig- 
orenko and Zhukov 2007a I. Then, neglecting some of 



the final state interactions (FSI), the three-body Green's 
function (GF) can be constructed in an analytical form: 



G(+)(XY,X'Y') = ^ &G^+^(X,X') 



E 



(Y,Y'). (28) 



The operator G^j^'' is the ordinary two-body GF for the 
respective subsystem. This method allows to take into 
account one (out of three present) FSI exactly in dipro- 
ton model. In the direct decay model two FSIs can be 
treated exactly but only in the limit of infinitely heavy 
core, see Figs. 26 (b,c). For heavy 2p emitters this is a 



good approximation. The availability of approximate but 
analytic GF allows to determine decay properties without 
further approximations. 



1. Direct decay model 

The decay amplitude can be factorized in the "V" co- 
ordinate system [Fig. [26|c)| if we neglect p-p interaction 
and also assume an infinitely heavy core. The dominating 
contribution to the width can be obtained in the single- 
pole approximation for two-body GFs where only the low- 
est states in the subsystems core+p with resonance ener- 
gies Ep-^ and Ep^ are considered (Galitsky and Cheltsov 



45 




FIG. 29 (Color online) Momentum density distribution on the kinematical plane {e,cos{dk)} for *^Fe. Diproton model (a,b), 
direct decay model (c,d), three-body model(e,f), and experimental distributio n (g,h). Correlation patt erns are provided in the 



"T" (a,c,e,g) and "Y" (c,d,f,h) Jacobi systems. The calculations (e,f) are from ( Grigorenko et al. 2010 1 and data (g,h) are from 



(Miernik et al. 2007c I 



1964): 



Et {Va\ 
27r 



^P^{eET) 



de 

isET-Ep,)^ + rp,{eET)yi 

(29) 



((1 - s)Et - Ep^Y + Tp, ((1 - e)ETY/A 



It was shown in ( Grigorenko and Zhukov) 2007b) that 
the matrix element (V3) can be well approximated as 

{V:if ^D^iET-Ep^^Ep^f ■ (30) 

The parameter is a constant. D3 ss 1.0 — 1.5. The Tp. 
is the width of the ground state resonance in the core+p 
subsystem. It can be expressed using the R-matrix ap- 
proach to nuclear reactions ( [Lane and Thomas 1958 ) : 



TpXE)^2-i^PiXE,R,Z^Z2 



(31) 



The penetrability Pi is defined via the Coulomb functions 



irregular at the origin: 



F and G ([ Abramowitz and Stegun[ |1968|, regular and 

(32) 



kR 



In the above rj = ZiZ2aM/k is the Sommerfeld parame- 
ter, a is the fine structure constant. The 7^ is the reduced 



width which is factorized into the "Wigner limit" j'^ri^: 

1 



2Afi?2 



(33) 



and the dimensionless reduced width 6^ (spectroscopic 
factor). The parameter R is the "channel radius", typi- 
cally taken as: 



i? = ro(Acore + l)'/', 



(34) 



where tq = 1.4 fm. 

The results of calculations in the direct decay approx- 
imation, Eq. (29), for the case of "'^Fe and for the values 



of angular momentum I — 1 and / — 3 are presented in 
Fig. |28[ These are likely to provide upper and lower 
width limits for given decay energy which are consis- 
tent with the limits obtained in the complete three-body 
model, see Sec. |VII.D.1| 

Eq. ( 29 1 has a similar structure and can be consid- 



ered as a refined version of the original expression (27) 



proposed by Goldansky. Compared to the latter it pro- 
vides quantitative results (the preexponent is defined) 
and brings the explicit dependence on the property of 
the core+^» subsystem. The energy distribution between 
the core and one of the protons predicted by the direct 
decay approximation has a narrow bell-shape peaked at 
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FIG. 30 An example of a two-proton decay event of ^**Ni 
recorded with the optical time projection chamber described 
in section [lII.C.l| Top: the image recorded by the CCD cam- 
era showing a long track of the **Ni ion entering the chamber 
from below and the two bright, short tracks of protons emit- 
ted 1.576 ms after the implantation. Bottom: a part of the 
time profile of the total light intensity measured by the PMT 
showing in detail the 2p emission. From (Pomorski et al. 
|2011b|). 



e — 1/2. which reflects the symmetry between emitted 
protons, see Fig. |32[ With increase of the atomic num- 
ber this distribution becomes narrower. The bell-shape 
profile of this inclusive distributions is by now well con- 



firmed by experiment (Blank et al. 2011 Grigorenko 
et al.\ |2009b| |Miernik et al.\ |2007c[ ) However, looking 



at the complete distribution (see Fig. 29 that shows the 
distribution for the case of angular coupling [p^y2]o) the 
deficiency of the model becomes evident. Qualitatively, 
the direct decay model angular distribution always has 
backward-forward symmetry, while there is a strong an- 
gular asymmetry observed in experiment. 



Using the fact that the numerator in Eq. ( 29 1 is sharply 
peaked at e/2 and the denominator for the true 2p decays 
is a smooth function within the decay window, it can be 
further approximated: 



rdir(£^T 



EtD^{Et — Epi — E, 



p2) 



2n{ET/2-Ep,)^ET/2-Ep2) 



X [ deTp,{eET)rp,{{l-e)ET) , (35) 
Jo 

This approximation allows to estimate how the decay 
width depends on the position of the lowest resonances 
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FIG. 31 (Colour online) Complete correlation picture for 
^^Mg g.s. decay, presented in "T" (a) and "Y" (b) Jacobi sys- 
tems. Comparison of two inclusive distributions projected 
on a plane is given in the lower row ( Mukha et al. 2008b I . 



The angles dp_p and 0^ are between projected momenta of 
protons (c) and projected Jacobi momenta in the "T" system 
(d). Solid, dashed, and doted curves correspond to three- 
body model, diproton model, and phase space simulations. 
Data are from (jMukha et al. 2008b I and calculations are 
from ( [Grigorenko et ttZ.||2010p . 



in the core-p subsystem. For the ground state true 2p 



emitters E, 



■pi 



Epi = E2r and we get : 

Tdir ^ [Et/2 — E2r) ^ 



(36) 



We note that a somewhat different expression for the 
width, in the analogous approximation, was introduced 
in Refs. (lAzhari et aZ.I|l998l|Kryger et aZ.I|l995| for ^^Q: 



rkr = 



2^ 



de 



Tp,{£ET)TpM-e)ET) 
{eET-Ep,Y+Tp,{eETY/A 



(37) 



The expression is given here in our notation. It was intro- 
duced by analogy with two-body R-matrix expressions. 



Using the same approximation as in Eq. ( 35 1 we obtain 



Tdir 
Ekr 



' {ET/2-Ep,nET/2-Ep,Y ■ ^ ' 



For the ground state true 2p emitters C « 4. In the series 



of papers (Barker 1999 2001 2002 2003 Brown and 



Barker 2003 Brown et al. 2002 ) formulae equivalent to 



Eq. (37 1 were used. The results obtained in these works 
should include the factor C varying between 2 and 4 to 
be consistent with the direct decay approximation. 
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FIG. 32 (Color online) Systematics of energy distributions in 
the "Y" Jacobi system (core-p channel) calculated in the direct 
decay model of Eq. (291. All distributions are normalized to 



unity at the maximum. 



(a) 






"0 / \ 











0.0 0.2 0.4 0.6 0.8 1.0 




FIG. 33 (Color online) (a) Classical one-proton flight path Rf 
for ^^O and ®Be connected with the time delay in the core+p 



subsystem calculated by Eq. (39 1, (b) Probability of d ifferent 
Rf values for the decay of ^^O estimated by Eq. (291. 



2. Simultaneous vs. sequential decay 

The true two-proton emission process is sometimes in- 
terpreted as "sequential decay via tails of higher-lying 
resonances". Indeed. Eq. (36 1 shows that such resonances 



are very important for the true two-proton emission pro- 
cess. However, we shall argue that the sequential inter- 
pretation may be misleading. It can be shown that the 
resonance scattering can be interpreted in terms of the 



time delay Ti{E,R) [e.g. (Baz 
TiiE,R) 



19671]: 



T{E)/4 



(39) 



Let us estimate the distance which the "first" emitted 
proton can travel while the "second" is "confined" in the 
tail of the resonance. Among the known and prospective 
two-proton emitters only the lightest ^Be and ^^O have 
such fiight path comparable or exceeding ~ 1 fm. It 
can be seen in Fig. 33 (a) that for ^^O the estimated 
flight path can achieve 8 fm, but the contribution of such 
situations to the total decay probability is minor, Fig. |33| 
(b). Furthermore, having a typical scattering length in 
the nucleon-nucleon channel around 20 fm even for ^^O 
we can not maintain that a reliable spatial separation of 
the core+p subsystem is present. 

So, the poles of the Greens' functions for two-body 
subsystems strongly increase the probability of true two- 
proton decay when they come close to the three-body 
decay window. However, only for one nuclear system 
(^^O) we can see considerable sequential separation of 
subsystems in space and even in that case the estimated 
contribution of sequential decay mechanism to the width 
is minor. 

Important insight about the transition from true three- 
body to sequential decay can be obtained using the direct 
decay model. In Fig. [34] the results of calculations for 
'^Be, 1^0, and ^^Mg with help of Eq. ^ are shown for a 
fixed energy i?2r of the resonance in the core+p subsys- 
tem as a function of the three-body energy Et- We see 



a difference between the democratic decays (broad states 
in the subsystems, overlapping with the decay window) 
and a radioactive decay (narrow state in the subsystem) . 
The curve for ^^Mg in Fig. 34 ^a) shows two components: 
three-body regime (true three-body decay) and two-body 
regime (sequential decay) with a narrow transition zone 
in between (a kink in the line). In the cases of ^Be and 
^■^O (democratic decays) there is no sharp difference vis- 
ible between the two decay regimes and the transition is 
smooth. 

The transition from the three-body to the two-body 
regime can be also illustrated by the energy distribution 
in the core-p channel. Again, there is a qualitative dif- 
ference between democratic and radioactive decays. In 
the former case, when the two-body resonance "enters" 
the three-body energy window, the width of the distri- 
bution suddenly becomes smaller, see Fig. 34 'b). On the 
other hand, if the state in the subsystem is narrow. Fig. 
]34[ c) , the two "horns" appear in the energy distribution 
corresponding to the proton energies defined by the in- 
termediate state. As long as the difference Ex — i?2r is 
small, these horns do not contribute significantly to the 
total width and the decay remains effectively of the three- 
body type. Gradually the horns become dominating, 
while the relative contribution of the three-body "bell" at 
e — Et/2 vanishes. It is reasonable to put the borderline 
between true three-body decay and sequential decay as 
the moment when the probabilities in the bell-like profile 
at e = Et/2 and in the sequential horns become equal. 
From this one can derive a criterion for the three-body 
decay mechanism: 



EqEt < E2 



(40) 



where for radioactive decays Eq ~ 0.76 — 0.84 and for 
democratic decays £0 ~ 0.3 — 0.5. Some uncertainty of 
the estimated Sq values here is connected with particular 
selection of parameter set for the system: Et, A, Z, and 
I. 

We note that the condition for the three-body decay 
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FIG. 34 (Color online) Transition from the true three-body to the sequential decay for ^Be, ^^O, and ^®Mg in the direct decay 
approximation Eq. (291. The width vs. decay energy Et (a), evolution of energy distributions in the "Y" Jacobi system (core-p 



channel) for the case of democratic decay (b), and for radioactive decay (c). The curves in (b,c) are normalized to unity for 
the maximum of the "bell" profile at £ ~ 0.5. 



proposed by Goldansky ( 1960 ) was more restrictive 



was taken i?o ~ 50— 100 keV (Goldansky 1961 Janecke 



Et + r2r/2 < i?2r 



(41) 



19651. The diproton model was often used in a trivial 



See also (Alvarez-Rodriguez et al. 2008a I for another Ormand 19971. The results with E^p = for different 



form with E'o = (Brown 1991 Nazarewicz et al. 1996 



formulation of how sequential and direct three-body de- 
cays may be distinguished. 



3. Diproton model 



Historically the first idea of the diproton approxima- 



ro values for^Fe are given in Fig. 28 (a). These results 



Later, it was investigated in more detail (Goldansky 



strongly overestimate 2p decay width and could be con- 
sidered as its strict upper limit. Eq. ( [45| corresponds to 
the "Coulomb-corrected phase volume". A radius of the 
p-p channel of Rpp ^ 1 — 2 fm should be chosen here to 
reproduce reasonably the low-energy p-p scattering pro- 
porties. A treatment of the density in the spirit of the 



tion appeared in the original work of |Goldansky| (|1960|) Migdal- Watson approximation, Eq. (|46|, was proposed 



1961 ) and gradually it has become one of the major ideas 



about the mechanism of the true 2p decay, frequently 
used as a synonym of it. In this approach a single "dipro- 
ton" particle (two protons in ^ = 0, 5 = state) is 
tunneling through the barrier along the straight classi- 
cal trajectory. The width in the "T" system is written 
as: 



by Barker (2001j). Here S(j{E) is the phase shift for the s- 
wave p-p scattering. Results of the approaches (45 1 and 



( 46 1 can be well approximated by expression Eq. ( 44 1 



with an appropriate choice of Eq . 



The diproton model in the form (42 1 has been inferred 



it has never been derived. To check the validity of this ap- 
proach the "dynamic" diproton model was developed by 



Jo 



Grigorenko and Zhukov ( 2007a I using the analytical GF 



Eq. ( 28 I of the simplified three-body Hamiltonian with a 



(42) 



diproton Ansatz. It was demonstrated that starting from 
the Hamiltonian level, the results cannot be reduced to 



where the typical definitions of the diproton channel ra- Eq. (|42| with (|43|. Thus, the diproton model is typically 



dius (Barker 2001) is 



R. 



dp 



ro(Ay^ + 2i/3), 



(43) 



used in an inconsistent way. Correctly applied, the dipro- 
ton model provides too small values of the width, see Fig. 



with ro = 1.45 (Barker 2001 1. Value 7dp is defined in 
an ordinary R-matrix form ( 33 1 . The function p{E) is 
the "density of diproton states" typically used in several 
forms: 



28 a), "diproton GF" curve. Moreover, Eq. (42 1 leads to 
a single, narrow, low-energy peak in the energy distri- 
bution in the "T" system, which differs from the exper- 



imental distributions, see Figs. 27 c), 31 |c), and 29 ^g). 
In addition, Eq. (42 1 provides an angle-independent de- 



p{E) = E6{E-Eo), 
p{E)^ Po{E,Rpp,l), 



P{E) 



ME)]- 



(44) 
(45) 
(46) 



Eq. (44 1 corresponds to a "fixed-energy diproton". It is 



evident that some energy should be contained in the in- 
ternal motion of the p-p subsystem. In the early works it 



cay probability in the "T" system in contradiction to the 
experimental findings. 



D. Three-body model of 2p radioactivity 

In the three-cluster model of 2p radioactivity intro- 



duced by Grigorenko et al. (2000 2002) the three-body 



49 



Schrodinger equation with the complex energy 



{H,-E + zr/2) = 



(47) 



is solved using hyperspherical harmonics (HH) method. 
The solution with outgoing boundary condition is found 
in this method as 



(+) 

3>t,JM 



{P, ^5 



P 



-5/2 



Y.K,^^^^K,{p)Ji:^m- (48) 



The hypermomentum x = \/2ME is the analogue of the 
ordinary momentum in the HH three-body approach and 
the hyperradius p is a collective variable describing the 
"breathing" motion of the system 



A1A2A3 



'■23 



•■31 



A1+A2 + A3 V A3 Ai A2 



(49) 



The hyperspherical harmonics J^K-y are functions of the 
5-dimensional "solid angle" fl^ = {6p,ilx,^y}. Here ft^ 
and riy are ordinary solid angles of the Jacobi vectors X 
and Y and tan(6'p) = ^/MJMyX/Y. 

For radioactivity problems with \Et\ ^ F the following 
procedure was found to be reliable. First, the discrete 
spectrum problem 



iH3 - £;t) *3b = 



(50) 



is considered with some box boundary conditions (e.g. 
zero or quasistationary) yielding the "box" WF ^'sb and 
the value of the real resonant energy Ex- Then the in- 
homogeneous equation 



{H3 - Et) ^^3^^ = -ir/2^. 



3b 



(51) 



is solved with arbitrary F and the actual width is defined 
afterwards as the flux through a hypersphere of large 
radius divided by the normalization within this radius: 



Im 



F= ^ = 
N 



P=Pn 



(+) 



(52) 

Formula (52) is intuitive. Still, some formal issues con- 



nected with its derivation can be found in Ref. (Grig- 



orenko and Zhukov 2007a) 



The asymptotic form of the three-body potentials in 
the hyperspherical harmonics method is 



(P) 



K 



OK'y,K'Y 



Mp ' 

(53) 

where the multiindex {K^} = {K, L, S,lx,ly, s^} is a 
complete set of quantum numbers for three clusters. The 
matrix UK-y.K'-y' arises due to contributions from the 
short-range nuclear forces, and Njcy^K'-y' ^ 3 are inte- 
gers: the effective contribution of the short-range forces 



decreases as p~^ or faster in hypersherical space. The 
diagonal centrifugal term Ck = ^k{jC-k + 1) depends 
on the "effective angular momentum" Ck ~ K + 3/2. 
Coulomb pairwise potentials generate the long-range 
part of the hyperspherical potentials behaving as p~^. 
From the technical side, these three-body hyperspher- 
ical "Coulomb like" potentials cause problems due to 
long-range channel coupling (nondiagonal "Sommerfeld 
parameters" rjx-y.K'Y) that does not allow one to de- 
couple the HH equations in the asymptotic region. To 
deal with this problem, the finite-size potential matrix 
(in truncated hyperspherical basis) can be diagonalized 
with respect to the long-range term by the orthogonal 
transformation V = A^VA: 



V, 



ip) 



p- 



Mp ' 



JK-/,K'i' 



(54) 

This potential includes nondiagonal "centrifugal" terms 
Ckj.k'j' and, to achieve the asymptotics in the diago- 
nalized representation, we still need to go very far in p, 
where the terms ^ p^^ become negligible compared to 
those with ^ p^^. At such values of p, the hyperradial 
part of the asymptotic solution with pure outgoing na- 
ture can be constructed using the Coulomb functions: 



K'-y' 



X [GcoiVK'Y,xp)+iFco{VK'Y^'^P)] -(55) 

These boundary conditions are exact in the truncated HH 
basis at a hypersphere of very large radius, depending 
on energy and basis size. For example, for ^^Fe with 
Et = 1.154 MeV and i^max — 20, radii pmax between 
500 and 2000 fm are needed to get reasonable solutions. 

The Coulomb three-body decays of nuclei have also 
been studied by other HH based theoretical methods. 
The S-matrix approach was used in calculations of ^Be 
g.s. in Refs. (Danilin and Zhukov 1993 Descouvemont 



et a?.[|2006[|Vasilevsky et al.\\200i\ and ^^O g.s. in Ref 



(Grigorenko et al. 2002). WKB calculations with adia- 



batic HH potentials were done for the width o f the first 
excited state of ^^Ne in (Garrido et al. 2008). Having 



comparable quality for internal structure, these methods 
may have deficiencies dealing specifically with radioactiv- 
ity problems: S-matrix calculations are difficult for very 
narrow states and WKB can not provide momentum dis- 
tributions. We note that other HH based methods have 
been used to calculate the momentum distributions from 
three-body decay of several excited sta tes in light nuclei 
dAlvarez-Rodn'guez et am2010[ [2008bl. 



1. Lifetimes 

The Coulomb and p-p potentials are the known ingredi- 
ents of the three-cluster calculations. In the light nuclei. 
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the core-p potentials are fixed by fitting to the exper- 
imental single-particle spectra of the core-p subsystem. 
In most cases this allows to obtain a reasonable descrip- 
tion of the energy and structure of the ground and low- 
est excited states. For heavier nuclei, where the single 
particle description is not well justified and spectra of 
subsystems are often not known experimentally, certain 
systematic guidelines were proposed by |Grigorenko et oT] 
(2003a I. Choosing the depths of potential components 



with different I values we can vary the structure of the 
2p emitter for a fixed proton decay energy E2r in the 
coie-p subsystem. A short-range (in the hyperradius) 
three-cluster potential is used to vary the decay energy 
Ex- 

For the majority of prospective 2p emitters the de- 
cay energies are not known. Therefore, the lifetime pre- 
dictions are provided in terms of possible lifetime bands 
bound from above and below by calculations with pure 
configurations. The calculations for the p-f shell nucleus 
^^Fe for different cases of configuration mixing are pro- 
vided in Fig. [28] It can be seen that within the three- 
body model the precise lifetime data can be used to ex- 
tract structural information about 2p emitters. 

A broader view of the true 2p decay phenomenon is 
provided in Fig. [35] For the light 2p emitters specific 
lifetimes are predicted (gray circles). For the heavier 
2p emitters the predicted lifetime "bands" are in a good 
agreement with the experimental data. The lifetime 
range for the known true 2p emitters spans about 18 
orders of the magnitude. This plot emphasizes the com- 
plexity of the problem requiring a variety of experimental 
methods to cover the possible lifetime range. 
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FIG. 35 (Color online) Lifetime vs. decay energy systematics 
for several known and prospective true 2p emitters calculated 
in three-body model. Hatching indicates lifetime ranges ac- 
cessible to diflerent experimental techniques. Experimental 
results are shown by diamonds. Gray circles show specific 
predictions, where available. 



2. Spatial correlations 

The evolution of spatial correlations in the three-body 
model is illustrated by the case of the ^^Fe structure with 
Wip^) = 98%, see Fig. 



36 



In the Jacobi "T" system, this 
case of almost pure [/^Jo" configuration has a very emi- 
nent correlation pattern with four peaks in the internal 
region. Such correlations can be related to the so-called 
Pauli focusing - counting of excitation quanta in the "T" 
system. 

It can be seen in Fig. [36] that under the barrier the 
four-hump structure is dissolved and a different correla- 
tion pattern is formed while transition to the asymptotic 
region takes place. Finally, the momentum distribution 
with one peak in energy between two protons is formed. 



see Fig. 37 ^b). This happens due to p-p interaction in the 
subbarrier region ( Grigorenko et al. 2000 ) . It is evident 



here that the penetration process infiuences the distri- 
bution strongly and the information about correlations 
in the nuclear interior is present in the final momentum 
distributions in a very indirect way. Detailed theoretical 
calculations of the subbarrier propagation process are re- 



quired. 



3. Momentum correlations and nuclear structure 



The correlation pictures calculated in the three-body 



model are provided in Figs. [27][a,b), )3T[a,b), and 29'e,f) 
for ^Be, ^^Mg and ^^Fe, respectively. Despite the strong 
differences there are certain common features in these 
distributions. 

(i) There are Coulomb suppression regions due to p-p 
repulsion at e ~ in the "T" system and at e ^ 0.5, 
cos{9k) — 1 in the "Y" system. Analogous regions due 
to core-p repulsion appear at e ~ 0.5, cos(9k) ^ ±1 in 
the "T" system and at e ~ 0, e ~ 1 in the "Y" system. 
These regions grow drastically as we move from light to 
heavier 2p emitters. 

(ii) In the "Y" system the particles are concentrated 
around e ^ 0.5. This effect becomes more pronounced as 
the core charge increases. This particular aspect of the 
correlation picture is consistent with the original predic- 
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FIG. 36 (Color online) Spatial correlations in the "^^Fe WF 
in the "T" system for W^(p^) = 2%. Panels (a) and (b) shows 
different radial ranges. Logarithmic scale: two (a) and one 
(b) contours per order of the magnitude. 



tions of Goldansky and with calculations performed in 
a simplified direct decay model, see Fig. [32] In the "T" 
system this effect appears as a concentration of particles 
in the T-like configuration with cos(0fc) ^ 0. 
(iii) The p-p final-state interaction gives rise to enhance- 
ments at small values of e in the "T" system. This en- 
hancement varies from minor to a very expressed one 
depending on the particular dynamics of the system, see 
also Fig. [37] 

An important question about 2p radioactivity asked 
from the earliest days of this research is "Which kind 
of information about internal structure can be extracted 
from the 2p decay data?" It was not clear in advance 
that information about the nuclear interior would sur- 
vive in the process of penetration of protons through the 
Coulomb barrier. The calculations in the three-body 
model have demonstrated that correlations in the p-p 
channel are sufficiently sensitive to nuclear structure to 
be discussed as prospective instrument for nuclear spec- 
troscopy in both s-d and p-f shell nuclei, see Fig. [37] 
These correlations are manifested by the e variable in 
the "T" system and the cos{9k) variable in the "Y" sys- 



tem (Figs. 27 and 29 1 



4. Long-range character of Coulomb interaction in the 
three-body continuum 

The formulation of approximate boundary conditions 
Eq. (55) is satisfactory for distances around 10'^ fm. Such 



distances appear to be insufficient to get converged mo- 
mentum distributions for heavy 2p emitters. Extrapola- 
tion by means of classical trajectories was used by |Grig- 
orenko et al. (20101 to improve the momentum distribu- 



tions. Classical extrapolation modifies mainly the cos(6'fc) 
distribution in "T" system and the e distribution in "Y" 
system, see Fig. [38{^a,b). The effect of classical extrapo- 
lation is important already on the current level of exper- 
imental precision. Fig. [38[c). A huge range is required 
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FIG. 37 (Color online) The energy distributions between two 
protons ("T" Jacobi system) in different s-d shell nuclei (a) 
and with different assumptions about the internal structure 
of "^Fe (b). W{\?-) is the weight of the \(^\ configuration in 
the nuclear interior. All distributions are normalized to unity 
maximum value. From (Grigorenko and Zhukov 2003b I. 



both for the extrapolation range (^ lO'^ fm) and for the 
starting point of the classical procedure (~^ 10"^ fm) under 
typical 2p decay conditions. The maximal extrapolation 
radius is comparable to atomic distances and sensitivity 
of 2p momentum distributions to atomic screening of nu- 
clear Coulomb potential can be found. This brings the 
2p decays to the borderline with atomic phenomena. 



E. Two-proton radioactivity and many-body nuclear structure 

A very important aspect of 2p radioactivity studies 
is the understanding of the connection between specific 
asymptotic observables (2p correlations) and the WF 
structure in the internal region (nuclear spectroscopy). 
In Sec. IVII.D.3l such a connection was demonstrated for 
the three-body model of true 2p decay. However, the 
discussed three-body approach should be regarded as a 
first approximation model of this process. By analogy to 
two-body decay, the exponential component of the width 
is expected to be treated appropriately in this method. 
The three-cluster model accounts well for single-particle 
degrees of freedom in the core-p channel. Therefore, 
calculation of the "preexponent" in three-cluster model 
is well justified only for the light closed-shell (closed- 
subshell) systems or systems with closed-shell (closed- 
subshell) core. With increasing mass number fewer sys- 
tems can be found along the dripline whose dynamics can 
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FIG. 38 Contour maps of the momentum density distribu- 
tion on the kinematical plane {e, cos(Sfe)} for *^Fe in the "T" 
Jacobi coordinate system without (a) and with (b) classical 
extrapolation. Panel (c) shows comparison with experimental 
data of ( [Miernik et al.\ |2007c[ | . Experimental and theoretical 
events are convoluted with the experimental resolution. From 
(Grigorenko et aH[20Tol). 



be reasonably described in the three-cluster approxima- 
tion. Taking many-body effects into account in the 2p 
decay calculations is a natural next step in the studies of 
the phenomenon. Here, we briefly sketch several possible 
approaches to this problem which are being considered 
presently. 



1. Spectroscopic information in R-matrix approaches 

The information about many body structure is embed- 
ded in the simplified models of two proton radioactivity 
via spectroscopic factors. The single particle proton fac- 
tors 9"^ Eq. (31 1 were used for this purpose within the 



direct decay model Eq. (37 


) (Barker 


1999 


2003 


Brown 


et al. 2002 


Kryger et al. 


19951. 



The diproton spectroscopic factors [for use within the 
diproton model Eq. ( [42| ] are obtained by projecting WFs 
of valence protons on the "diproton quantum numbers" 
for two protons (relative angular momentum zero and 
total spin zero). It is also typically assumed that the 
relative motion of the protons is represented by the lowest 
oscillator. Calculations of this type can be found in Refs. 



Barker 


2003 Brown 


1991 Brown et al. 


2002 



2. Three-body model plus RMF amplitudes 

A simple method of structure treatment in the three- 



cluster model was proposed by Fomichev et al. (20111. 
This approach is equivalent to a standard method used in 
the R-matrix description of two-body decays. The major 
components of the three-body cluster WFs with J'^ = 0+ 
can be written in a schematic spectroscopic notation as 



(+) 



X,, 



(56) 



The values which can be put in correspondence with the 



components of the three-body WF Eq. ( 56 1 , to take into 
account the many-body structure, are overlaps of the 
many-body WFs of the precursor-daughter pair multi- 
plied by a combinatorial term. They can be written in 
the same spectroscopic notation: 



2!(A-2)! 



1/2 



' ^ 1. 



(57) 



For WFs normalized to unity in the internal region the 
amplitudes of overlaps do not sum to unity in the general 
case in contrast with amplitudes of components: 



1 



2p 



(58) 



Therefore the three-body width as a function of weights 
for all considered single-particle configurations {Xf} — 
{Xf, . . . , Xj^} should be renormalized as: 



r({xf}) ^ N,j,r{{xf/N2j,}). 



(59) 



For systems which are not closed-shell/closed-subshell, 
the coefficient N2p is considerably different from unity 
already in the Hartree approximation. The higher com- 
plexity of the three-body decays is reffected in the fact 
that inclusion of structure information can not be re- 
duced to multiplication by some factor, but should pre- 
cede three-body calculations. 

An advantage of the method is the opportunity to 
use relatively standard input from many-body theory, 
Eq. (57 1. For example, 'Fomichev et al. (20111 calcu- 



lated these amplitudes in the relativistic mean field model 



(Litvinova et al. 2008). A disadvantage is the different 
treatment of the internal structure and the decay process. 
The same problem, however, exists also in the conven- 
tional two-body R-matrix phenomenology. The center 
of mass is treated appropriately in the three-body decay 
calculations, but only effectively in the many-body part. 



3. Shell model approaches 

Several shell-model approaches have been developed 
for studies of two-nucleon decays. 

The Gamow shell model is based on the utilization of 



the so-called Berggren basis (Michel et al. 2003|. This 
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is a basis composed of Gamow states - single-particle 



states with decay asymptotic 



For homogeneous 



Schrodinger equations such states should have complex 
energies. The eigenenergy of the shell-model Hamilto- 
nian diagonalized in such a basis is naturally complex 
as well, with the imaginary part interpreted as a half- 
width. Technicalities make application of this method 
problematic for extremely narrow states. So far, it has 



been applied to light neutron-rich systems ( Michel et al. 



2010a and Refs. therein). 



A possible reason could be that the diproton approxi- 
mation tend to underestimate width, see the discussion 
in Section [VlI.C.3| The direct decay lifetimes calculated 
for ^^Fe are 20-50 times larger than experimental one, see 
Fig.|28| 

A disadvantage of the shell-model methods is that no 
matching to the three-body boundary conditions was re- 
alized so far and calculations of momentum distributions 
for decay products are not possible. The necessity to 
incorporate three-body asymptotics in the SMEC calcu- 



The real-energy continuum shell model was employed 
for studies of narrow states decaying by two-nucleon 



emission by|Rotureau et al. (2005 20061 and Volya and posed in (Rotureau et al. 



Zelevinsky (^20061. In this approach the Hamiltonian is 



lations was emphasized in (Blank and Ploszajczak 2008 



Rotureau et al. 20061. The formalism for this was pro- 



split in two parts 



H ^ H 



2006|) using hyperspherical de- 
The ob- 



composition (|48|) and boundary conditions (|55 1 
stacle, which 



Rotureau et al. (2006) faced here is in the 



SQ 



QV 



H 



VQ 



vv- 



(60) 



Functional subspaces Q and V can be identified as "in- 
ternal region" with discretized (shell-model) basis (pi and 
"external region" with continuum basis (pE defined by 
projection operators: 



contact character of the residual interaction in SMEC. 
Delta- function interactions lead to collapse of three-body 
calculations. Use of finite-range residual interaction is 



technically complicated in SMEC (Rotureau et al. 2006 1. 



V= / dE\cpE){^E\ 



(61) 



Elimination of the subspace V can be done by the Green's 
function methods leading to effective (energy dependent 
and non-Hermitian) Hamiltonian % in Q 



T~Lqq{E) — Hqq + Hq-pG''p\e)H-pq 



(62) 



which allows straightforward application of powerful 
shell-model methods. An advantage of the continuum 
shell model is the use of extensive shell model expertise 
and applicability to corresponding broad range of nuclear 
systems. 



The continuum shell model method of (Volya and 



Zelevinsky 2006 1 was applied in extensive studies of the 



neutron-rich helium and oxygen isotopes providing good 
description of the known cases of 2n decay and predic- 
tions for several unknown. In this approach the residual 
interaction in subspace V is neglected. Therefore, in the 
sense of the three-body continuum dynamics this model 
is analogous to the direct decay model. Sec. |VII.C.1[ It 
can be shown that the expression for the 2n width from 



An important conceptual problem exists that should 
be common to all shell-model approaches. The asymp- 
totic interaction in the p-p channel (in the V subspace) 
should be the "vacuum" p-p potential. This should be rec- 
onciled in the internal region with the shell-model resid- 
ual interaction which has effective character, reflecting 
many-body dynamics of nuclear interior. This reconcil- 
iation should take place somewhere close to the nuclear 
surface. The observables for 2p decay are especially sen- 
sitive to this radial range, see Sec. |VII.D.2| making it an 
important issue. 



4. Microscopic cluster models 

Microscopic cluster models do not have the conceptual 
problem mentioned in the end of the previous section uti- 
lizing the same (quasi) realistic nucleon-nucleon interac- 
tions in all space. Microscopic resonating group methods 
(RGM) have been applied to the lightest true 2p emitter 
^Be. 



Eq. (35) 



The complex scaling method of (Csoto 19941 was the 
first example of such studies. In this method the bound- 
ary conditions are not explicitly formulated and therefore 
only resonance widths can be calculated. 

Algebraic version of RGM (Nesterov et al. 2010 



in the form 



( [Volya and Zelevinsky[ |2006[ ) can be approximated by |Vasilevsky et al.\ |2001[) for three-cluster systems use WF 



The shell model embedded into continuum (SMEC) 
was applied by Rotureau et al. (2006) to studies of 2p 



^{A) = A [*(^i)*(A2)*(A3)*3^(p, n^)] 



(63) 



radioactivity using either the diproton or direct decay 
Ansatz. The later is called "sequential 2p emission" hav- 
ing in mind "emission via tails of higher-lying states", see 
also the discussion of Sec. IVII.C.2l to avoid confusion of 
terminology. The calculated diproton lifetimes for "^^Fe, 
■^^Ni, and ^-^Zn are typically exceeding the experimental 
values by a few times, see, e.g.. Fig. 28 for ''^Fe results. 



where '^[Ai) are cluster WFs, A = A1+A2+AZ, and 
is three-body scattering WF in the "S-matrix represen- 
tation" (three-body plane wave plus three-body outgoing 
wave) . The derived S-matrix is diagonalized and the ma- 
trix elements with the most expressed resonance behavior 
("eigenphases") are used to extract resonance energy and 
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width: 



VIII. EMISSION OF NEUTRONS 



(f5{E)/dE^ 



= 0, 



r = 2 [d5{E)/dE]- 



\E=Et 



A hyperspherical decomposition was used for WF 'i'^^ 
as especially suitable for imposing long-range boundary 
conditions for true three-body decays. A somewhat sim- 



pler version of the boundary conditions than ( 55 1 was 



realized in (Vasilevsky et al. 20011 with diagonal terms 



only retained in the Coulomb term of Eq. ( 53 1 . The alge 



braic method employs two different hyperspherical basis 
sets for the nuclear interior and asymptotic providing a 
natural and controllable treatment of these two regions. 
The disadvantage is that the calculations already for ^Be 
seem to be not sufficiently converged at the present level 
of sophistication compared to fully converged pure three- 



body calculations ( Grigorenko et al. 2009a I 



A general challenge for microscopic cluster models is 
that realistic calculations can be performed only for the 
lightest systems with effective nucleon-nucleon interac- 
tions typically employed. However, modern trend is to 
move to heavier systems and realistic interactions. 



F. Three and four proton emission 

The beta-delayed three-proton decay of ^^Fe was re- 
ported by Miernik et al. (2007b I. Most likely this is a 



sequential decay and it has been discussed already in Sec 
IV Three-proton decay of highly excited states of ^^Ne 
populated in fragmentation reaction, via narrow inter- 
mediate states in ^^F and ^"'O, was observed by Mukha 



et al. ( 2009 1 . Theoretical prospects of studies of this de- 



cay branch are so far unclear. 

Following the analogy with the true 2p emission {Ip 
emission is energetically prohibited) we can define the 
true Ap emission (emission of i protons is energetically 
prohibited for all i < A). The only candidate to fulfil this 
condition studied so far is *C. It was shown recently that 
it is not a true Ap emitter (Charity et al. 20101. This 



nucleus undergoes a sequence of true 2p emissions decay- 
ing via the ^Be ground state. The next candidate for the 
true Ap emitter is 21Si but no information is available for 
this isotope yet. 



With the progress in reaching experimentally the neu- 
tron drip-line, the interest to study nuclei beyond this 
limit is rising. Emission of protons beyond the pro- 
ton drip-line has an analogue in neutron(s) emission be- 
yond the neutron drip- line. In this Section we want to 
comment on the possibility of neutron(s) emission which 
may take the form of neutron/few-neutron radioactivity. 
Some aspects of such processes were discussed by |Thoen-| 
(2004|. 



nessen 



For illustration we consider theoretical estimates for 
the two pairs: ^^S-^^0 and They are isobaric 

partners but not mirror nuclei (the {Z, N} values are 
{16, 10}-{8, 18} and {15, 10}-{8, 17}). They have similar 
mass/charge ratios, opened s-d shell, and should presum- 
ably have comparable structural and radial characteris- 
tics. ^^S is a candidate to be a true 2p emitter and ^^O 
could be a narrow two-neutron resonance. 

For the one-neutron emission, simple estimates can be 



obtained by the standard R-matrix expressions, Eqs. (31 1 



and ( 32 1 . The results for protons and neutrons are shown 
in Fig. |39Ka). Assuming a possible s- or d-wave ground 
state for A — 25, one can see that the decay energy win- 
dow corresponding to proton radioactivity ranges from 
~ 50 to ~ 200 keV. In contrast, for the neutron-emitting 
partner (^^O) to be classified as radioactive, the decay en- 
ergy, even for the d-wave, would have to be smaller than 
1 keV. It is highly improbable that such a fine-tuned en- 
ergy is actually found. It seems that a realistic chance to 
observe one-neutron radioactivity may appear only for /- 
wave and higher-Z states. As long as the heavier neutron 
drip-line nuclei are not known, long-lived neutron emit- 
ters cannot be excluded. However, it is likely to happen 
beyond the s-d shell. 

The two-neutron radioactivity, or true two-neutron 
(2n) decay, is defined in analogy to the true 2p decay (see 
Sec. VILA I as a simultaneous emission of two neutrons 



where one-neutron emission is energetically prohibited. 



The possible 2n radioactivity, shown in Fig. 39 b), has 
a few important differences in comparison to the one- 
neutron radioactivity. 

(i) Low-energy s-wave neutron emission could take place 
in the form of a virtual state, which can not be inter- 
preted in terms of width (the neutron s-wave curve is 
missing in Fig. |39ja)). For the 2n emission the phase 
space for a few-body final state creates additional effec- 
tive centrifugal barrier, so that a narrow resonance state 
is formed even for the decay of a [s^] configuration. The 
possibility to observe narrow "three-body virtual states" 



build on \s ] configurations was discussed by Grigorenko 



and Zhukov (20081 for the case of "He. The search for 



the 2n ground-state decay has indicated only the exis- 
tence of quite broad (T > 1 MeV) states in ^H, ^°He, and 



131 



Li (Golovkov et al. 2005 



Korsheninnikov et al. 19941 



2009 Johansson et al. 2010 
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FIG. 39 (Color online) Estimated widths for the neutron (a) and the two-neutron emission (b) compared with widths for the 
proton and the two-proton emission. Panel (c) shows estimates for the four-neutron emission. Hatched areas indicate lifetime 
ranges accessible by different experimental techniques, see also Fig. |35[ 



(ii) Similarly to the p-2p decay comparison, the widths 
for the true two-neutron emission are much smaller than 
for the one-neutron emission with the same decay energy. 
The energy window for 2p radioactivity of ^^S extends 
up to 500-1700 keV, thus is about an order of magni- 
tude larger than for the proton radioactivity of ^^P. The 
estimated relative increase of the energy window for 2n 
radioactivity compared to In radioactivity is even larger 
(about two orders of magnitude). For example, for the 
[cP] and the [p] configurations, the true 2n decays would 
be classified as radioactive for decay energies ranging up 
to ~ 200 keV and ^ 600 keV, respectively. Such broad 
ranges make the search for 2n radioactivity much more 
promissing. 

(iii) In contrast to the In situation, the 2n estimates in 
Fig. |39[ b) should be interpreted as lifetime limits due to 
the possibility of configuration mixing. The [s^] and [p^] 
curves are likely to provide lower lifetime limits for s-d 
and p-f configurations, respectively. The [(fi] and [p] 
curves provide upper lifetime limits for them. From the 
experience collected in the 2p decay studies we may argue 
that the realistic results are located between the limiting 
curves, but closer to that for the lower I. 

The systems decaying by a simultaneous emission of 
more than two neutrons should have longer lifetimes for 
the given decay energy. The estimates of the true four- 
neutron decay width can be performed in the approxima- 



tion generalizing the direct decay model of Eq. (29 1. For 



a true four-neutron decay this model gives (Grigorenko 



et al. 20111: 



l-£l-£2 



JQ 



62) Et- The and r„. are energy and width of the 
lowest resonance in the core and i-th neutron subsystem. 

The decay of is strongly suspected to proceed by 
a true four-neutron emission. The decay energy is un- 



certain, but seems to be smaller than 2 MeV (Caamano 



eram2008| [Korsheninnikov et a?. | [20031 [Nikolskii et al. 
2010 1. Th e ground-sta t e dec ay energi es of the subsys 
tems: "^H (iTilley et al.\ |1992| and (iGolovkov et al. 



20051 are around 2 MeV. There is only un upper limit 



of 1 ns for the lifetime of (Golovkov et al.\ 12004 1. 



The nearest heavier candidates for the true four-neutron 
emitters are ^'^Be and ^^O. 



In Fig. 39 c) the estimates for the true four-neutron 
emission lifetime are shown for ^H, ^^Be, and ^^O. For 
simplicity we used = 1.5 MeV in all cases. The other 
parameters are chosen according to guidelines from Sec. 
|VII.C.1[ The orbital configurations were chosen to yield 
results closer to the upper limit for the width, there- 
fore, the provided estimates are conservative. These es- 
timates indicate that long-living states decaying by true 
four-neutron emission are possible, with decay energies 
up to hundreds of keV, even for p and s-d shell nuclei. 

In summary, the observation of neutron radioactivity 
in s-d shell nuclei seems unrealistic. It is more probable 
that this process occurs in heavier {p-f shell) systems. 
The discovery of 2n or 4n radioactivity is much more 
probable since the energy windows corresponding to the 
radioactive timescale are much broader. However, both 
2n and 4n lifetimes depend strongly on the nuclear struc- 
ture and quantitative predictions require further investi- 
gations. 



=1 (£^»-^^«.)2 + r„^(£,)2/4 



Here Ei = EiEt for i = 1, 2, 3 and E4 — {1 — ei — 82 — 
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IX. CONCLUSIONS 

The exploration of the driplines of the nuclear chart 

have made great progress during the last decade. Our 
knowledge on nuclei along and beyond the proton dripline 
has increased significantly and studies of the neutron 
dripline has now reached the Ne-Mg region. New ra- 
dioactive decay modes appear here, namely one-proton 
and two-proton emission from long-lived states as well 
as many beta-delayed particle emission branches. This 
review has given an overview of the present state of the 
field. 

Beta-delayed particle emission will be prominent at the 
driplines and many modes will mainly occur there, the ex- 
ception being beta-delayed one- and multi-neutron emis- 
sion that is expected to play a major role for neutron-rich 
nuclei, in particular in the still unexplored part of the nu- 
clear chart. 

The single-proton emission, be it as radioactivity or 
following a beta-decay, is now well understood and estab- 
lished as a valuable probe of nuclear structure evolution, 
in close analogy to the role that alpha decay has played. 

Two-proton radioactivity is the most recently discov- 
ered decay mode. It is as well sensitive to nuclear struc- 
ture, but there are interesting physics qiiestions in the 
emission process itself that are not fully clarified and 
are intimately related to the more general problem of 
three-body break-up. To some extent this also holds 
for beta-delayed multi-particle decays, where only a few 
cases have been thoroughly explored experimentally. 

In all the cases mentioned theoretical and experimen- 
tal activities are ongoing in order to tackle the important 
questions that await clarification. An essential ingredient 
in answering these questions is the present rapid devel- 
opments in the production of secondary beams of near- 
dripline isotopes and the corresponding continuing evo- 
lution of experimental techniques. Along with dedicated 
theoretical efforts this should ensure that the progress 
continues. 
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